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Abstract 

r • x P« riment entitled the Middeck Active 

Centro! Experiment (MACE) proposed by the Space 
Engineering Research Center (SERC) at the 
Massachusetts Institute of Technology is desenbed. The 
th i eCtI Zi i° f thj r Program is to investigate and validate 

ftariS^tS?® ? f j he dyn f“ C J > ° f an acrive ’y controlled 
flexible, articulating, multibody platfonn free floating in 

zero gravity. A rationale and experimental approach for 
the program are presented. The rationale shows that on- 
orbit testing, coupled with ground testing and a strong 
analyticalprogram, is necessary in order to fully 
understand both how flexibility of the platform affecu 
;£* ^ * problem, as well as how gravity perturbs 

thus structural flexibility causing deviations between 1- 
and 0-gravity behavior. The experimental approach 
physics of multibody platforms, by 
identifying the appropriate attributes, tests, and 
performance metrics of the test article, and defines the 

framework 11 "^ t0 successfulI >' validate the analytical 
Introduction 

ecraft have become the foeus of intense 
scientific and engineering research in recent years. 
Beginning with Space Station Freedom, spacecraft will 
re8 t™t.d to the mass and volume limits 
imposed by a single launch vehicle. Instead, they can be 
assembled on-orbit, using astronauts on EVA or 

(L<?^nV qU1?ment ' T !' eS * Urge Space Structures 
, . will have increasingly stringent pointing or shape 

requirements, and a simultaneous reduction in their 
structural mass and associated stiffness. It will often be 
impossible to fully test such a structure on the ground 
before its operational deployment. Therefore, the 
question that must be dealt with is how can confidence in 
the designer, ability to predict the on^rb.t structural 
dynamic behavior of such spacecraft be increased in 
order to assure operational success of its mission? The 

effi^nr I«Q 168 dev,,0 P ment ^ an effective and 
efficient LSS qualification procedure. 

Kgure 1 shows the possible qualification tests that 
can tie performed on a spacecraft in order to characterize 
its on-orbit open and closed-loop structural dynamic 
behavior. Each cell in the matrix corresponds to a test or 
accurate analysis which the designer must be capable of 
performing before qualifying the operational vehicle 
The least expensive test, a scaled test of a spacecraft 
component performed in atmosphere on the ground, is in 
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corner ™J" rfo f lnad m “ in ^wer right 

sequence of temt^ j , d * V * op ^s most *fficient 
U« q A. ? u- , U ^° d anal >' m# *' this matrix, which 

the^peratioiud £££* *“ ** P * rfo ™“ C * 

The goal of the Middeck 0-gravity Dynamics 

tn P uZT nt i M0D Fl fa F' ly of “®ht i experiment 
R^iX r° n ? UCUd by th ? MIT Spac « Engineering 
rfrtT ’ '* t0 * XP u lor * tha gravity-dependent 
aspects of this matnx in order to develop the analytical 

and j* 9 * ^uences necessary to conduct an 
tv, j e^^ent spacecraft qualification procedure. 
Tins development is accomplished through a series of 

conHi!rr«^f “Fk* y inexpan9iv « flight •xperiments 
conducted m the interactive shirtsleeve environment of 
the Space Shuttle Middeck. 
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MODE -1, which is to ba flown in August 1991 will 

fl°tS^ a ks S T l 0f0 k i r n * l0 ? P dynamic «*P«riments on 
fluid tanks, deployable and erectable trusses, flexible 

appendages, and rotary joints typical of future 
spacecraft. MODE-1 is funded under the NASA OAET 
ln*SUp program. 

MODE-2, which is to be flown in September 1993 is 

u* °a thC M 5 DE ‘ 1 instr unientation with the new 
Middeck Active Control Experiment (MACE) test 
artic.e. It will continue the work begun in MODE-1 but 
will extend it to include structures utilizing Controlled 
Structures Technology (CST). CST uses active control 
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Overview 

The Space Engineering Research Center (SERC) at MIT, started in 
Tlllv 1988 hXromrieted two years of research. The Center is approaching 
thl operational phare of its first testbed is midway through the 
construction of a second testbed and is in the design phase of a third We 
presently have seven participating faculty, four participating staff 
members^ ten graduate students and numerous undergraduates. 

This report reviews the testbed programs, individual graduate 
research other SERC activities not funded by the Center, interaction vnth 
non- MIT organizations, and SERC milestones. Publis e papers m 
possible by SERC funding are included at the end of the report. 
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Testbed Program 


Optical iNimtr^KQMETER Tpyrnm 

cr _„. first of the MIT SERC testbeds is based on a design for a 35 meter 
IKSSijnjf ^ op . tlc f 1 interferometer, chosen for the stringent alignment and 
m^Sn Ura ^ l control requirements placed on such a structure byTcientific 

Str a hpdr«j eC f - f At th u .T of the la8t re P°rt. a scaled 3 f mete? 

tr , uss / ra , me had been assembled and suspended in thl 
aboratory, and optical component design was underway for an internal 
laser metrology system capable of measuring structural^salieniSent^to 
an accuracy of 10 nanometers. At the time of this™ ting XcSSction 
and integration of testbed optics, sensors, real time computer and support 
18 and preliminary tests have WpKfon 

Th t, reader ls JfI erred h) a full description of the state of the 
“ the Publication MIT s Interferometer Testbed (included at the 

theN^T* w/ff ’ P re P arad in August for inclusion in the proceedings of 
S!,w: ^ 0rk J h ?P °“ Technologies for Space Interferometry, which two 

the areas Srtty’ for ffe teft^ = f 8ummary “ Panted here to highlight 

to ^ ^easure^nter^f^at^e^i^h'changes 

SmSK?? 1 fle i ablhty - A 670 P Watt HP and beam steeriiiXJ 
?h?pf 1 ?/^?- ted a * trus « vertex and Powers optical measurement^es to 
three (cat s eye) retroreflectors at distant locations in the truss. Amfient 
laboratory disturbances measured by one optical leg are approximately 20 
“S™ broadband) which is 2-3 times below tL deffred XedJoop 

performa^crLlwc Pre88nCe ° f 8caled d,8turba nce8) dictated by ou? 

the tn^s^J^^^tro^on^ate^wn^ra^Acdve^piezoelectric^inounte are 

under construction for the three cat’s eye retroreflectors, which will be used 

tol ott^S4;(er 5 M^ W hTni 8 i^3LT^ 

_ S A 8 o tem I dent *fi ca tion: A multichannel spectrum analyzer and a 
nerihrm 32 acc f ler0 P* etei ;s has been integrated to the testbed and used to 
fndSlSc system identification of the "naked truss" frame The data 

occuring S in G c?umps. damPin ^ U% ^ l0W6r) and ^ htly spaced^ 

frpnnpn^^ EI S m ? nt . A . 228 node finite element model predicts a 

•W of the first mode that is in 10% error with results of the system ID 

, The “? dal 18 being updated, and is being utilized for controllaWlitv 
and observability studies used for sensor and actuator placement. V 

• fK Passive Damping: Several truss members have been constructed 
with viscoelastic material and will be installed in the truss tn ec * c ,ui: D u. 
nominal level of damping in the testbed ^e eimJs areTin^ Llted in » 
axial component tester constructed during the summer. S a 



Control Ex^ me “ t8; feedb8Ck 

have been dosed arouna acuve uicm 

During the next si^ months^ the mterferomet^^stbe gtu( j ents _ 

h P e addressed: 

* - i- nf fVio laflftr mel 


1) 

2 ) 

3 ) 

4 ) 

5 ) 

6 ) 


• tne iouowiug «auvu 

completion of all six legs of the laser metrology system 
integration and testing of active mirror mounts 
real time computer software 
integration of a scaled disturbance source 

zss «a t^rtisrsssssz of the 

completed truss, 
passive damping treatment 


MULTIBODY TESTBED 

Work is proceed^ on ^p f'^etailsf’see^he next^ection. This 
Control Expenment (MACE). ^o^deUuw.^ gtmctures that uti lize active 
experiment is designed to stu y d whose structural characteristics 

change S ^theCen^sSg 

sSuA? ^ ^ * mght Mper,ment on the 

M ' d structure: This Engineering Model ^^‘“^^“other^est article 
flexible Lexan, tubular ^‘"‘^'““"“rimbSfed payload will be located 
elements at five rigid, metallic nodi es. i (pitch and yaw but no 

at each end of the bus. The gm 9 ^ d body contro ] i s supplied by three 
roll about payload ^j n f* 0 ^ S l g wWls located at the center of the bus. The 
:Sbftursu» pended from three CSA zero-g suspens.on devrces. 
The test article is shown in the figure below. , o • 

Sensors: Baseline sensors located ^l^/i^ts^elch^f which can 
rate gyro packages and three tr V a * a £, c ® strain gauges will be bonded to 
be attached to anv node of ^VtAt article and each of 

g* s‘^W“p"^ be attached to the g.mhal 

m0t0r9 Types ofcontrol ^ Tesuonthe earing Porting 

fZ^f'acMuSoa . of ^ 

reacting against the jle^b i e s reference scanning profile. Multiple 

i'nteract^n a us^ a involve tl sim^ultaneo r us er p 0 inting and scanning of both 

Hardware status: All hard . w ^ the S excepUon a of^ th^actile'bending 
used on the MACE EM have arrived with the excepuo 
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segment, the gimbals, and the CSA suspension devices Tho 
segment will consist of a square Lexan rod with niP 7 nAlor*tm/< * e 

mounted on its edges. It is in the final design stages and will replace oneTof 
the paBswe segments and provide localized bendfng in the strarture The 
first gimbal is currently being assembled at LSMC and e 

dehvery in mid-December. LjfMC is pro“dtag sup^rt to Ae I^Ce1,v« 

tha n6 / two fiscal years in the development or the EM and control 
methodologies, and in the redesigns of the gimbal actuX-s The CSA 
suspense deuces are n the final stages of the manufacturer”' testing and 
are also expected for delivery in mid-December. 8Ung a d 

■ Instrumentation and signal conditioning equipment have been 

and are m t ^ e final stages of integration with the AC- 100 Real 
Time computer we have purchased arnf received The AC-100 in 
conjunction with the CSA suspension devices mil serve as a zero ? 
simulation and real-time digital control farilitv Q f+ QT . f*„ a ze .^?'f 
responsibility of supporting the MACE program .. 7 ft * initia ^ 

with SrSETliL ISlV 18 , COntinuing 6 } ud y contro1 issues involves 

£°b?' en rd V-S* been develop 

for their ^llys"s anVfn“St e Pr ° blem ^ i “ Ue ^ t0 various CSI 

aeveioped. a ruJly representative model is expected to be comiileted hv tb* 
that pofnTcTo^dr„r^tm^ h u e b^n ed comi>,etion ofthe testbed At 


y 



Multibody testbed test article 
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Research Programs 


This section gives a discussion of f lhe gradua^U research 

programs that are funded throug of ma j e during the last six 

“s' °£id eSd ^ 

?„ d ?„dicate tEe^otivation for * the research. The phrase general affiliation 
implies that the research is not tied to a specific testbed. 


DESIGN OF STRUCTURES FOR CONTROL 
Mr. Robert N. Jacques and Dr. David W. Miller 
(Interferometer Testbed) 

spacecraf^as^een^done 8 i^egl^^i^y?Jw^h|con^^des^n^^n^^f^e 

late in the design process. • JJZnrr^utrementrbecome tighter, many 
larger, and pointing and alignme ^ q r encie8 i n8 ide the needed 
flexible modes of the spacecra interaction of the structure with 

bandwidth of the controller. The fi t designing these subsystems 

the control is very strong and to to use a computer 

separately. Currently, the app PP : op i v (defined set of structural and 
program which iterates over a .P r CT ec ' S l y p ^ controller 

control design vanables Ce g. sizmg t ne or se veral objectives such 

gains) to find a design which The difficulty with this 

Is minimum structural mass or RMS moboa l^annc^ ^ . g yery 

approach is that even when an P . understanding is crucial in the 
little insight into why it is optimal. This ^underatanmng is prohibi t 

very early, stages of design when ^ y P°f^^m^tobe used in the 
formal optimization. Even when ? P . neP ded in selecting the general 
design oi the system, en jp 1 ?£ er \ I 'f - 1 Variables The work done here is an 
layout of the system and ^design optimization of these systems 

fnd”evek,p t tme Wes'of thSb" which can be used in the preliminary 
design of a controlled structure. 

These were simple spnng-mass-dashpo if/S Action of a controller and 
simple model is that it is analogous interaction ot a u 

' disturbance with a single mode o . . ’ c iosed form. This made it 

possible to derive the optimal contro l! performance of the system. 

ISS/S 

numerically optimal designs. 
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Some of the conclusions revealed by this approach have been quite 
interesting. The first concerns the role of damping in a controlled 
structure. It was found that in a perfectly modelled system, passive 
damping is very important when no control is used, however as the level of 
control is increased, the need for passive damping eventually vanishes 
This was because at higher levels of control, the damping tended to resist 
the control force instead of aiding it. However, when the control is designed 
tor a structure that has unmodelled dynamics, it was found that the 
damping was extremely important in the unmodelled modes, particularly 
at high gain. A second conclusion was that the most common way to 
optimize a structure was to make it less sensitive to disturbances as 
opposed to letting the disturbance into the structure and trying to control it 
°u ^ iyr en th ? disturbances were step loads, this implied that the structure 
should be made very stiff to reduce its displacement due to these loads and 
when the disturbances were impulses or white noise, the structure should 
be made to have high inertia where the disturbances are applied. Third it 
was found, that the design strategy one should use can depend on the 
amount of control effort available. Other methods for improving* the 
controlled performance, such as making the structure more sensitive to 
actuation were found to help only in very special cases. 

In the future, several areas still need to be investigated. The beam 
example used was very simple. It remains to be shown that the rules of 
thimab can be used on a more complex system. The SERC interferometer 
testbed is a good candidate for this analysis. Computer models of the 
structure have been developed. Also, it will be possible to perform 
experiments on the actual testbed to verify the results of computer analysis 
Another area that needs to be investigated is the effect of controller type on 
design .strategy. So far, the controllers used in this research have been 
optimal LQR/LQG. These controllers work very well in perfectly modelled 
systems, but have been shown to have some serious short comings in the 
way of robustness. Actual MIMO controllers used on space structures will 

^vi -l .y bo based on H<x> methods. Their influence on structural design 
will be important. 6 


A Statistical Modelling Approach for Broadband 
Control of uncertain Structures 

Mr. Douglas G. MacMartin and Prof. Steven R. Hall 

(Interferometer Testbed) 


Background 

The goal of this research is to develop a methodology for designing 
optimal active control systems for broadband control of uncertain flexible 
structures. Thus a control system is required which will give good 
performance over many modes, while being robust to parametric 
uncertainties in these modes. Rather than modelling the detailed modal 
behaviour of the structure using state space based methods, it may be more 
useful to model only some statistical aspects of the response. One approach 
that uses this philosophy is Statistical Energy Analysis, and the modelling 
assumptions from SEA will also be used in this research. It is assumed 
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that the individual moded ^ 

modal en ^®f:^j e r? g t i^se energies, and the energy can be related to the 
rierTnp^tlTt“patto 'usinf conservation of energy. 

In order to apply this, the poweH dowf bedevebped 

ttSSSL 

local dynamicsof ^ e f a ^ c ‘,“ re e d loopVnamics via a transfer function H, 

a charactenz^ion of the closed l^ p yn ^ reflccted int o the structure, and 

defined such that H H is tnereia p transfer function H is 

ait rSBPX«^ i sX3^t%!S 

wave of unit power. . • ctatia+iral 

Eneigy* Analyst 

true performance metric of interest is 

J = i:.C(a»fr{[(/ * ff ]}n„toWo) 

bS M<rn*rmeaSMes°tL a3 powe^in^tVom ^ma" 

each frequency band. 

Cost Functional Minimization , 

This cost functional combines aspects i of “ ^p?oach for 

Much of the research cost functional. Preliminary 

SSSSeMl Xch is included at the back of this report. 

One approach to optimizing the cost ^ctjorml ^ase ^ ° n Th e cost 

the necessary conditions using a Ri^tfand^i single Lyapunov equation, 
can be evaluated using a si g , Riccati and Lyapunov equations 

For a given, fixed .compensator order the Riccati ana Lagra ng e 

required to evaluate the cost ca Pp au _ ente( j cost set to zero. The 

multipliers, and the first van , . . ■ £_ om this approach, but they are 

necessary conditions are easily ob LQG problem, the 

difficult to Simplify and understand^ In greatly by 
conditions for a fixed-order compens ^ ig d^cuit to identify for 
introducing a projection operator [ ]. ne ^ e ssary conditions, resulting in 

^uati® th. Lyapunov equation that is used to evaluate 

the C ° There are a number of s 

problems, and the solution o *. problems were investigated; 

current optimization P r “J’f m r ; x ed h orde r approximation, estimation and 
^rtrof^obl^r'o^^tU'ensatoftor ^t, motor, or approximate, is 
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desired, to be applied to a finite set of models. This approach can be annlied 
to m y e robustness to parametric uncertainty, by choosing a finite number of 
different models to represent the range of parameter variation It is also 
applicable to certam fault tolerance applications, where there may actually 
be a finite set of models which are possible. The multi-model 
approxmatmn and estimation problems have been solved by identifying ^ 

ofir tk projectl ? n °Pf rator * and in both of these, there is a mSfmiSJ 
order for the approximant or estimator that is optimal [21 (also included in 
the back of this report). For the optimal control problem, there is no a nrlori 
the order of a compensator which is optimal, and this is directlv 
related to the difficulty in analyzing the necessary conditions This result is 
also (apparently) true for our problem. inis result is 

Another related problem is that of finding a stable compensator that 
optimizes an #2 performance. (Note that the optimal compensator for a 
given system might not be stable.) This problem can be formulated so that it 
has the same double A c “ structure asthe mul^moder?on^l problem 
and the cost functional in Equation (2). A variety of researchers havp 
investigated this problem in the past, without obtaininl a particularlv 

finVtp^n^ 17 a g + °- nt i i 0ne c l oncluslon of this body of research fs that for a 
finite order, rational plant, the optimal compensator may be irrational A 

order.) •' &. compeZt^rt ot 

plant, and with almost the same cost as the optimal infinitP „J£! 

compe n sa tor The rationale for this belief co^ 'from numricll 

RntW>f Q tl0nS i f that j a i ve b t e , n carried out using a quasi-Newton algorithm 

l^L^ U - m0del P r T obl , em - and th e XyXoo cost functional have been 

2 ti 2 S “ anner - In both cases, and for several different plantothe 

wthiifa few'npr^nt C pfV? e £ Sat ° r ,i degree equal t0 the P lant ^ ad a cost 

within a few percent of higher order compensators. With a conmensator of 

possible ° r5er ' tHere W3S VirtUally n ° improvemeT^ the cost 

? ay a,8 ° be possible to minimize the cost by representing it in 
« n °‘k* r . fo ™ “ore amenable to different techniques. Along tlhf line 

wo^ ofzhou h p7oZ b m n w n h 0t€d ^ etWee 2 th i 8 ***- C08t functional, and the 
worK ot Zhou et al. [4], where the combined problem was studied from 

mdU u C6d norm P° int of view - The cost in [4J is the induced 
^“i 0f a P fa ^ with two inputs. If the second input is restricted to be 

Bernstein^ al iPtl* equ, f val f ent the ^ 2 /Z problem studied b£ 
oernstein et al. [5J. If the restriction is removed, then the Droblem is 

equivalent to our XuH* problem. Zhou et al. have not solved this case but 

this representation of the problem may lead to solution techniques 

The cost in Equation (2) also has a time domain formulation in terms 
of a stochastic Stackelberg non-zero sum differential game The control u 
minimizes the two norm of z under the influence of a white noiseTnput co 
and a deterministic but unknown worst case noise cuj. Once the control 

a a mp!?-fi h 2 Sen, * the A deteri ^ inistic noise uses fhis information and optimizes 
co f* A set of equations identical to those obtained Tom the 
instead 6 multipher a PProach can be obtained using this framework 
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COST FUNCTIONAL INTERPRETATION 

— « * * _ u » « A J 1 1 A 


The dereverberation procedure «m be integrated ^a^ra^ng all 

possible plants over some class of > a the average plant is the 

S-JS3SSS3 

iTolesire^en™ fts^" but j5t all the way to the real axis. This is 
the 

Krbe%h“cht a si5 a control law , and 

Satt^^rSeqrency is lifted in the Action which increases 
the cost. 

Future Work 

Ultimately, an experimental t^t of ^ 

desired. This will be done £ the MI T ^ d ^^?n!xmnining the use 
^“trec^a. !ferenu\ h thl a tss y forTa, control. The first 
experiment will be to perform rate feedback. , . . i 

An understanding of the optimal solution^of The 

control, and stable _ ^2 ° r P ^ to be improved in order to be a 

SKS^ ^ : 

allow m. 

SSS to'tim^d^lays^act^ato^an^sen^r d^am.cfand so 
forth. 
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Systems Engineering Approaches to the Design of cst 

Spacecraft 

Prof. Joseph F. Shea 

Srim 3°.?,. P T tjng r f equ i, remen ts Places constrain7s on the spa™ craft 

tools^sed for disturbance 8 T 

spacecraft " 8 thC systcms engineering approaches to the design of°CST 


A Systems engineering approach to Disturbance 
Minimization for Spacecraft Utilizing CST 

Mr. Christopher E Eyerman and Prof. Joseph F. Shea 
(Interferometer Testbed) 

bnng the system response to disturbances wit^rperfo^nce hZts * 

vehicle design of a Space-based Optical interferometer 

Mr. Andrew M. Nisbet and Prof. Joseph F. Shea 
(Interferometer Testbed) 

As a follow-on to the previously described research 
complement the development of the SERC interferometer testbed aliew 
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•nmiprt has been started to complete a full design of an interferometric 
?narPCTaft accepting only the tetrahedral truss structure as a given. Work 
tn has been in imderstanding the basic science requirements of the 

to da methods of disturBance min imization available. The next 

St b t e he to C «« 

t^e ^sTlega o^sp^eCTaSt^ttfti^e^a^ustment 8 ^ ^sed^ to producnthe 

science time and their cost in terms of mass, power, and equipment. The 

of th« tettg 

remaining spacecraft subsystems will come after this initial trade. 


ROBUST MULTIVARIABLE CONTROL STUDIES 
Prof. Michael Athans 
( Interferometer Testbed) 

This research deals with the development of concepts, theories, 

carrying out the research described in the sequel. 
graduate 6 8tudents 8l MOearch ^assistants' 

working toward their MS thesis. 

Mr Lublin continued his research related to the deveiopment of 
nominal models, and quantification of the 

isjsssr. is^-assas SjS MffS . wa 

finite-element model of the interferometer by explicitly modeling the twei e 

strut^The^'omplexity 1 ' fo^he^ res^tant^^rfite-element problpErequred the 

simpler 8 b^am^modef "de^Foped* (iuring^e^m^m^ 6 eporthfg^period, at 
least for the low frequency modes. 

We seek models of reasonable accuracy for initiating different 
multivariable control system designs; however, of equal 
Elly understand the nature of, and bound on, mwleling jsrrer.. These 
errors are both parametric due to uncertainties in damping, stillness, 
leometrv etc as well as dynamic in terms of poorly known high-frequency 
modes and other parasitic dynamic phenomena. However, one o y 

fi^w durine ^ this reporting period was that significant modeling errors 
ca^afso Sn tL'Toddfng" of the actuators and sens o^wtaA wg 
impart the nominal values and directions of the multivar 
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traiismissi 011 zeros of the transfer function matrix used to model the 
interferometer structure in control synthesis studies. 

Errors associated with the location and directional properties of 
transmission zeros may have serious impact upon the performance of the 
J 0 " ™ , system for disturbance rejection, especially if these are lightly 
damped zeros, forang a low bandwidth control system implementation 7 

fi^fr?n 188Ue + ^ 6re 1 / lve . 8t g j^d by Mr. Lublin, who compared the transfer 
function matrices of a two-dimensional truss structure with that obtained 
^a^app^mation and the models obtained using fiSte 
element methods. Although the modes of both systems were very close d 
difference between point force and moment controls vs then- physical 
distributed implementation resulted in a very significant modeling P eiror fn 
the transmission zeros and their directions. It is becoming more and more 
apparent that finite-element models may not be sufficiently accurate to 
implement a high bandwidth control system. Also, modeling ^or?fn both 
poles and zeroes may require larger damping than one may have initially 
suspected. Thus, these models must be refined by empirical data. For these 

interferometer^ testbed 6r identification of the 

_ .a P 0u 6^ a8 devoted most of his time in researching design 
methodologies for stability-robustness and performance-robustness of 

SSna y8t ^ 8 S1 ^ ifica ? t Parametric uncertainty. During this 

period we were able to derive some new results that can be used to 
design multivanable control systems using full state feedback strategies 

*ed fatelv as „ e8 rnch° r a \ 8 V mple " a «-»Pring aysum, which £sT een 
sea lately as a benchmark for assessing performance robustness for 

diverse control synthesis methods (Bernstein- Wei in 1990 ACC). We found 

nreLn^ f n^Wfi 0n \ esults ' n remarkable performance robustness in the 
° f ® 1 . gnificant uncertainty in the potential energy stored by the 

S™ g /T£ ertam spn r ng con J 8t | n J t) - state variable feedback gains are 
computed by means of a modified algebraic Riccati equation in which the 
nature and size of the parameter uncertainties are manifested in two extra 
terms: one that looks lie the term that appears in the formation 6 and 
seems to explain the fact that parameter errors can be partially modeled as 

w the n 1 °™ inal design; the other term modifies the 
8tate ,Y anabl ® penalties, and tries to force the system to move 
o?m!^ 8 °P aratin £ regimes that minimize the impact of uncertain energy 
stored m the system. These encouraging results have motivated us S 

thflt n thp° Ur theor . etlcal ^ es ® a 1 rch along these directions. It should be noted 
toterferemete? 86111 methodol °^ cannot be directly applied to the 

testbed, because we do not measure enough state variables We must 
extend the currently available results in a very non-triviaTway ^foreThe^ 

interferometer ^X d Ct,0n wlth the Posseted sensors and actuators in the 

• t . w ® bava also continued our theoritical effort that may provide us 
IwliSP 110 * Stance regarding performance limitations, say in terms of 
lnce *^fJ ectl j )n ’ for , dynamic systems characterized by lightly- 
damped stable poles and minimum phase zeros. Such performance 
are und f rs _ tood for unstable and/or nonminimum phase 

fntelSls U e S tc ? Fo e r S thPt^ 0m a / vanc . ed complex variable theory (pfisson 
integrals, etc.). For the types of structures (e.g. the interferometer testbed) 
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, .. fKp absence of modeling errors and other nonlinear 

under consideration, in the s a systems are “stably invertible so in 

actuator/sensor effects? 1 ^ht^Lcellent disturbance-rejection performance 
principle it is possible to j However, 

Ling dynamic multa^nab^ wm^ensators^WoVer^n^neaj^tigs, reduced 

from a pragmatic point decentralized implementation will impose a 
complexity controllers and decent alizea ? ormance while maintaining 

limit on achievable disturbance rejec p (hopefully) will make these 
stability. We seek a theoretical em)r desi 

tradeoffs transparent, and .reduc< e currently execute. Although the 

iterations that control ^sign ^g insight gained by von Flotow and 

theory is still at its infancy, it -confirm* i the insignt gainea ^ minimum 

hi, indents on singU-,nput «£«»»• “^SlbiutV if there exist, 

pole-zero pair, bear the .magmary 

axis. 


Mr. Gary Blackwood and Prof. A. von Flotow 
(Interferometer Testbed) 

The optical interferometer testbed has s t ^ e ™Ls P of?he Entire 

components that have . s “ a , C ^LLtain the relative pathlength 

structure. The control task is to ^a^erance of SO nanometers rms 
differences between these co “P an , e .^® ce9 (primarily a reaction wheel 
in the presence of scaled ujfh b?iadbanrand ( Larrowband components), 
disturbance model, with both b r ° aa ° a ?~ b ed opt ical systems under 
The testbed is representative of a ulf ^elativif ^leCTee^ of^passive and active 
consideration by astronomers, obiectives^ is still an open question, 

control necessary to ™ eet . pe r fo ^^ rat ^j J ; n MS thesis the feasibility of 
In related work. Jim Garcia de®°" 8trB “ d mounted to a flexible 

the modal 

damping ratios. 

mouiiU^l^be u^eTto'po^ti^nthe °mulU- 

order control systems that are d gn gn analytical framework, 

flexible mode, of the testbed 8t ™ ct “ r > e „.^' t " i " 1 , , be^nvestigated for 
dynamic modification, in performance 

their effects in reducing «!tudied both for the optical output 

improvement. These approaches mil be ' foL the basis 

location as well as for fbstmbanre '^“.V ^eaent time, three three-axis 

lm^^ntfto a ?heT:tV e d n ^rcur over the next few mfnths, and open 
loop input-output testing will begin. 


1 A 
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Realization of Impedance Matching for Uncertain 
Structures Using Passive Elements 

Mr. Ron Spangler and Prof. Steven Hall 
(Interferometer Testbed) 

The idea on which this research is based is that it should ho nnociKi 
realizable)*, gj, 

connection between the frequencv-domain ^ I f nuaa P r . ovide8 » 

system and a etate-spice-dolE ^?ep?SlktaUon BrieHv ", p08, * ive rea ! 
system has a transfer function matrix GKs) with the follomng properties? 8 ' 

- G(s) has no poles in Re[s] > 0 

- GT*(s) + G(s) is positive semidefinite hermitian in Re[s]>0 

realizaTion rfoSPcSfit (A B CD?‘JhSJI.‘^ t . fi)r ?" y «tate-space 

symmetric positive deHmte such ?haf e ” 8 ‘ matnCes P ' L ' and W P 

PA + ATP = - LL t 
PB = C T - LW 


WTW = D + D T (l) 

closed-loop system, or impedance compensator, to be positive real! Sd thuf 
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'Tbkm^s^tTere aK ^number^fwlSM^P^rW tTttaelfuatkins 

compensator may not be the best optnnum. 

. In rtaflemma A e n q tToLu^ed h fn ‘the 

XT^tor^ 

resistors and capacthjrs (or dasbpoU and *P™§^ 1 £3 8 g ome knowledge 
ground. I seems possible that the positive real which WQuld gerve 

of network synthesis would be to replace the ground with a 

tdr^ngToSr^s W rro^ranS-'e^rcb P .„ this area is 
expected to continue in the next six months. 

More recently however, another approach ‘“^^^^atTon 
passive controllers has keen formu . matrix necessary 

approach uses a cost-functional approach stiffness 

“” h ,n 

either case tie masses in the network are given a fixed vaue. 

Q to approximate total energy of vibration. 



passive 

network 


Figure 1: Teat Case for numerical Optimization of passive controllers. 

The plant model is that of an undamped beam in modal 

beam tip position an<f v velocity ' and l its i output is a matrices of 
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includes both the open loop controller and plant states, 
input has closed-loop influence matrix Bq 


The disturbance 


. ,T° £- nd A® necessary conditions for optimalitv of the 

controller K and C matrices write the closed-loop cost as y 01 the 

J = tr(PV) 

PA + A T P + Q = 0 (2) 

where A is the closed-loop system matrix and V = E{d 2 )B c B c ' T Adioinin? 
nHt2 pUn ° V equation to the cost in (2) via Lagrange multiplier nStri^^f 
and taking matrix partial derivatives with respect to P H K and C vields 
the necessary conditions. The A matrix is a 4x4 block matrix ’with C S,dJr 

dir?™?™ ° f -ll 8 blocks ’ 80 the necessary conditions resulting from the 
derivatives with respect to these matrices are in terms of cowespondin* 
blocks of H and P, and are hideously complex. Thus there is no closed-form 
solution for the optimal K and C matrices Iorm 

So finding the optimal values is a numerical optimization Droblem 
requiring the gradient of the cost with respect to each of the network 

for pI?ame e te?p SPnng damperX U Can be shown that thi ® gradient is, 


= 2rr \hP — ] 

dp [ dp J 


(3) 


multiplier referred* \o\bove^sofves° V eqUaU ° n in (2> and »• the L “^e 

AH + HA^ + V = 0 (4) 

To find the gradient of A with respect to p reauires SDerifvintr the 
topology of the controller network, so that the control?er C and K matrices 
as functions of the individual parameters are known. 

• f . ^ bus ^ ar t ™,° cases of ^e above controllers have been investigated 
with one mass and two masses, each of value m. The single mass case has 
only two parameters to be optimized (a single spring and damper) and the 
two mass case has six, as seen below. All optimizations were performed for 
a value of m/beam mass - 0 1. Using a matlab f^ction wWdf pe^ms a 

Ken perfomed Slate p i rameter8 8everal “Pt^tfonsTave 

^lThe g » tlVe ’ t hu8 / insu f i 1 ng . Passive definite K and C. In most Sses there 

Jrfwhfrh ^fh ber (p l SS i b1 ^ mfir \i te) of lo cal minima of the cost, any number 
of which can be reached given the proper initial parameter set. numDer 
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Figure 2: The most general 2 mass passive coinpenMtor«(riicture. 

velodty!^h^8 the^es^t^ft^ t^timlfir^ 

a ten-mode beam tothe classical vibration 

conditio^fo^optimafity^ere differs from the classical problem’s minimum 

of maximum response (really an ?&o criterion). , 

A two mass controller has also been optimi ^j^^Ce teen 

ten mode beam model. In the t which one or more of the 

found which are on the > consi trainti has ^ s0 been f 0U nd using 

parameter values are zero. A loc ^ “ f parameters (k 2 and c 2 ) have 
different initial conditions in which two .of ^the para 1 , d> \ t ifl believe d 

no effect, remaining at wha t^e ver lnx s i x -dimensional parameter 

there is a better minimum in the ^ ^ sharp, and will be 

space which has not yet been found It is probablyvery snarp „ d 

reached from a relatively ' small set of In the ten 
schemes for identifying this op 1 constraints, with k 2 - 0. 

^ Sw3£S^2wSS«SS^ 

strongly by the first mode res ?°.“”' • others vields the optimum. In 

job possible on this mode ana ignon & _ , , two-mass damper are not 
this case, all degrees of frredom afforded by a ‘^ d ma |^““^ ngle mas8 
needed, rather it is the increased massthaUs desire^ lurn^g^ ^ mu]u . 

m^rmXl^aTshoVn that the ' closed-loop ’ cost is still much higher than 
the more general optima discussed above. 

cost 

resulting passive compensators. The ““Xsed loop system matrix, but 
gained by minimizing the 9U norm ot the iciosea mop j numerical 

Iradients of this norm are , d '^ c n f ‘‘“^^.ngcamedoSt currently, 
procedure. There is a P*®* d ® ^ control in which a performance 

criterion i^ndnimized ^hiTeriL^rou^y^tisfying an ^ norm-bound 
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constraint. Many of these problems become the problem in the limit of a 
scalar parameter in the cost functional, but do possess calculable 
gradients for numerical optimization purposes. 
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Figure 3: Disturbance to tip velocity transfer functions of the undamped 
beam (solid) and the passively-compensated beam using a single mass. 


HH 1 - 

10 * 



Jstar=9604.3, k2star=0 (may no£be global opt) 

Figure 5: Transfer functions of the undam ped beam (dashed) and the 2 
mass compensator with the best optimal cost yet found. 




POINTING CONTROL FOR PAYLOADS USING VARIOUS 

measurement Topologies 

Mr. Mark E. Campbell and Prof. Edward Crawley 
(Multibody Testbed) 


Objective 


The objective of this research is to quantify, through analysis 

e r?u nmen ^ improvements made by successive relaxation 

of tho MiH S ffl?i n A S +° n ^P 0 . 1 ^ and successive additions to the hardware 

of the Middeck Active Control Experiment(MACE) model. naraware 


Approach 


The approach will be as follows: 


1 . 

2 . 


3 . 


Develop a simple two dimensional model of MACE and 
compare it to typical section models. 

Develop control schemes for the model by working through a 

C T b £ mg tbe relaxation of the constraints on the 
topolop and adding hardware to the system. Some of the 

finite clement “** — 

After selecting a few of the elements of the matrix as the best 

t£ nt i r< K a PPr° ac h es > these will be implemented and tested in 
the lab on the actual MACE test article. 


Work to Date 

planar M modefof C MAci an FW ‘lif i " vestl « atio " »f the dynamics of the 

v lb cat,IrmX e of D 4 F b U s th^^l 

-Ki as t .? ul ^f ted do j vn to a three D.O.F. model: the rotation of the bi!l 

problem. * ™ b f the nominal 8tate was the MACE sample 

mnri a iTflu n * ew m ? del ^ as then compared to the Non-CG mount flexihl* 
model of the typical section models written by David Miller in Linear Closed 
Loop Analysis, version 1.0. This showed the two models to £ ve^ siiSlar 
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The immediate future of the research is to begin a linear control 
analysis of the three D.O.F. model. Classical sequential loop closinres will 
be investigated initially and some MIMO control techniques will be looked 
• * ?! The nlan is to begin taking the linear model through the 

ESS MAC? sample* problem. The -suits of tins 

simplified system could then be used for comparison to the more 
complicated systems, such as the finite element model. 

The order or size of most of the nonlinear terms in the six D.O.F. 
model ^ll isoVe investigated, just to know what is being thrown away 
when the nonlinear system is linearized. 

Planned Results 

The planned results are that we will fully understand the advantages 
and disadvantages of the matrix control options of the planar model and 

^dfbet^ & m'jf the P?™del 

between ‘S'mojrcompUcated models’o" MACE* and the simple planar 
model. 


CLOSED-LOOP stability of multibody systems with 

DISTRIBUTED FLEXIBILITY 
Mr. Carlos E. Padilla and Prof. Andreas H. von Flotow 
(Multibody Testbed) 

In the past few months we have been investigating £ he .PF®^?* 
closed-loop stability of multibody systems with distributed flcxibilityj^g., 
Xtic man'"ulators, MACE test article). Using exmtmg ! resvdto for the 
control of rigid manipulators as a starting point, I am trying to obtain exact 
stabihty^ rSSts fo? systems whose models are known exactly and for which 
ah the states are available for feedback. (Now that the concept of state 
feedback for a distribuWd parameter system o ee d s to ^ e clanfie*defined by 
thiq research) It can be shown that for independent joint P1D control or 
flexible manipulators stability is guaranteed. 

considerations independently of the specific model for the distributed 
parameter system. 

Recent literature results attempt to show exponential stability of the 
tracki^ contrJl for multi-joint flexible-link manipulators A 
generalization of the computed torques approach for flexible mampulators, 
together with the Positive-Real Lemma, are used to show that stebility of the 
closed loop system is guaranteed when a class of positive real (passi ) 
controllers* are used as joint controllers 1 . However^ this residt depends 
both the joint controller and the structure being sufficiently stiff and havi g 
enough damping. Further, in order to take into account structural 
damping to increase the implied performance guarantees, the observability 
of the inodes used to model the flexible members from the joint sensor must 
be con si d ere d . This indicates that both stability and performance may 
depend on the choice of mode shapes: a modelling choice I 
the one hand, to generalize this result to the case when full state feedback is 
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available for control. On the other hand, I would like to investigate further 
the modelling issues affecting the stability proof and if possible extend the 
results to practical cases (i.e, non-vanishing regions of convergence). 

* 88 H es have become evident in my stability analvsis of thp 
MACE test article for ”not-quite-independent" joint PD control^ Stability 
results depend on the model and in particular stability bounds seem to 
depend on the particular choice of mode shapes used to discretize the 
distributed flexibility I am at present looking at these and other 
modelling for control design issues and my efforts will be directed 
towards completing a paper to appear in the 32nd SDM conference: Further 
Approximations in Flexible Multibody Dynamics. Specifically, I want to 
investigate how far a given set of dynamic equations, for the type of systems 
under consideration, can be simplified given a set of operational constraints 
m terms of maximum rigid body displacements and rates. It is foreseen 
that the degree of simplification will be tied again to the modelling choice of 
mode shapes. This suggests that some mode shapes capture the "essence" 
of the physical system better than others. I am interested in finding 
mndamentm conditions (e.g., the satisfaction of certain energy integrals) 
^ the b , est choice of mode shapes given that only a finite number 

oi them will be used to model the distributed parameter system. 
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Gravity-Suspension System Effects 

Mr. Daniel A. Rey and Prof. Harold Alexander 
(Multibody Testbed) 

Significant progress has been made since the last report towards a 

S rehensive understanding of the differential stiffening effects of axial 
ng and shear loads on a beam in the presence of axial, bending or 
snear deformations. It was found that axial loads are the only load type 
which contribute homogenepus stiffening while all other load types are 
coupled to off-diagonal elastic degrees of freedom. An example of such a 
coupled stiffening effect to the suspended MACE test article is the increased 
horizontal plane bending stiffness which arises when a vertical bending 
moment is applied to the torsionally deformed test article as a payload out of 

inJ!^f C ?i! *1“? 8 if w , 8 t0 5 rard ® the longitudinal axis and away from the 
center of the test article. Contributions to the shear mode of deformation 
have yet to be determined. This analytical understanding will be 
complemented by future research which will investigate the extent to which 
existing finite element software packages capture these effects. A 
parametric variation analysis will then follow using analytical and 
numencal modelling techniques to determine how differential stiffening 
affects the system eigenstructure; tracking also the slope and curvature 
node positions for the first three flexible modes. 
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a • of-. Jw nf thp first order effects of a suspension system on a 

w Irtfcle hae identified the different static and dynamic 
suspended test article nas iaenu ^ tem a i te rs the test article 

mechanisms by which ^ uspe ^^^ eac h attachment point due to 


COMMAND SHAPING FOR THE MULTIBODY TESTBED 
Mr. James Hyde and Prof. Warren Seering 
(Multibody Testbed) 

Work in the last six months has concentrated on: 

1) Solidifying a M.S. thesis topic, 

2) Completing the MACE inertial wheels, 

3) Performing and presenting initial modeling and simulation 
results, 

4) Identifying areas for future work. 

The theeie title has narrowed to: "Using Input .Command l Shaping^ 
located in the MIT Artificial Intelligence La ^°J^7i^ c h f a V b e or C ’?^ ?he 

support will augment this cooperation. 

A particular recent design efTort has been the MAOT inertial wheels, 

§li^ 

“ the AI Lab, mid controllf? development and construction was completed 
by SERC personnel. 

Soecific preliminary thesis research began in February 1990 . ine 
first research locus was understanding Singers Impulse Sequenc utinlc 
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yield multiple mode shapers as a direct solution to a set of constraint 
equations. The new sequences will be shorter than the convolved 
sequences, thus minimizing system response delays. The major problem 
with this approach is that the constraint equations are a set of 
simultaneous, non-linear, transcendental equations, which are very 
difficult to solve. 

Non-linear equation solvers perform well when they are provided 
with good initial guesses, hence our next step was to linearize the 
equations, and through an optimization process, generate approximate 
multiple mode shapers. These approximate shapers were interpreted and 
fed into a standard non-linear equation solver as initial guesses for the 
exact solution of the full equations. The exact sequence, theoretically, could 
be used to eliminate multiple mode residual vibration in a flexible system. 

To • validate this concept, we employed a finite element MatLab model 
of MACE. After identifying the system eigenvalues and constructing an 
input shaping sequence, we conducted simulations of response to shaped 
and unshaped slewing maneuvers. The input shaper removed about 90% of 
the vibration in the system, identifying serious potential for the shaping 
method. b 

xtaoa'X 6 ^ similar tests were collected and presented at the 1990 
NASA/DoD Controls-Structures Interaction Conference, and also 
submitted to the 1991 IEEE Conference on Robotics and Automation and the 
1991 Automatic Control Conference (see the attached paper for details). The 
response at the CSI conference was favorable, and several attendees stated 
intentions of using the technique on some of their test articles - a success 
considering our major aim of raising the technique's level of familiarity 
within the aerospace community. 

Current work is aimed at increasing the number of vibrational 
modes that can be eliminated by a single input shaping sequence. We have 
h&d success at raising this number to five modes for some of our current 
models; we hope to develop shapers capable of eliminating vibration in up to 
ten or fifteen modes. The major obstacle at this point is the robustness of 
the current non-linear equation solver, and we are going to try various 
other solvers in the coming weeks. 

Future activity will include the use of non-linear models of MACE. 
The equations defining the shapers are derived from linear flexible system 
and the exact effects of system non-linearity on the input shaper’s 
effectiveness are unclear. The new models will allow us to identify some of 
the problems caused by non-linearities, and hopefully assist in developing 
approaches for dealing with those problems. Sometime in March or April 
we hope to be conducting tests on the physical MACE experiment, still 
under construction at SERC. 
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NONLINEAR ANALYSIS OF COMMAND SHAPING 
Mr. Ken Chang and Prof. Warren Seering 
(Multibody Testbed) 

Tt was determined that Jim Hyde’s work on input shaping required a 
"TtSfne MSCOS o^SUN S^a^ wor£^ la “jhe^d^nd 

^lab'sf 1 upgraded ^UNIX^ystem * Th^ progranf is™^ working ’end’ resides 
on the hard disk. 

The next logical step in this non-linear modeling work is to include 
iU • ? eUonor in thp model This requires first finding the natural 

$£ VfpSyeaSfy tojfemeAtfl by Lnipufating/shaping the input 
before it goes into the DISCOS program. 


High bandwidth control of low Area density 

DEFORMABLE MIRRORS 
Mr. Eric Anderson, Mr. Jonathan How, 

Prof. Steven Hall and Dr. David Miller 
(Possible Third SERC Testbed) 


Introduction 

A Droiect to investigate the applicability of Controlled Structures 
The’objective "worffs step" controf' "technolo^ 

class of 

applications. 


Problem Statement 

A hrief survey of the available literature shows that designs to control 
mirror^ surfaces haveexisted for over twenty years*. A little more analysis 
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by controlling the surface in its quasi-static regime. The flexible modes are 
tound at a high frequencies, beyond the required bandwidth. Consequently, 
tor the purposes of the control design, it is possible to treat the mirror as a 
collection of rigid panels and ignore the flexible dynamics. Many of these 
designs have been developed using a combination of wavefront and edge 
sensors, piezoelectric actuators, and a static figure control feedback. This 
approach assumes that the control roll-off would be such that the flexible 
modes of the stiff structure would not be destabilized. Very few of the 
available control approaches discuss the topic of accounting for the effects of 
control spillover. However, it is recognized that the mass per unit surface 
area of the mirrors must be reduced from about 200 kg/m 2 to 50 kg/m 2 
before the much larger designs envisioned for future space missions 
become feasible. Weight reduction will result in a much more flexible 
structure. Future control designs must account for this increased 
flexibility to avoid spillover destabilization of these modes. 

For this work, we will specifically be looking at reducing the weight 
of the structure even further with an emphasis on reducing the 
manufacturing (e.g. face grinding) costs by actively compensating for long 
wavelength flaws in the mirror surface. The objective is to use a much 
cheaper high precision active structure for the mirror combined with a 
large amount of active control to achieve performance comparable to the 
previous designs. The resulting mirror will then be better in terms of both 
cost and weight. 

Control systems for current mirrors typically use actuators which 
push off a backplane which is of comparable stiffness to the mirror. This 
backplane typically constitutes a significant contribution to the total mass 
of the mirror structure. Our aim is to develop a self-straining structure 
(which eliminates the need for the backplane), and the corresponding 
control architecture to improve the shape performance of a very flexible 
primary mirror over a wide range of frequencies in the presence of 
disturbances and modelling errors. The controller will be based on 
feedback of a combination of structural and optical measurements. A 
typical representation of the relative control authority of these two feedback 
loops is shown in Figure 1. Note that there are (at least) two main sets of 
sensors, but probably only one set of actuators. As shown, the feedback 
from the structural measurements will be used to control the mirror over a 
wide frequency range. In particular, it is important to be able to minimize 
the bandwidth f\. This will allow for a significant overlap between the two 
control approaches even when the target object is a faint star and photon 
limitations result in a very low wavefront sensor update rate. At the higher 
frequency end, the aim will be to improve the shape performance of* the 
modes which will contribute to the RMS ripple error but typically have a low 
damping ratio. 


Overall Approach 

•ii l This research will develop analytical solutions to thE problem which 
will be verified on an experimental testbed to be designed and constructed 
Ihe experimental issues will deal with the final design and lay-up of the 
test structure, the identification of its properties, and tne implementation of 
the controllers. This work can also be separated into the design of the 
controller and the structure, but it is recognized that the boundary between 
the two is not distinct since both will have to deal with issues such as the 


sensor end actuator type and architecture The approaches that will be 
taken are discussed in tne following two sections. 


APPROACH TO THE STRUCTURAL DESIGN 




UAVU — _ _ 

rm, a a^miptnral design will proceed in close coordination with the 

complicated flexible controlled structures, 
struct^ 

active structures are best suited established based 

conventional approach of „ rnC essine requirements, reliability, 

and^imDUcit^of^^i^and^anu^ctiTre'Hpe^nnance objectives for 
precision^ control will be incorporated in the structural design 

and distribution of active components in th “ t ™^ r . e _ diaturbanci 
on the proposed control atcbiUcture “ d Kacy of sensors 

Actuators ZSlXTvahXed. with allowance for some redundancy and 
fault tolerance. , A 

The structure will be Sd 

SSSMXg W^e JoredTn conjunction with ongoing efforts in 
SERC 

In parallel with the analytical tural 0 se naor and actuator 

presence of disturbances to the best precision possible. 


APPROACH TO THE CONTROL DESIGN 


aw ***- , 4 , 

Thp control design for this testbed will consist of determining the 
solution'll ~m 6 plea problems 

aC ^\ t0r VnTiXrm ta an 0n a'ndniO evaluati^ of the best alternatives by 
will be to perform an an J implementation costs versus the 

performance 3 Ond^nTinT in^perimOntal comparison of the Anal 

designs. 1 

There are several factors about this problem that will strongly 

Srn SrVr^ ar^aTracTion of the 
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I. 

II. 

III. 


wavelength of light. This typically will require a very dense array of both 
sensors and actuators. Consequently, the number of measurements 
available to the controller will be very large for any reasonably sized 
problem, so the control architecture must be designed to efficiently handle 
this information. Since the aim is also to design a controller for a laree 
frequency range (see Figure 1), either a good model must be developed for 
the higher frequency regions or the controller must be designed for 
robustness to parameter uncertainty. Another factor in the tradeoff 
analysis is that the control commands will be calculated based on the 
measurements from at least two sets of different sensors, so the controller 
must be able to combine these in an intelligent fashion. 

approach that will be taken in this work will be to investigate 
various control architectures, and in the process develop a list of the 
advantages and disadvantages of each in terms of the issues given above. 
Uther cntena, such as the communication requirements, the flexibility of 
redesign after failure, and the performance with noisy measurements will 
also be addressed. Some of the control architectures that will be 
investigated are: 

Centralized design 

Fully decentralized design 

Decentralized design with communication allowed between 
controllers 

IV. Multi-level design 

where the decentralization will probably be at the mirror segment level. 
Clearly, each of these designs has some obvious advantages and 
disadvantages. A centralized design would provide very good overall 
performance, but will be hard to implement due to the amount of 
information available. Fully decentralized designs are much simpler to 
lmpfoment, but have performance problems. Multi-level designs combine 
many of the advantages of these approaches, but there are potential stability 

§1!?? * m ?- C * 0r \ t i r0 • 8 P lllov ? r - It is expected that some form of 

decentralized controller m a multi-level architecture will provide the best 
overall performance for this type of problem with a large amount of sensor 

§*T P rev \ ou ! design studies have certainly indicated 
“V 8 the case 5 - 6 ’ 7 . An analysis based on these points should eliminate 
Potential designs, and further computer simulations will indicate 
which designs are worth implementing on the final testbed. 

.. Son^e of the other issues to be addressed during the initial stages of 
the control design include a determination of the anticipated noise levels in 
the measurements, since this might require that the controllers include a 
. kI'jV 11 a ^ e braic expression for the performance objective will also have 
to be developed. This could be m the form of compensating the modes of the 

°u by c ° ntrol bng the optical modes of the mirror through 
* u . 8e °C t b e Zermke polynomials. 8 Once the final designs are selected, a 
implement the controllers (i.e., digital, analog, or digital 
analog) will 1 have to be made. The controllers will then have to be 
redesigned for this specialized implementation case and to include the 
sensor and actuator dynamics. 
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Near-term Plans 


Important tasks to be undertaken early in the project: 

I. 


VUAA V — — - 

Modelling issues such as the feasibility of accurately modeling 
the high frequency modes or the need for the controller to be 
explicitly designed for robustness will be addressed (mode vs. 
wave modelling, active vs. passive damping). 

II The type, location for good performance, dynamic range, 
resolution and dynamics of the sensors and actuators will be 
evaluated, and the appropriate sensor-actuator matching tor 
efficient control will be studied. 

III Efficient combination of low and high bandwidth control, using 
wavefront and structural control sensors, will be investigated 
(where does the break in frequency (Figure 1) occur/). 

IV. An algebraic representation of the performance criterion will 
be expressed. 

V. Several control architectures employing dynamic or static 
feedback approaches will be compared. 

VI. The ID testbed will be constructed using a self- straining beam 
with distributed sensing and actuation. 

VII Mounting of the beam, as well as optimal mass, stiffness, and 
damping distributions will be addressed. 

VIII An optical reference system will be designed to acquire 
displacement measurements to judge the control performance 
(discrete points or continuous surface). 


references 

1) Creedon, J. and Lindgren, A, "Control of the Optical Surface of a Thin, 
Deformable Primary Mirror with Application to an Orbiting Astronomical 
Observatory," Automatics, Vol. 6, 1970, pp. 643-660. 

2) Fanson, J., Anderson, E„ and Rapp, D., "Active Structures for Control of Large 
Precision Optical Systems, "Optical Engineering, Nov. 1990. 

3) Wada, B. (Editor), "Adaptive Structures," Winter Annual Meeting of the 
American Society of Mechanical Engineers, San Francisco, California, Volume 
AD-15, December 1989. 

4) N W Hagood, E. F. Crawley, J. de Luis, and E. H. Anderson, Development of 
Integrated Components for Control of Intelligent Structures, Smart Materials, 
Structures, and Mathematical Issues, pages 80-104, Technomic, 1989. 

5) Aubrun, J., "Theory of the Control of Structures by Controllers, 

AIAA Journal of Guidance and Control, Vol. 3, No. 5., 78-1689R, September- 
October 1980, pp. 444-45 1. 

6) Hall, S., Crawley, E., How, J., and Ward, B„ "A Hierarchic Control Architecture 
for Intelligent Structures," To be published in the AIAA Journal Guidance, 
Control, and Dynamics, (Space Systems Laboratory Report (419-88), November 

1988. 


30 


7) How, J., Local Control Design Methodologies for a Hierarchic Control 
Architecture, SM Thesis, MIT Department of Aeronautics and Astronautics, (Space 
Systems Laboratory Report / # 14-89), November 1989. 

8) Serra, K., Adaptive Optics Control using Zemike Modes, SM Thesis, MIT 
Department of Electrical Engineering and Computer Science, February 1977. 

Objectives and Issues 

Objective: Control to X/20 at 0.6 pm with the lowest area density. 


Wave front 



Static and stiff 


Static and dynamic 
modes important 
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ROBUST PERFORMANCE FOR SYSTEMS WITH PARAMETERIC 

UNCERTAINTY 

Mr. Nesbitt Hagood and Prof. Edward Crawley 
(General Affiliation) 

There has been much renewed interest in the problem of designing 
robust control systems for plants which have highly structured, real 
parameter" uncertainty. This type of uncertainty can be found in flexible 
stmcTures wTre higher modes of the plant can have poorly known natural 
freauendes and dfmping. Since the uncertain modes are usually in the 
critical rolloff region of the structural controller, they can !ead Jto 
degradation of system stability and performance. Such is the case for the 
interferometer structure being developed as a testbed for structural 
controllers In complicated structures like the interferometer the fimte 
element models rarely represent the structural dynamics t<> accuracy 
needed in conventional Linear Quadratic Gaussian (LQG) design 
methodologies. The existing design methodologies must be modified and 
reevaluated in the light of the difficult structural control problem especially 
in the area of robustness to uncertainty. 

In the area of robustness to parametric uncertainty, the f 
technique chosen for evaluation is based on considenng the average 
Derformance of the plant over a set described by the uncertain parameters. 
It has been shown that if the average performance is finite, then no plant in 
the set can have an unstable response. This is intuitrve since if any plan 
was unstable it would drive up the average response. This idea can be used 
to design controllers for uncertain sets of systems by finding controllers 
which Minimize the average quadratic cost. If such a controller can be 
found it will both guarantee stability and increase th ® average system 
nerformance. The idea of using the average has been used in the past by 
Bryson Skelton and others to design robust controllers for aerospace 
systems. The present work has concentrated on providing a formalism to 
average cost design. 

The difficulty with designing controllers to minimize the average cost 
lies in computing the average. For systems with many ^ u nc erta intie s , 
computing the average amounts to performing a high dimensiona 
integration. There are some established ways for approaching thus problem 
such^as Monte-Carlo or Fast Probability Integration but they are to 
computationally intensive for use in control synthesis. The impetus of the 
present work has been to develop efficient ways to compute approximations 
to the average which can then be incorporated into the design. One such 
efficient ^ppfoximation to the average cost, known as the Bouret equation 
was founcfin the area of random wave propagation, turbulence 
and Quantum mechanics. It is ideally suited for calculation of the 
approximate average cost because of its linearity and capability to handle 
large numbers of uncertainty. 

The past six months have been spend developing control 
methodologies base on the Bouret approximations to the average cost f* d a f 
related equation which overbounds the average cost. While minimization of 
the Bouret approximation will provide robustness but not guarantee system 
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stability a priori minimization of the overbound does guarantee a priori 
stability. While it seems that such guaranteed stability would be desirable 
for flexible structures, the guarantee caries with it much lower 

E erformance and high required control effort. The past six months have 
een spent Quantifying the tradeoff between those controllers which 
guarantee stability and those which provide robustness but not necessarily 
stability throughout the set of possible systems. In the coming six months, 
this trade-off will become more clear as the average based controllers are 
first simulated and finaly implemented on the SERC interferometer testbed. 
The two competing controllers, the Bouret approximate average and the 
average bounding controller will be implemented and compared to existing 
robust control methodologies such as and LQG/loop transfer recovery. 
The controllers will be compared on the basis of robustness offered, control 
cost, and performance. 


Area- Averaging Sensors 

Mr. Simon Collins, Dr. David W. Miller and Dr. Marthinus van Schoor 

(General Affiliation) 

Area-averaging sensors which exhibit magnitude rolloff without 
associated phase lag have been experimentally demonstrated. This work 
was published in a paper entitled “Development of Spatially Convolving 
Sensors for Structural Control Applications” which was presented at the 
1990 AIAA Structures, Dynamics and Materials Conference in Long 
Beach, CA. Achieving magnitude rolloff without phase lag enables gain 
stabilization without sacrificing phase margin. Future work in this area 
will first involve demonstrating the advantages of zero phase lag sensor 
over a point sensor, with or without rolloff dynamics, in an actual closed- 
loop experiment. 

Work during the last six months has focused on the spatial nature of 
optimal control solutions. In this effort, full state Linear Quadratic 
Regulator solutions to discrete approximations of distributed parameter 
systems is formulated. It is realized that the finite set of feedback gains, 
associated with measurements of structural motion at various locations! 
are simplv discrete approximations of the continuously distributed feedback 
kernel which convolves with the state function to generate the control 
command. By observing the discrete gains, this kernel can be inferred. It 
is then shown that displacement and rotation feedback kernels can be 
transformed into equivalent curvature feedback kernels. This facilitates 
implemention using the area-averaging sensor concept. The paper entitled 
' Formulation of Full State Feedback for Infinite Order Structural Systems" 
will be presented at the First Joint U.S./Japan Conference on Adaptive 
Structures in Maui, Hawaii. Future work will involve the experimental 
implementation of the optimal feedback kernel. 
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the Need for passive damping in feedback controlled 

STRUCTURE 

Prof. Andreas .H. von Flotow 
(General Affiliation) 

It is well established that an infinite dimension^ mafJem^icaZ 

ms»s SMUssKfi 

control, t Q structural gain at each resonance. Thus, even a 

characteristics of such a problem structure. 

This naper reports upon the enabling effect of passive damping in the 

terms of classical ideas of gain and phase stabilization .The .paper derives 
compensation. 

1990 ) Tn^hictf ^Mceptabfe°interaction e of^a* «S« Vt^SugTE? .j£ 

required t0 
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Other SERC Activities 


Middeck o- Gravity dynamics experiment(MODE) 

Prof. Edward Crawley, Dr. Marthinus C. van Schoor, Dr. David W. Miller 
MODF is an STS middeck flight experiment funded under the NASA 

SEs =;= s=“~ 

structure. . . 

The MODE program completed its Critical Design Review ( A 
March 1990^ Construction of flight hardware is proceeding with completion 

in January, 1991. 


Middeck Active control experiment (MACE) 

Prof. Edward Crawley and Dr. David W. Miller 
MACF was a Phase B study, funded by NASA Langley Research 

SERC multibody testbed. 
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Government and Industry Interaction 


ITEK INVOLVEMENT 

As a result of the Center, SERC has become aware of the research 
and development in active shape control being performed at ITEK and ITEK 
has become familiar with the research effort in the Center. This has led to 
the evolution of a joint research program in the development, of lightweight, 
active mirror technology. An actively deformable, precision surface is the 
third SERC testbed ana is presently in the conceptual phase. 


JPL Involvement 

SERC and JPL have continued their close relationship by working on 
research into active struts, fostered by SERC graduate students residing at 
JPL during the past year, and through an effort to understand the 
relationship between certain similar research results achieved at SERC 
and JPL. 

w . „ JPL, through Dr. Robert Laskin, has started to participate in the 
Middeck Active Control Experiment (MACE). 


NASA LANGLEY RESEARCH CENTER 

While coordination with efforts at NASA LaRC has been good in the 
past, it was recently suggested that such coordination could be 
strengthened. The present plan is to better coordinate SERC's multibody 
LaRC s evolutionary CSI testbed. The August 1990 Summer 
SERC Symposium was held at NASA LaRC. 


McDonnell Douglas Space Systems Company 

MDSSC is continuing its participation in the Middeck 0-Gravity 
Dynamics Experiment (MODE) and is part of the Middeck Active Control 
Experiment (MACE) team. 


Lockheed missiles and Space company 

LMSC has joined the Middeck Active Control Experiment (MACE) 
$ e »^r^ an “ 8U PP° rte d several undergraduate and graduate students at 
LMSC during the summer of 1990 to work on MACE component 
technologies. LMSC plans to significantly increase involvement in MACE 
in calendar year 1991. 
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Milestones 


PAST 


Event 

Objective 

Date 

SERC “ 
Symposium 

Presented 
research (held at 
NASA LaRC). 

August 1990 

Multibody 
Testbed Critical 
Design Review 

Formalize design 
of multibody 
testbed. 

February 1990 


FUTURE 


Event 

Objective 

Date 

Steering 
Committee and 
Technical 
Representative 
Meeting 

Review the SERC 
program 

January 1991 

Science Advisory 

Committee 

Meeting 

Review the 

interferometer 

testbed. 

Review science 
issues for 
multibody 
testbed. 

Spring 1991 

SERC 

Symposium 

Present research 
(to be held at 
M.I.T.). 

June 1991 
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Abstract 

Spacecraft, space-borne robotic systems, and manufacturing 
equipment often utilize lightweight materials and configurations that 
give rise to vibration problems. Prior research has led to the _ 
development of input command pre-shapers that can sl g^jfic y 
reduce residual vibration. These shapers exhibit marked insensitivity 
to errors in natural frequency estimates and can becombined to 
minimize vibration at more than one frequency. This paper a 

method for the development of multiple mode input shaoers which 
are simpler to implement than previous designs and produce smaller 
system response delays. The new technique involves the solution of a 

poup of simultaneous non-linear impulse "an^T/^SA^ 

resulting shapers were tested on a model of MACE, an MIT/NAbA 

experimental flexible structure. 


Introduction 

Space-borne robotic systems and vehicles often employ lightweight materials 
and configurations that result in a high degree of system flexibility. The 
system's light weight facilitates launching, but chronic vibration problems are 
a common result. Manufacturing equipment also increasingly utilizes lighter 
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structural elements, a main objective being improving the speed of 
automated assembly. The combination of a lightweight structure with high 
performance requirements often leads to serious vibration problems. The 
growing demand for high accuracy manipulation is in no way aided by these 
simultaneous attempts to increase speed and decrease weight. Partial or 
complete suppression of system vibration can improve spacecraft durability 

and performance, and would allow manufacturing systems to operate faster 
and more economically. 

Attempts to decrease the vibration inherent in flexible systems have 
enjoyed varied success over the past decade. Cannon and Schmitz [1] 
experimented with the non-colocated feedback control of a flexible beam. 
Through the use of accurate system models and optical tip position sensing 
they achieved significant vibration reduction in their planar test article. 

Yurkovich and Tzes [81 reduced vibration in the presence of unknown 
and/or varying payloads by employing on line system identification and 
controller tuning. By using frequency domain techniques to examine the 
system response following a sample input, enough information was gained to 
adjust the controller gain scheduling to compensate for vibration problems. 

Wie [9] employed H« controllers to reduce vibration while providing 
robustness to modelling errors. This technique displayed solid performance, 
but was relatively difficult to implement 

Input command shaping is an attractive vibration reduction method 
because it is essentially "hands off;" inputs can be fed through a shaper and 
into the system, and ideally the resulting system output will be vibration free. 
Shapers also usually reside completely outside of a given control system and 
are thus easily compatible with other vibration schemes (see figure 1). Smith 
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Input 


Gosed Loop System 



Output 


Figure It Shaper Position in Control System 


16] conducted early shaping investigations, largely through the use of posicast 
control. 

Meckl [3] examined the use of shaped force profiles to reduce vibration 
in manufacturing systems. Meckl created profiles by using a versine ( 1 - 
cosine ) function to modify force commands. When integrated twice, these 
force profiles became input trajectories that reduced system vibration at a 

structure’s first natural frequency. 

A major problem with command shaping is that its success usually 

depends on solid prior knowledge of plant dynamics. Many attempts at input 

shaping have been criticized because the shapers exhibited significant 

dependence on precise system models. 

Singer [4] presented a simple shaping algorithm that demonstrated 
strong insensitivity to modelling errors. The shapers were assembled from 
impulse sequences and produced only small delays in system response times, 
on the order of one period of a system’s natural frequency. This technique 
performed notable vibration reduction in tests of a full scale mockup of the 
Space Shuttle Robotic Manipulator System, conducted at NASA s 
Manipulator Development Facility at the Johnson Space Right Center. 
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Tzes, Englehart, and Yurkovich [7] studied the effects of combining 
input shaping with a closed loop, acceleration feedback controller. They 
conducted tests on a flexible beam and proved that each technique can 
complement the other, resulting in enhanced vibration reduction. This work 
supports the assertion that input command shaping can be used concurrently 
with other vibration suppression schemes. 

Singer [4] originally assembled shapers designed to cancel single mode 
vibration and later expanded the algorithm to handle multiple mode 
problems. The initial multiple mode technique was somewhat cumbersome, 
however, and the main purpose of this paper is to present an improved 
method for developing multiple mode shapers. Simpler impulse trains can 
be assembled by directly solving a full set of multiple mode vibration 
equations. These new shapers have all the vibration reduction capabilities of 
the original shapers, and yet exhibit savings in implementation complexity 
and response time. We present an approach for solving the vibration 
equations and offer evidence of the new shapers' potential through tests 
conducted on a model of MACE, an MIT/NASA experimental flexible 
structure. 

Single Mode Shaping 

To develop a single mode input shaper, we first note that the second 
order system response to an impulse input is described by: 

yitt = Ai sin((/-ti)coV l-^ 2 ) (1) 

where yi(t) is the output, Ai is the impulse amplitude and ti is the time at 
which the impulse occurs. The system's vibration frequency is co, with 
damping If the system is linear, its total response to a series of N impulses 
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can be expressed as a sum of the the responses to each impulse "i." The 
gnitude of the total response following the Nth impulse is given by: 


ma« 


Amp = Ai sin(tid)V 1 - ? )) 2 + 

(S, A i cos(tiO)V 1 - C 2 )) 2 J 


1/2 


( 2 ) 


A train of properly arranged impulses can suppress residual vibration 
by forcing Amp to equal zero. This can only happen when both the sine and 
cosine terms in equation (2) independently equal zero: 

£ Aj e^ 0 * sin(tiCoV 1 - C ) = 0 

i=i 

N r 5" 

jT A, cos(tiCoV 1 - C ) = 0 

i=i 

To construct an impulse sequence that will act as a vibration reducing 
input shaper, we start by imposing two initial constraints: 


ti-0 

N 

X Ai-1 

i = 1 


(4) 

(5) 


The first is simply an origin specification, and the ¥ cond is a normalization 
constraint. Normalizing a shaper's impulse magnitudes ensures that a 
shaped input will not exceed limitations imposed on the original input, such 
as actuator or stress limits. We specify an arbitrary value for Al, and with N 
= 2, we can use equations (3) to solve for the time and amplitude of the 

second impulse in a two-impulse shaper. 

This shaper will completely cancel residual vibration in a single mode 

system, as long as the natural frequency and damping ratio are perfectly 
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known. To account for possible modelling inaccuracies, the shaper should 
exhibit some insensitivity to errors in natural frequency and damping ratio 
estimates. By differentiating equations (3) with respect to natural frequency, 
we generate two additional impulse constraints: 


1 / — 

X Aj tj e5“* sin(tjG)V 1 - £ 2 ) = 0 


N 

X A 


i=l 


i t, eS°*> cos(t,o)V 1 - £ 2 ) = o 


( 6 ) 


Setting the partial derivative with respect to natural frequency equal to zero 
also sets the partial derivative with respect to damping ratio equal to zero [4] 
These new constraints require the addition of a third impulse to our 
sequence; we have four equations, we need two unknown amplitudes and 
two unknown times. The three impulse sequence will force the residual 
vibration to be low even if the system parameters are not precisely known. 

The standard three impulse, single mode shaper features impulses 
with a 1-2-1 magnitude configuration and times that are equally spaced. A 
typical sequence is shown in figure 2. 


Amplitude 

(Normalized) 


Time (sec) 

Figure 2: Typical single mode three impulse shaper. 
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Note: these impulses are constrained to having positive amplitudes. 
By using negative impulses, the time of a series' final impulse can be 
decreased, but negative impulses tend to tax a system's actuators and 
introduce high stress levels. In the remainder of this paper, all shapers will 
utilize impulses with positive amplitudes. For a full derivation of the above 
equations, see Singer [4]. 


Adding Modes 

To cancel multiple mode vibration, we can convolve several single mode 
impulse sequences into longer trains. Convolution results in a sequence 
whose final impulse is located at a time equal to the sum of the damped 
periods of the cancelled modes. The value of the final impulse s time will be 
referred to as the shaper's ’length.'' The number of impulses in the 
convolved sequence is equal to 3 m , where m equals the number of modes. A 
standard three mode convolved shaper is shown in figure 3. This sequence 
was solved for a zero damping case, so the twenty-seven impulses are 
arranged symmetrically about the center of the pattern. 


Amplitude 

(Normalized) 



Figure 3: Three mode convolved impulse sequence. 
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The convolved multiple mode sequences are easily generated, but their 
failings become clear when the cancellation of higher mode vibration is 
desired. The number of impulses in a convolved shaper increases 
exponentially with added modes, and by the time the third or fourth mode is 
added, the sequence has become packed with impulses and can be difficult to 
implement in real time. Shapers with more impulses increase the time 
required to modify an input, and might force a decrease in servo rate. 

The solution to these problems is to build multiple mode sequences 
not through convolution of single mode sequences, but through a direct 
solution of the constraint equations, (3) and (6), as written to include an 
arbitrary number of modes: 

N , 

I A i e^sin(t l to j Vl-;J )-0 

I»1 

N y 

Z A ( eSw cos(tj£DjV 1 - Cj ) - 0 

iml (7a) 

N 

L Ai ti c Sw sin(t,Ci)jV 1 - £ 2 ) = o 

i-l 

Z Ai ti eCw* cos(t,o>jV 1 - ) a 0 

ial (7b) 

Repeating equations (7) for additional modes "j" generates a set of 
simultaneous non-linear impulse expressions. Solving these equations can 
yield shapers with shorter lengths than the convolved shapers. Shorter 
sequence lengths decrease the delay in system response caused by using the 
shaper. The direct solution sequences, moreover, use only (2*m) + 1 
impulses, m being the number of cancelled modes. This linearly increasing 
impulse population leads to vastly fewer impulses in higher mode shapers, 
reducing implementation time. 
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The savings in length and impulse density that the direct solution 
sequences support are offset by an increase in sequence generation complexity. 
The single mode shaper equations, given the constraints of positive impulse 
amplitudes and shortest possible overall length, had a closed form solution. 
As written for the multiple mode case, the shaper equations require a strict set 
of constraints just to limit their infinite solution space, and no general 
solution has been found. 

Solving the Equations 

The key to solving the multiple mode equations thus far has been to employ 
a linear approximation. Equations (7) are non-linear only in terms of 
impulse time. A straightforward approach is to pick a time for the sequence s 
final impulse, essentially defining a sequence length, and then divide the 
length into a fine time mesh. An impulse is placed at each time slot, with 
unknown amplitude but known time. The equations are now under 
constrained, but a linear approximation to the exact shaper sequence can be 
generated through optimization. 

The constraints for the optimization problem are the multiple mode 
equations (7) and the normalization requirement (5). The cost function is the 
sum of the second derivatives of equations (7a): 

Cost * | Z Z A i t? sin(titi)j V 1 - Cj ) + 

j-l i-l 

Z Z Ai t? cos(tiCi)jV 1 - Cj ) I 

J-l i-l (8) 

Minimizing the second derivative expressions forces the impulse sequence to 
be even more insensitive to modelling errors. With these guidelines, the 
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linear problem becomes "minimize the cost function subject to the stated 
constraints." 

Solution of the linear problem was accomplished using GAMS, a 
standard optimization package. GAMS utilizes a version of the primal 
simplex method to perform linear optimization. If M is the number of 
modes to be cancelled, o is the length of the sequence, and dt is the value of a 
single time mesh element, GAMS' constraint matrix consists of: 

rows: r = (4 * M) + 1 
columns: c = a / dt 

The variable vector is the series of impulse amplitudes, A[. The simplex 
method dictates that at least (c - r) amplitudes will equal zero, and additional 
impulses are occasionally set with zero amplitude. The optimized GAMS 
output yielded an impulse train with a number of impulses that was less than 
or equal to r. This train was a linear approximation to the exact multiple 
mode shaper. 

The second phase of the linear work was to find the feasible solution 
with the smallest possible final impulse time. This was achieved through 
multiple GAMS runs, systematically reducing the time of the final impulse 
and using a binary search algorithm that recognized when GAMS returned an 
infeasible solution, meaning that the time had been reduced too far. This 
technique ensured that the final GAMS output was the shortest possible 
approximation. Figure 4 shows a typical final GAMS result from a three 
mode problem. 
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Amplitude 

(Normalized) 



Figure 4: GAMS output for a standard three mode system. 


In many cases, this output could stand alone as an effective shaper 
sequence, especially if the time mesh was set at a digital controller's servo 
rate. The raw GAMS output, however, had about twice as many impulses as 
the exact solution demanded. Our goal was to arrive at a sequence with as 
few imp uls es as possible, so we used the GAMS output as an initial guess in a 
non-linear equation solver. 

Given the large number of impulses in the final approximate sequence, 
the GAMS output had to be interpreted to obtain useful guesses of the exact 
solutions to equations (7). GAMS would often place impulses in adjacent 
time intervals; these impulses were replaced by single spikes that combined 
the amplitudes of the neighbors and adopted their exact average time. To 
further reduce the number of impulses, the interpretation algorithm sought 
out the closest non-adjacent neighbors. These pairs were combined by 
summing their amplitudes and taking a weighted average of their times. 

This set of techniques yielded a sequence whose number of impulses matched 
that required by the non-linear multiple mode shaper equations (7). 
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Mathematica™ was used as the non-linear equation solver, mainly 
because of its additional potential as a programming language that could 
envelop the entire computational side of the impulse sequence generation. 
The time and amplitude of the first (i = 1) impulse in equations (7) were held 
constant, matching the first impulse from the interpreted GAMS output. The 
remaining times and amplitudes were allowed to vary, with initial guesses of 
their values provided by the reduced GAMS sequence. 

Mathematica employs a Newtonian gradient search algorithm to arrive 
at its solutions. This method worked quite well, as long as our guesses were 
sufficiently close to the optimal solutions. As the non-linear equations are 
continuous and differentiable, adequate gradients are readily available, and 
points of singularity are usually easy to avoid. The resultant exact impulse 
sequence, after interpreting the GAMS output in Figure 4 to find initial 
guesses, is shown in Figure 5. The impulses in the Mathematica result were 
re-normalized to ensure that the constraint of equation (5) was upheld. 


Amplitude 

(Normalized) 



Figure 5: Exact impulse sequence solution, output from Mathematica. 
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We found that Mathematica could solve the equations only part of the 
time. A particular structure's absolute modal values and relative modal 
spacing could disrupt the GAMS program, or Mathematica, or both, resulting 
in a group of unsolved equations. A variety of different approaches have 
been used to increase the robustness of the solution process, the ultimate goal 
being the discovery of a closed form solution. While current work is 
continuing in this area, the linear approximation/non-linear solution 
algorithm has successfully generated workable impulse sequences for 
different groups of three and five modes. 

Modeling and Results 

The three mode case of particular interest involves a set of frequencies found 
using a model of an actual flexible system, the MACE test article. The MACE 
experiment is a joint MIT/NASA project designed to study methods for 
controlling flexible systems in micro or zero gravity fields. MACE is a flexible 
structure with two multi-axis pointing payloads residing on either end of a 
tubular bus. The system incorporates attitude control through a set of three- 
axis torque wheels, and utilizes inertial position sensing information gained 
from gyroscope packages mounted at the center of the bus and inside each 
payload. A simple system schematic is shown in Figure 6. 
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The MACE project consists of a Ground Test Article, to reside at MIT, 
and a Flight Test Article, scheduled for a Space Shuttle launch in 1993. The 
Ground Test Article is currently being assembled and should be available for 
testing in December of 1990. The ground article will be actively suspended to 
emulate a flexible spacecraft in a micro gravity field. This experiment is a 
prime candidate for practical validation of the command shaper techniques. 

While the physical MACE structure has been under construction, 
personnel at MIT's Space Engineering Research Center (SERC) have 
developed several computer models of the experiment. The frequencies used 
in the three mode case mentioned above were found using a linear finite 
element MatLab model of MACE. The planar model depicted the segmented 
bus and one of the pointing payloads, as shown in Figure 7. 


Rotating Cantilevered Mass 
(Payload) 



Segmented Beam 
(Bus) 






Payload Rotation 


Nodal Displacements and Rotations 





b. W 


w- 


F n M„ 


Figure 7: Finite element model of MACE. 
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The first three eigenvalues of the model were fed into the GAMS / 
Mathematica routine, generating the shaper sequence shown in Figure 5. 

Next, simple torque inputs were fed through the shapers and into the 
modelled payload's gimbal axis, producing the adjusted inputs shown in 
figure 8. The resulting translation of the beam element on the opposite end 
of the bus is shown in figure 9, and detailed views of the unshaped and 
shaped response are provided by figures 10 and 11, respectively. The model 
had a system of eight modes of vibration, and only the first three were used in 
forming the input shaper. It is clear from the figures, however, that 
cancelling these three modes was sufficient to suppress the majority of the 
structure’s vibration. 
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Figure 8: System inputs adjusted by the input shaper. 
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Figure 10: Response to unshaped input (detail). 
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Figure 11: Response to shaped input (detail). 


Conclusions 

These MatLab results are somewhat predictable. The input shapers are 
defined by equations that predict the response of linear systems, and the 
MatLab model was also linear. Cancelling the vibrations of the MatLab 
model, therefore, served mainly as a confirmation of the proper solution of 
the constraint equations, and allowed for concrete visualization of what a 
system experiences when the input shapers are employed. The lessons 
learned from this initial case will also be valuable When more complex, 
higher mode shapers are developed. 

The next step in this program is to employ more accurate models of 
MACE. The test article has been simulated non-linearly, using the DISCOS 
program. This model will likely predict some of the shaper’s failings in 
suppressing vibration in non-linear systems. 

The second major future task is to improve the equation solving 
algorithm to facilitate the construction of higher mode shapers. Sequences 
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that can cancel up to ten or fifteen modes are not out of the question. In 
addition to increasing the number of cancelled modes, we are devoting effort 
to decreasing the sequence generation time. We have shown that the direct 
solution sequences are easier to implement than the convolved sequences, 
but they are much more difficult to generate. These continued efforts to solve 
the constraint equations, coupled with the lessons learned from the DISCOS 
model, will aid in the generation of input shapers capable of effectively 
reducing vibrations in the actual MACE structure. 
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THE MIDDECK O-GRAVITY DYNAMICS EXPERIMENT 
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generalized coordinates. 

Derive governing differential equations by applying 
Lagrange's Principle. 

Solve nonlinear equations subject to harmonic excitation. 
Verify predictions with MODE flight and ground test results 
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THE MIDDECK ACTIVE CONTROL EXPERIMENT (MACE) 
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Substantial commonality of ESM hardware/software 
Significant savings in integration/certification process. 
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O 


tn 

0) 

£ 

y 

£ 

£ 

cd 

£ 


4-3 



© 

£ 

cd 

*- 

<D 



CAPTURING THE ESSENTIAL PHYSICS: TEST ARTICLE 
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sufficiently complex geometry such that the test article 
undergoes full 3-D kinematic and coupled flexible motion. 


GROUND BASED ENGINEERING MODEL TESTBED 
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5) More complex test article geometries can be used to further 
challenge ground testing. 

6) Fundamental to the flight experiment is the understanding or 
gravity influences on the problem. 
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Nonlinear Three-Dimensional Model 
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GRAVITY INFLUENCES 








THREE CLASSES OF CONTROL ALGORITHMS WILL BE 
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associated control algorithms. 

- These tests identify the ability to which the test article's 
control algorithm can be tuned on orbit. 



VALIDATION OF THE ANALYTICAL FRAMEWORK 




VA tjdation of the analytical framework 

• Ground testing will include suspension and gravity effects, 
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Independent of absolute pertormance ieve 
CST effectiveness. 

Sensor dynamic range of 60 dB is typical. 
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SUMMARY 

The MODE family of flight experiments is designed to verify 
analytical tools developed to predict the gravity dependent 
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Introduction 

The MIT Space Engineering Research Center (SERC) has developed a controlled structures 
technology (CST) testbed based on one design for a space-based optical interferometer. The role 
of the testbed is to provide a versatile platform for experimental investigation and discovery of CST 
approaches. In particular, it will serve as the focus for experimental verification of CSI 
methodologies and control strategies at SERC. The testbed program has an emphasis on 
experimental CST-incorporating a broad suite of actuators and sensors, active struts, system 
identification, passive damping, active mirror mounts, and precision component characterization. 

The SERC testbed represents a one-tenth scaled version of an optical interferometer concept 
based on an inherently rigid tetrahedral configuration with collecting apertures on one face. The 
testbed consists of six 3.5 meter long truss legs joined at four vertices and is suspended with 
attachment points at three vertices (Figure 1). Each aluminum leg has a 0.2m by 0.2m by 0 25m 
triangular cross-section. The structure has a first flexible mode at 3 1 Hz and has over ^0 global 
modes below 200 Hz. The stiff tetrahedral design differs from similar testbeds (such as the JPL 
Phase B) in that the structural topology is closed. The tetrahedral design minimizes structure 
deflections at the vertices (site of optical components for maximum baseline) resulting in reduce 
stroke requirements for isolation and pointing of optics. Typical total light path length stability 
goals are on the order of X/20, with a wavelength of light, x, of roughly 500 nanometers [1]. It is 
expected that active structural control will be necessary to achieve this goal in the presence ot 

disturbances. . _ .... _ . 

A unique feature of the SERC testbed is the implementation of a muln-axis laser metrology, 
incorporating complex bends in multiple beam path lengths. At three mock siderostat locations are 
precision three-axis active mirror mounts. The fourth vertex holds a laser head and other optics. 
These optical components provide laser interferometric displacement measurements for baseline 
metrology (six axes define the position of the mock collecting apertures relative to the fourth 
reference point). We are concerned that the testbed represents a scaled model of an actual scientific 
observatory as closely as possible. At the same time, we seek to perform CST research which is 

generic and applicable in different areas. , 

The structure is instrumented with accelerometers, load cells, strain gages, experimental 
piezoceramic and piezopolymer sensors, and (initially) three piezoceramic active strut members 
The stiffness of the active struts has been selected to approximately match the impedance of 
structure as seen by the actuator at the active strut mounting location, leading naturally to contro 
designs based on passive shunting, wave impedance, or balanced bridge feedback. 

A finite element model of the testbed was constructed and a conventional system 
identification using an external excitation source will be carried out. The results (frequencies, 
mode shapes) will be compared and the subsequent roles of each of these models in the control 
design determined. Because of inherent inaccuracies of the finite element model in representing 
lightly damped closely spaced modes, the experimentally determined modal model is preferred tor 
control design. Methods for generating uncertainty information from the system identification for 
application in robust control methodologies, and studies of model reduction techniques are 
planned. 

* Research Associate ** Graduate Research Assistant 
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Three-axis active mirror mounts have been designed which provide ±3.5 microns of stroke 
over a frequency range of 500 Hz. Two mounts employ conventional piezoelectric actuators; the 
third mount utilizes electrostrictive actuators that exhibit superior bidirectional repeatability, a 
result of greater linearity and reduced hysteresis as compared to piezoelectrics. The moving mirror 
mass has been sized to reflect the approximate scaled masses of siderostats of the proposed space- 
based optical interferometer. The actual moving mass of the the mirrors will be varied to determine 
the level at which interaction with the structural flexibility becomes significant. 

The remainder of the paper begins with a description of the optics portion of the testbed. 
Then the testbed CST program is reviewed with attention focussed in six areas: results from other 
research closely-related to the testbed, finite element modelling, system identification, passive 
damping, an axial component tester, and control experiments. 


Optics 

In this section, the optical components of the testbed are described. The focus here is on 
the implementation of the on-board metrology system. Functional explanations of space-based 
interferometry can be found elsewhere in this volume. 

Beam-combining coherence requirements for an actual space-based interferometer will 
require on-board sensing and correction mechanisms capable of controlling path lengths to X/20. 
Multi-aperture non-interferometric imaging instruments with similar baselines and operating 
wavelengths can have more demanding requirements. The sensing system for orienting the 
instrument relative to an external reference coordinate frame should have resolution and stability on 
par with the resolution and stability of the internal metrology system. Our immediate concern is the 
reduction in errors due to flexibility (Figure 2). 


I 
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Figure 2: Sources of Path Length Error 

The Interferometer CST Testbed under construction at MIT addresses the problem of most 
direct relevance to CST: control of the instrument geometry in order to control projected baselines 
and internal path lengths. The testbed control goal is to maintain fixed distances between points on 


the structure which represent collecting apertures (mock siderostats) and metrology nodes, since 
relative motion among these points changes both projected baselines and internal path lengths. 

A sample interferometer mission to image a tenth magnitude object at visible wavelengths 
with one milliarcsecond resolution using one meter apertures leads to path length stability 
requirements of approximately 80 nm rms. The basic testbed configuration is intended to include 
enough detail to be representative without being overly complex and costly. Many of the features 
may be applicable to other spacecraft requiring precision control. Sensing of the external (rigid 
body) orientation of the testbed and the science optics are not currently addressed, although 
metrology systems for both of these could be tied directly into the on-board baseline metrology 
system with little difficulty. Additionally, each mock siderostat mount includes provisions for 
mounting a small flat mirror with its reflecting surface coplanar with the vertex of the metrology 
system at that siderostat. Such a mirror might form part of a future science optics chain. 

Laser Metrology 

A six-axis laser metrology system forming an optical tetrahedron (Figure 3) will provide 
the primary measure of control effectiveness. One vertex is located at each of the three mock 
siderostats with the fourth vertex containing the out-of-plane reference point. The outputs of the 
near and far legs will yield relative displacements among the vertices with the minimum number of 
laser axes. The vertices of the optical and structural tetrahedrons typically do not coincide since the 
siderostat locations were chosen to represent non-redundant baselines without necessarily requiring 
rigid body tilting of the entire instrument. In the initial configuration, one siderostat plate will be 
located near one structural vertex; the others are roughly 1/2 and 1/3 of the distance down two 
different legs. The relative angles between the actively-mounted cat's eyes will be less than 96.5 
degrees, which is within the cone of operation. 



Figure 3: Illustration of the Structural and Optical Tetrahedral Trusses 

The power required to operate six axes, instead of the nine that would be needed to 
determine Ax,Ay,Az for the three siderostats, permitted the use of a commercially available laser 
measurement system using a single laser head mounted on the testbed. We are using a dual- 
frequency stabilized laser head (670 pW total power), detectors, and fringe counting electronics 
manufactured by Hewlett-Packard Corporation. A lens and 45-degree polarizer assembly plus a 
short length of optical fiber allow the detector electronics packages to be located out of the way of 



the measurement optics and associated mounting fixtures. The VME-based fringe counting 

electronics provide a seamless link to our real-time control computer. ....... UD 

Figure 4 details one measurement axis. The measurement resolution is limited by the Hr- 
supplied electronics to X/64 at X = 633 nm, or approximately 10 nm. Greater resolution can be 
obtained with alternate electronics, such as the VME modules developed by Mike Shao s group at 
JPL For operation in air without wavelength tracking over short time scales and in a laboratory 
disturbance environment, we feel that 10 nm resolution will be adequate. Our closed loop control 
frequency range is 2-200 Hz, so changes in the refractive index of air and other sources of error 
wfrh long time constants will not pose any problem. A preliminary error budget suggests that 
measurement resolution will be ~ 17 nm. 
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Figure 4: One Axis of the Laser Metrology System 


Cat's eye refforefirc tonf will be used ,o provide wide fields of- view -a. .he veni, ces of Jhe optical 
tetrahedron. These are similar to cats eyes used by C. Townes (UC Berkeley lUp 
interferometer) and D. Hutter (US Naval Observatory Astrometric Interferometer) dthough n this 
application there is no siderostat slew range to contend with The minimumsiz^ 
fora given amount of spherical aberration is a function of the laser beam diameter and the ^fract 
index of the cat's eve glass The metrology laser beam diameter of 6 mm at the laser head led to a 
c*“e TSS? ma™ which was unnecessarily cumbersome for taptaen^on ™ a 
moving platform. Reducing the beam diameter permits the cats eye size to be reduced while 
maintaining the same spherical aberration performance^ Lenses reduce the collimated beam 
diameter to 4 mm without reducing the available power. The cat s eye parameters are. 


glass index at 633 nm 

radius of small hemisphere 

radius of large hemisphere 

max. AOPL across beam cross section 

mass 

usable field of view: 


1.72 (Schott SF 10) 
25 mm 
34.7 mm 


A/10 
511 g 

+/- 60 degree cone (see Figure 5) 
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Figure 6: Layout of Fourth Vertex 





Modifying the radii to generate a small amount of focusing will help counteract the increased 
divergence of the smaller beam diameter and produce better overlap at the detectors. The cu ™f 
surface of the large hemisphere will have a silver reflective coating with a protective overcoat. The 
curved surface of the small hemisphere will be coated with a broadband anti-reflective coati g 
targetted to be the proper thickness for X = 633 nm at half the cone angle. Ana-reflective 
performance at other ingles will depend on the spectral response of the coating. The hemispheres 
will be aligned after coating and joined by optical contacting. 

0thgrC ^h^ re mainder of the optics for each measurement leg consists of a polarizing beamsplitter 
cube with crystal quartz quarter wave plates cemented to opposite faces, plus the associated feed 
ootics EaSbeamsplitter-waveplate assembly is mounted in a semi-custom mount which provides 
Kiu^memSees of freedom needed to align the measurement beam with respect to us 
retroreflector endpoints. Three of these mounts are rigidly attached to an open pyramidal bucket 
that is itself rigidly attached to the main fourth vertex optics plate (Figure 6). The remaining t ree 
mounts Se Sy attached to the siderostat optics plates (Figure 7) in the far leg measurement 
Tths These rigid mounts are designed to prevent motion of the beamsplitter optics from 
appearing as morion of the retroreflectors. The feed optics must maintain the orthogonality of the 
E^zations through complex bends in order to minimize errors due to polanzatton mixing. 
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Figure 7: Layout of Siderostat Plate 





Results from Recent Testbed-Related Research 

In conjunction with the development of the testbed, several other areas of research are 
being pursued. Three studies have been completed, and are documented elsewhere [1-3], Some 
relevant results are summarized below. 

Systems Level Dist urbance Minimization Using Controlled Structures Technology 

Disturbances present on a typical large space-based observatory are detailed. The spectrum 
of disturbances is divided into those which depend on the space (Earth orbital) environment and 
those which are internal. Various CST techniques for minimizing the effect of disturbances on 
mission requirements are reviewed. These include passive structural tailoring, passive damping, 
vibration isolation, and active structural control. The full-scale 35-meter baseline version of the 
interferometer testbed is used as a case study for evaluating the flowdown of systems level 
information to the structural requirements. The power, attitude control, and interferometer and 
metrology subsystems are discussed with respect to their role as disturbance sources. Finally, an 
approach for systems level disturbance minimization is outlined. 

Experim ental Characterization of Damping at Nanostrain Levels 

In light of the increasing trend towards nanometer-level requirements on structural stability, 
it was considered beneficial to characterize damping at extremely small displacement and strain 
levels. There has been discussion in the CSI community recently regarding dynamic behavior of 
structures at extremely low vibration levels. In particular, it was not known whether there was a 
radical change in properties below a particular vibration or displacement floor. In this study, 
damping was measured in aluminum and graphite/epoxy material specimens in air and in vacuum, 
and in the bare interferometer testbed truss. It was demonstrated that material damping was 
independent of strain from ten microstrain down to one nanostrain. Excellent correlation with 
thermoelastic material models was obtained. Damping in the testbed was found to be independent 
of strain below one microstrain (Figure 8). The linearity can be exploited by doing system 
identification at micron displacement levels instead of nanometer levels. The results were 
immediately instrumental in allowing the use of relatively inexpensive accelerometers for system 
identification on the testbed, rather than the extremely accurate high cost sensors. 
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Figure 8: Experimentally Determined Damping in 44 Hz Testbed Mode 
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Figure 9: Experimental Approach for Investigating Path Length Control 
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using a control strtegy that ignored the structural dynamics of the flexible base :snwnw. 25 ? 

approach was successful provided that the actuated mass was small ^P^ 

the structure. The effect of damping was investigated and quantified. An order of magnit 





reduction in vibration levels was demonstrated (Figure 10). This concept-implemented only for a 
single input single output case-will be extended to the interferometer testbed, where active mirror 
mounts will be used to position the cat's eye retroreflectors in three displacement degrees of 
freedom. Preliminary analysis of the finite element model suggests that the ratio of the moving 
mirror mass to the modal masses of the structure are small enough to allow the design of a high 
performance stable controller without further considerations of the structural dynamics. Research 
into active isolation will focus on the extension of this approach to cases involving noncollocation, 
multiple flexible modes, and multi input multi output systems. 



Motion (rms) 
(0-20.000 Hz) 

Aiieauanon 

Motion (peak) 
(0-20,000 Hz) 

Fixed Block 
Uncontrol led 

17 nm 


16 nm 

Controlled 

0.77 nm 

-27 dB 

5 nm (Si piezo mode 

Collocated 




Uncontrolled 

77 nm 


250 nm @ 1st mode 

Control led 

3. 1 nm 

-26.8dB 

11 nm @ piezo mode 

Non-Collocated 




Uncontrolled 

45 n m 


190 nm @ 1st mode 

Controlled 

3.9 nm 

-21.2 dB 

10 nm @ piezo mode 


Figure 10: Results from Path Length Control Experiment 


Finite Element Model 

The purpose of the finite element model is to provide a basis for analytical studies of 
structural modification, and to serve as one basis for control design. The accuracy of the finite 
element model is verified by comparison of frequencies and mode shapes with an experimentally 
derived modal model. It is not likely that the model will be used for control design if experimental 
models arc available. The effort in finite element modelling is outlined in Figure 1 1. 

Two finite element models have been constructed using ADINA: a continuum beam model 
and a model which contains separate elements for each strut. The continuum model has sufficient 
accuracy to make it useful for examining various approaches to control. Some features of the 
models arc described below. 

Continuum Model 

• Equivalent continuum cross-sectional properties for each leg of the truss were 

derived. (The six legs have identical cross-sections.) 

• Each of the six legs was then modeled with 14 Timoshenko beam elements. 

• The first flexible mode is at 38 Hz. 

• The low mode shapes are characterized by 1st and 2nd bending and the torsion of 

individual legs. 
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Figure 11: Near Term Efforts on Testbed Finite Element Model 


Full Struts Model . . 

• There are 228 nodes representing the aluminum joints. 

• Each of the 696 struts is modeled with a Timoshenko beam. 

• The first flexible mode is at 34.56 Hz. 

• The model runs in under 2 minutes on the Cray II. . 

• The low mode shapes are characterized by 1st and 2nd bending and torsion o 

individual legs. 

• There are 35 flexible modes below 200 Hz. 



Figure 12: Frequency Distribution of First 39 Finite Element Modes 










Distributio n of Modes 

Figure 12 shows the frequency vs. mode number for the first 39 modes based on the full 
struts model. Because of the inherent symmetry of the structure, there are repeated eigenvalues 
(multiplicity 2 or 3) present. Further, there is a separation between clumps of modes from 54-98 
Hz and from 142-194 Hz. The repeated roots and clumping of modes will disappear once 
concentrated mass of the three siderostats and the fourth vertex are added. Also, the added mass 
will drop the frequencies further so that there will be more than 50 modes below 200 Hz. 

Matlab Postprocessing 

The eigenvectors from the full struts model have been used to calculate strain energy 
distributions for each mode. The elements can be ranked from most to least strain energy by mode 
or sum of modes. This information will be used to choose passive damping element locations, and 
later as an initial criterion for active member location selection. With an improved model (including 
optics) we will be able to calculate a rough optical path performance metric to rank locations on an 
‘open loop’ controllability basis ( i.e . without simulating performance of the closed loop system). 

In separate work, a two-dimensional truss model has been used as a sample problem to 
develop necessary tools for control based on state-space models from the ADINA output. 
Implications of the close modal spacing and light damping are being studied. 

Addition of Damping 

The finite element model will be augmented with experimentally determined damping 
values. In addition, damping will be added in select elements in conjunction with the viscoelastic 
struts experiment. A more careful study optimizing passive damping locations will be carried out 
later. 

Role of the Finite Element Model 

The finite element model in its current form is a useful tool for parametric studies, mode 
shape visualization, calculation of strain energy distributions, and selection of system ED 
accelerometer locations. We have improved the accuracy of the model, but it is still not perfect. 
With initial ID data we will be able to make a direct comparison between finite element and 
experimental data. 

At some point a broader discussion of the role of finite element models in CST may be in 
order. There are several points which must be addressed. From the academic perspective at 
SERC, these include: 

• The ‘need’ to develop a highly accurate finite element model because it is standard 

practice in industry. 

• The value of a finite element model for laying out identification and control 

architectures 

• The inadequacy of finite element models as a basis for control in a complex lightly 

damped structure 

• The relative value of finite element models and experimental system identification 

• The realistic potential for on-orbit system identification 

• The need for an accurate finite element model if system identification is not 

possible 

• The role of a hybrid approach which could include subscale and component 

identification 



System Identification 

A system identification is mosTcomplex structure to be 

to control experiments. To date, the interfero re m ID will serve two purposes. 

for control design. - . h test5ed wi n not reach its ‘final’ configuration 

Because of the complexity of the °P^.“£*?, D m ^ performed in order to provide 
for some time (late fall 1990). However, an structure P and model will be sufficiently 

verification of the naked truss finite elem J'J t ^ difficult In addition, the initial ID will 

complicated later M that pacing sources of er^ will w M£ j hardware systems. The ID 

— -“ ,ers - Laler ’ acuve 

members installed in the truss will te ure resd^The first flexible mode was measured at 31.34 
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Figure 13: Near Term Efforts on Passive Damping Augmentation 


A passive-as-possible approach 
stringent shape or pointing requirements 


referred in achieving a spacecraft which must ^meet 
)wever, in an environment which will include several 


potential disturbance sources, some son of vibration alleviation will also be necessarv The 
introduction of an active control system can greatly improve performance at the expense of 
complexity, cost, and the possibility of instability. Passive damping augmentation is a far less 
glamorous, but nevertheless effective alternative. When an active structural control system is 
considered necessary, passive damping can only be beneficial. It does not make seise to 

™?i eiT f m ag f essive st ™ ctural control on a plant with only 0.05% inherent damping. Our initial 
goal is to conduct enough tests to establish the basis for later comprehensive experimental studies 
of passive damping schemes to be carried out in the future. The program has been broken into 
four areas representing different approaches to passive damping. These 8 are shown in Figure 13. 

Viscoelastic 

Constrained layer viscoelastic struts have been tested in a small cantilevered truss Poron 
and Scotchdamp materials were compared, and Scotchdamp was found to be more effective The 
effectiveness of different Viscoelastic layer thicknesses has been judged based on ringdown 
experiments in a first bending mode. A significant component of the strain energy of the structure 
is in the damper strut. This allowed high loss factors (25 %), and large drops in frequency 
, f . Twelve of simple highly effective struts have been manufactured. With information 

^?n KHw° nent T? er ’ f W1 1 56 ^ m? 1 ! to model the struts with equivalent axial stiffness and 
oredtt^ information will be integrated into the finite element model where a 

prediction of the added global damping due to several viscoelastic struts is possible. A repeat of 
the ID experiments will yield a measured value for damping. 

Shunted Piezoelectrics 

A resistively-shunted strut [4] was built and tested in a small cantilevered truss. The initial 

«f S thI S clH e y e dl uT ging ’ Wlth dampi " g values be,ow those expected. Subsequent re-engineering 
Sr.5 -f 1 ™ 1 yi u ded no lr yP ro yement. The use of a commercial Physik Instrumente actuator in the 
S f U if^ Ve n ° bener results. Although the experiment was designed to concentrate a large amount 
ot strain energy in the piezoceramic material, this was apparently not the case. A careful test of the 
strut in the component tester will provide an accurate accounting of strain energy distribution. The 
active member actuator stiffnesses were selected with consideration of appropriate stiffness 
properties for the shunting application in the large testbed. 

Concerns 

ctnitc /*oa\ dif (j cult P u°j fac ^ in adding significant global damping is the large number of 
struts (696) m the testbed. The damping members must be selectively placed, perhaps near critical 
payloads. At this point, the constrained layer Viscoelastic struts are by far the least expensive and 
easiest to make. Drawbacks include the frequency-dependent loss factor and material 
propertytemperature sensitivity. The shunted piezoceramics are potentially more effective than we 
have demonstrateci to date, but are expensive. A device based on the Honeywell D-Strut design 
which is capable of broadband viscous damping and is relatively temperature-independent is 
desirable, but is at this point prohibitively expensive to incorporate into the testbed. 


Component Tester 

, ... An axial component tester has been constructed and is operational on an optics bench. This 
tacility includes a Physik Instrumente piezoceramic strut to drive various test articles which 
represent subcomponents of the testbed. Mainly, these are passive or active replacements for the 
aluminum struts. Load and displacement are measured, the latter with a Zygo Axiom 2/20 
interferometer system. The tester will be used in the 0.1-200 Hz frequency range, with 
displacements from 1 nm to 60 pm. Initial measurements to be conducted are: 


• stiffness of truss longerons and diagonals 

• stiffness of active struts 

• voltage/deflection plots of active struts 

• viscoelastic strut characterization 



The facility will be available in the future for characterization of other passive or active 
components. 


Control Experiments 

We do not forsee having the capability to do absolute shape control in the near future, since 
that requires rigid body control of the testbed. The initial effort involves separation of the 
structural control and optical metrology path length control loops. Capabilities will be established 
in each through simple closed-loop experiments. Figure 14 shows the near term goals for control 
experiments. 
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Figure 14: Near Term Efforts on Control Experiments 


Model Basis for Control Design 

There are several methods for generating a model which is a suitable basis tor control 
design. These include finite element models (usually augmented with experimental damping 
values), measured models based on modal models from system identification, and measured 
models based on direct information from input-output actuator-sensor pairs. The third approach is 
preferred if the proper measurements can be made. For all these designs model reduction may be 
necessary in the plant and controller. 


Real Time Software , . , „ . ~ 

The software to do linear, constant coefficient, digital control is functionally complete, i ne 
code is called MatCon for matrix control. The user interface is through Matlab, where a typical 
continuous control design is discretized. The discrete matrices and some other constants (number 
of inputs, outputs, and states, scaling factors, and the sampling period) are saved in a standard 
Matlab .mat file. The real-time computer then reads this data and starts the controller. The 
following algorithm is used. 

input vector y from A/D 
x . = F..x. + F l2 y n 


' n +1 
U. 


11 

= F x 
21 J 


output vector u to D /A 
wait for next sample time 
























While the controller is running the user can stop and start the controller, record states, inputs, and 
outputs, and scale inputs and outputs. The data file of input vectors is stored in on-board memory, 
until a set number of samples has been saved. The controller is then stopped and the data 
transferred to the hard disk on a Sun Sparcstation, where it can be read back into Matlab. States 
and outputs can be reconstructed from the saved input data for full analysis. There is a direct 
interface to the six HP laser measurement boards. Four-pole Bessel anti-alias filters with a comer 
frequency set by digital input-output from the real-time computer are used. The filter cards also 
provide a digitally programmed gain of 1,2, 4, 8, or 16 to help amplify low-level sensor signals. 
We will have the capability to process 16 inputs, 10 outputs, 32 states at 1kHz. The control 
bandwidth is not expected to exceed 150 Hz. 

Active Struts 

The active strut design is shown in Figure 15. In addition to the load cell and internal strain 
gage measurements, two accelerometers are mounted to the strut to provide an inertial collocated 
measurement and to permit system identification using the active struts. Three struts are currently 
available, and an additional homemade unit will initially be used as a disturbance source generator. 
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Figure 15: Active Strut Configuration 




Active Mirr or Mounts 

Active mirror mounts will be used to maintain to X/20 the linear positions of the cat's eye 
retroreflectors, which are located at the three mock siderostat locations shown in Figure 3. Output 
position control will be achieved by moving the cat's eye and mounting table using three 
microactuators: 0.7" piezoelectric stacks for two of the active mirror mounts and 0.4 

electrostrictive stacks for the third, as shown in Figure 16. The actuators will be run in common 
mode to actuate piston, or z, motion of the point M of the cat’s eye. In differential mode, the cat s 
eye and table will be tilted; resulting in x and y displacements through the lever arm and flexure 
assembly. The rotations and lateral displacements cannot be controlled independently, but this 
constraint will not be a problem for the envisioned set of control experiments in the near to medium 
term. Simultaneous displacements of +/- 3.5 pm can be achieved in all three directions. The 
mirror mount design includes the flexibility to introduce additional mass to simulate the scaled 
mass of the retroreflectors. Additionally, the mounts can later be modified to incorporate mass 
reactuation, where the he effect of moving the mass of the cat's eye is reduced or even canceUed. 
The result will be a reduction in the interaction between the mirror control system and the truss 

flexibility. 
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Figure 16: Active Mirror Mount Functional Drawing 


The active mirror mount is a small stroke device intended to control only path length errors 
in the flexible truss. These errors will result from disturbance sources that are introduced 
intentionally to simulate space disturbances and from disturbances present in the ambient noise 
environment of the laboratory. Figure 17 shows the ambient acceleration power spectral density 
(PSD) in the worst-case direction measured by a triax of moderate-sensitivity acceleromete s 
(lV/g) at a proposed active mirror mount location on the truss. In this very preliminary study, 
structural dynamic response, starting at 30 Hz, is also corrupted by electrical noise and vmous 
lower frequency suspension modes. A displacement PSD is calculated by scaling the acceleration 
PSD by 1/cd 4 , which leads to an estimate of rms displacement of 22 nm in the frequency ban o 
20-100 Hz. Assuming that a point corresponding to a siderostat on another leg experience 
same disturbance and vibrates out of phase with the first point throughout this frequency range, an 
ambient path length error of 44 nm rms can be expected. The actual error may be less once 
electrical noise is removed and the additional mass of the active mirror mounts is added to the 


truss. 





In separate work [5], electrostrictive ceramic PMN:BA, a material of interest to structural 
control engineers, was characterized for test parameters of frequency, amplitude, and temperature. 
Results indicate that at room temperature the material strain response is quite linear with almost 
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Figure 17 Acceleration Power Spectral Density in Worst-case Direction at Proposed Active 

Mirror Mount Location 



0 10 20 30 40 50 60 70 80 90 100 

Hz 

Figure 18 Displacement Power Spectral Density in Worst-case Direction at Proposed Active 
Mirror Mount Location (Derived from Acceleration Measurement) 


no phase due to hysteresis, unlike piezoelectrics, and is constant with frequency. However, the 
induced strain sensitivity is highly dependent on temperature, and hysteresis increases rapidly 
below room temperature. Since electrostrictive actuators will be used in one of the three active 



mirror mounts, the operating temperature of the actuators will need to be monitored during 
calibration and usage. 


Summary 

The SERC interferometer CST testbed will soon be fully operational. The facility will 
address concerns regarding extremely tight constraints imposed on structural motion in future 
space observatories. At the same time, the testbed will serve as a platform for exploration o a 
broad range of controlled structure technologies and approaches. 

Work described in this paper was funded by NASA Grant NAGW-1335. 
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Abstract 

A cost functional i. proposed and inveatigated which U mot,- 
voted by minimising the energy in a structure using only collo- 
cated feedback. Defined for an X.-norm bounded system, hi. 
cost functional also overbound, the X, cost. So^e pr°Pert.« of 
this cost functional axe given, and preliminary mult* or » the pro- 
cedure for minimixing it are presented. The frequency 
functional is shown to have a time domain representation 
of a Stackelberg non-xero sum differential game. 

Introduction 

This paper examines the properties, evaluation algorithm and 
an optimization approach for a cost functional for combined 
Xs/W. control. Combined Xj and He control is of interest since 
it combine, the problems of nominal performance and robust sta- 
bility. Related work include. X, optimization with an X. con- 
straint (1-4), minimum entropy X. control [4,5], and mixed Xj 
and X. control [6, 7]. The cost functional of interest to us is 

defined a* followi. M „ nnA _ 

Definition 1 Consider a system H(s) ={Ho{*) Hi(s)] an 
number 7 £ 1R, with Ho € 77Xj, Hi € 7 ZHe, and ilffill. < 7- 
Then the cost I(H,7) « defined by 


L(H, 7) ~^J trtc *{(/ -7 _ ’HiH,*)' ''H 0 H;}d« 


(1) 


The specific form of this cost is motivated by minimizing the to- 
tal vibrational energy of a structure with only a model of the local 
dynamics near an actuator and collocated sensor Previous work 
with this type of model has u*d X, [8] and X. [9] optimization, 
of the power flow. Briefly, the fraction of the input power flow 
that is reflected into the structure at the actuator location is a 
quadratic at each frequency, and can be reprewnted by a transfer 
function HH‘. The fraction of the power that « dissipated is then 
(/ - HH'), and the total power dissipated is (/ - HH )t where 
E is the structural energy as a function of frequency. If the power 
flow into the structure from external disturbance sources 1. given 
by *(}w) then a power balance yields that the total energy in the 
structure is given by L([H H#],l). A more detailed explanation 
may be found in [10]. 

In [71 a framework for mixed X*/X» control problems is con- 
sidered. There the cost functional is motivated in an input/output 
sense. The system is subject to two inputs, one of bounded spec- 
trum, and the other with bounded power. For the case where 
the first input signal is white and the second is causal, necessary 
and sufficient conditions are given for the existence of a controller 
which minimizes the cost. The non-white and non-causal caw is 
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described but not solved (see Section 3.3.) This case is, however, 
of particular interest as the cost theu equals L{H, 7), revealing a 
close relationship between the present approach and the approach 
taken in 17], This connection is currently under investigation. 

The final section of this paper gives a third interpretation of 
this cost in terms of a Stackelberg non-zero sum differential game. 

Properties 

The following basic properties of L{H,i) will be stated without 
proof, and can be easily shown to hold. 

Proposition 2 Let H{s) and 7 satisfy the conditions in Dtfini- 
t\on l. Then 

(i) L(H,7) 13 we H defined. 

(ii) L(H,i) > 0, and L{H,i) = 0 <=> Ho = 0 

(iii) L{UHV,j) = L{H,i) for any U,V 6 1ZC m with U'U = 

/, VV = I. 

In the case where X, = H 0 , further properties of the cost 
L( H, 7) can be established by relating it to the entropy I(H, 7) ot 
s. system defined, for example, in Reference [5]. 

Definition 3 For H € 77X,. 7 6 1R, and ||H|, < 7, the en- 
tropy At infinity u defined 6y 

~ OO 

Also let C{H) be the usual X, cost associated with the system H, 

OO 

C(H) := — j trace (3) 

-OD 

ProDoeition 4 For H = [Ho H 0 ] , with H 0 and 7 satisfying the 
conditions in Definition 3, consider the cost L{H,l), the entropy 
I(Hq } 7), And the H j coat C{H 0 ). Define { - T , 

(i) l(H, 7) = £({'(#0,7)) 

(ii) 

Proof: The first assertion follow, directly from the proof of 

Proposition 2.3.2 in Reference [5). The first inequality in (11) fol- 
lows^om L{H, 7) = /(ft. 7) + $(/(Ho, 7)) “ d the *“ uU £r0m 
Proposition 2.3.2 in [5] that £ (/(Ho. 7)) * «• The 
ity is obtained from the result that /(H 0 , 7) itself bounds the X, 

C °That L( H, 7) overbounds an X, cost can also be shown to hold 
for the caae Hi ^ Ho- 


Propoiition 5 L{H, 7 ) ^ C{Hq). 



# 5 11 - < . 7l (/ - < 1 Md (/ - 

7 ,. < 1- The retult then follows directly from the defi- 

mtion of L(H, 7) m Equation (1). Q 

Finally, note that relaxing the *..norm bound completely re- 
covers tne 7 t 3 cost. 

Proposition 6 Km L(H,y) = C(H 0 ). 

L r0 °, f: m' f ? 11 ° W V direCtly from the definition of L(H,y) in 

Equation (1) and the Dominated Convergence Theorem. ^ □ 

Evaluation of the Cost 

uSbIT^X] 5P<lCe repre,eaUtion for 4 5tnctl V Pr°P« sys- 



H = 


A | (?o B, 

IT 0 0 


= C(s/- A)-' [Bo B x 


W 


Figure 2: Feedback System 


Optimization 

fnr T . h ,' g °*}. of thi * >» Present *n spprosch for tolvin 

he &T controlJer th4 ‘ minimizes * cost functionil of" 
ne lorm (1). Linear time- invariant controller* w :n k , 

^.u,. by, thi, f, ra , aot 

r” “»*““• tk “ “opting » 

P 2* £**'“ C4n be described by the block diagram in Figure 2 
^ csn be realized in state space as 6 8 ' i 


P' X” eV4lUatC L[H ' l] iD teraU ° f the 5Ute 5 P 4ce data. 

.tnrtll h n 4 D0n ' zer ° tenn ,.f> couJd be included. H x is made 
strictly proper only to simplify the results. 

•^llTlll <\ H Thti^ ^ 6t *' V ‘ n h EqUat '° n Uh 7 € JR - 

= trace {C^C 7 } (5) 

«*e« P, Q satisfy ( A + -,-'B x BjP) stable and 

PA + A t P + -,- 3 PB x B?P + C t C = 0 (6) 

(A + 7 r'B x Bj P)Q + <?(X + ■r'B.BjPf + B o5 7 = 0 (7) 

Proof: Since < 7 ’/ Vw, then 3Af*’ 6 glV en by 

M-M = B*(/ - ‘r~ i H x H’y i H 0 (8) 

ri,.', 14 ^- r“ L repre *“ Uti ° n for c “ be found by noting 

mFi^Lfl “gf* tr “ ,f " functlon of th « feedbac k system show£ 


P = 


' A 

Bo 

B x 

Bs 1 

,c 7 \ 

0 

D so 

0 

D i, 

~TKT 

0 


(11’ 



A 

T*B x Bf 

B 0 ' 

M'M = 

-C T C 

-A T 

0 


L o 

B I 

0 


- k-iausiormanon, 

H-P(P.K), I *,»/>„] + P„K{I - P„K)-‘ )P K p„] (, 2 . 

W rT,° Tk™ V‘“ ,h0 " »“* “*«>“ P. -cl 

satisfy < 7 . The problem statement is then 

min{£(P't 7 ) : K admissible} ( 13 ) 

By^a scaling of H, without loss of generality consider the cast 

isJ hC ^ I 1 * 1 ' Problem is examined first, with normal- 

ised control weighting, so that Cj = [C 7 0) and Df t = (0 /] 

Theorem 8 Consider the problem statement (IS), with Ac l — 
then:^ ^ “ * ***** f"<&ack matrix that solves Equation (IS), 

F = ~ B I(PQ + PQ)(Q + Q)~ l (i4[ 

P, Q, p t and q satufy ^ = {Acl + BiB r p) ataiu ^ 

(is; 

(is: 

(17; 
(18 


0) 


a J - f 

PAcl + Al L P + PB x BjP + C T C + F t F = o 
A^Q + QA + BoBf = 0 
+ aL,p + C T C + f t f = 0 

^<5 + dAl^ + QPB x BJ + B x BjPQ = 0 


Proof: The closed loop system it H = 


With P given by Equation (6), then 


Act 

[Bo 

Bx I 

(? 

0 

0 - 

F 

0 

0 


Fror 


iy-p+7-*BiA r P|a.j 


( 10 ) 

;:^n 4b ! ef4C f 0r ^ l»J- Substituting Equation (8) 

. ® JJJ* “ " cle4r from (*) th4t ^e cost L(H,j) i, then given 
»»y ||J»/||„ where ||Af|| 3 = trace {CQC*} and Q satisfies the Ly* 
punov equation (7) [13]. } y 


- A)-* | ~ j(-s; - ^)-c T l — fjf]_ 




Figure 1. Block Diagram for M'M 


,h T ,, 1 ! ere( 

l“Jn h fl (16)> “ d P ,4tiifie# ‘be Ricca 

equation (15)^ Appending these two equation, to the cost a. co. 

stramta with Lagrange multipliers P and 6 respectively yield. rK 

aT ,o ' p ' iMF «<»• 

Preliminary result, indicate that an iterative approach to sol’ 

F 8 Gi^. « U4t, i n *| C ° nVer ? e ' r * pidly to the °P‘ im * 1 ^back la 
F. Given an uu ,a^ guess for F («y, from the minimum entror 

«ntml problem (5|), P,Q, B, end <} can be computed ^uenti J 
the wlution of Riecati and Lyapunov equations. Equation (1 
c*n then be evaluated for F, and the ptock. repeated 

dJ°rL Yhe ,tem ° f ° rder 14 “ d * &ted order w^^ensator of c 
der n., the neceasary conditions for the optimum can be foun 

M n using a Lagrange multiplier approach, in terms of 4 ore 

+ n e mAtruc equation* iimilAr to Equation* (15W18) Work 

cunotly in progress to simplify and interpret thesTresults for t 

dynamic compensation problem. Note that there is no a pn. 

£“V° * 3tp * et th ** no improvement in the coet can be acluev 
ior rig > ft. 
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Time Domain Interpretation 

The form of the augmented co*t for the linear control problem 
lead* to an interesting differential game* interpretation. It is well 
known that the central controller in the problem can be found 
as the solution to a zero-sum differential game [14], where for min- 
imizing IIT^H,., the control ti and noise w solve the optimization 
problems: 

u = argmin [ z T z - 'y 7 w T w dt (19) 

Jo 

w = argmin / — + ^w^uj dt (20) 

Jo 

u has some information y about the state, and u; has full informa- 
tion. . 

With the current cost functional, and under the assumption ot 
linear feedback, the optimization problem is again equivalent to a 
differential game, but it is no longer a zero-sum game. Whether 
the two problems are equivalent when both are allowed nonlinear 
feedback is unknown. 

Proposition 9 // an optimal linear compensator exists for prob- 
lem (13), then it is the same as that of a Stackelberg differential 
game with u as leader, as follower and wo as unit intensity 
white noise , where u and ti* solve the following optimization prob- 
lems: 

u = argmin Um E (^) 

u/i = argmin lim E |-z T z + Vwftui} (22) 

u has some information y about the state , and w has full informa- 
tion. 

Proof: Assuming a linear control law for u, the optimization 

problem for w x is easily solved with a single Riccati equation 
(which is Equation (15) for the state feedback case.) Append- 
ing this as a constraint for the optimization problem (21) results 
in an identical problem formulation to that of problem (13). ^ 

This game seems to be a more natural problem to pose than 
the pure H m differential game, since the control does not benefit 
from the use of noise, but instead optimizes an Hi type of cost 
functional, while the deterministic noise w x solves the some opti- 
mization problem as before. In addition, the plant is subject to a 
white noise input u? 0 . This looks similar to the framework of [6,7] 
since a single output is minimized in the presence of two distur- 
bance inputs, one of which is associated with the Hi nature of the 
problem while the other is associated with the H m Q*ture. 

Note that for a non-zero sum differential game, the solution 
depends on how the optimality is defined. For the Stackelberg 
or leader-follower solution [15-17], one player (here the control u) 
acts as leader and announces a strategy, and knowing this strategy, 
the follower (here the noise w x ) solves its optimization problem. 
Also note that in general, the optimal control for the Stackelberg 
problem is known to be nonlinear (17). Similar equations to (15)- 
(18) have been reported in [16], where the optimal linear state 
feedback law for a Stackelberg problem was found. The nonlinear, 
team optimal strategy obtained, for example, in [17] does not ap- 
ply to this problem since the leader u cannot increase the follower 
uVs cost indefinitely, and therefore cannot induce w x to follow a 
strategy desirable to u. 

The differential games representation of problem (13) allows 
the matrices of Equations (15H 18 ) to be P vca “ interpretation. 
— P and +P correspond to the optimal cost-to-go for the costs 
associated with w x and u respectively, Q is the co variance of the 
state, and is the sensitivity of the cost for u to changes in the 
cost for w\. 
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Abstract 


In certain applications modeling uncertainty can be represented by a finite number 
of plant models. This paper considers the problems of determining a feedback con- 
troller or estimator that optimizes an 7i 2 performance criterion involving a collection 
of plant models. The approach is based upon fixed-structure optimization in which 
the estimator or controller order are fixed prior to the development of optimality 
conditions. 

•% 

1 Introduction 

The goal of robust control design is to obtain controllers that maintain desirable per- 
formance in the face of modeling uncertainty. In certain cases modeling un certainty 
can be adequately represented by means of a finite number of plant models. This 
multi-model problem arises for example, if the plant can undergo sensor or actuator 
failure modes. A finite set of models has also been used to design for robustness to 
an infinite set of models, as in the case of parametric uncertainty [1], high frequency 
uncertainty [2], or parameter variations (e.g. for different flight regimes) [3]. 

A fundamental issue in multi-model problems is the simultaneous or Tellable sta- 
bilization problem. Here the goal is to design controllers that stabilize each model in 
a finite collection of plant models. Considerable progress has been made in solving 
this problem [4-10]. 

The goal of the present paper is to consider a multi-model optimization problem. 
Specifically, we consider a quadratic {'H 2 ) performance criterion involving a collection 
of plant models controlled by a single feedback compensator. The approach we take 
involves fixing the order of the compensator and optimizing over the feedback gains. 
This approach is similar to that of [11] where static output feedback controllers were 
considered. 

One of our principal objectives in considering the Multi-Model Control Problem 
is to examine the issue of compensator order. In [6] it is shown that simultaneous 


1 



stabilization of a pair of plants of bounded degree may require a compensator of 
arbitrarily high order. In the present paper we show how this issue manifests itself 
in the structure of the necessary conditions for optimality. 

To further elucidate the role of compensator order we also consider two related 
problems that are simpler in structure but that involve analogous issues. The ob- 
jective of the Multi-Model Approximation Problem is to determine a single model 
that simultaneously approximates a finite collection of models. For a collection of r 
models each of order n { , t = 1, . . . ,r, the maximal-order solution is given by a model 
of order £ n<, which is larger than the order of each of the given plant models. 

In a related vein we also consider the Multi-Model Estimation Problem wherein 
we seek an estimator for each model in a given collection of plant models. As in the 
Multi-Model Approximation Problem the maximal-order solution has order greater 
than the individual plant models. 

The fixed-structure approach applied to the multi-model problems is a direct 
extension of the technique utilized in [12-14], Indeed, by specializing these results to 
the case of a single model, the results of [12-14] are immediately recovered. 


2 Multi-Model Approximation 


Consider the following problem. 


Problem 1 (Optimal Multi-Model Approximation Problem) Given a set of 
r controllable and observable systems Hi, t = 1 . . .r, with state space representations 


* Ai 

Bi 

Ci 

0 


Ci(sl - AiY'Bi 


( 1 ) 


and a set ofr numbers ca € H, a, > 0, i = 1 . . . r, find a single approximation model 
of fixed order n m , with state space representation 


H m = 


A m 

B m 

C m 

0 


( 2 ) 


2 


that minimizes the weighted 7Y 2 model- approximation criterion > 

= (3) 

l-l 

To guarantee finite cost J, assume that each Hi is stable, and also restrict the 
optimization to the set of stable approximation models H m . Furthermore, since the 
value of J is independent of the internal realization of H m , assume that the realization 
in Equation (2) is controllable and observable. Thus require that c m ) e n 

where 

^ — {( > -^m ) ) * Am stO/hlc^ (j4m> ^m) COTliro/lflWc, (Amy Cm) obsCWClblc } 

Without loss of generality, the weightings a* can be assumed to be normalized so that 

2°^ = 1 
1=1 

With this normalization, the weighting can then be associated with the probability 
that H{ accurately models the dynamics of the system. 

The necessary conditions for an optimal solution to this problem are given in 
Theorem 4. The approach used to obtain them is presented briefly here for comparison 
with the approach required for the multi-model estimation and control problems in 
Sections 3 and 4, and for the single model problems in [12-14]. 

The model approximation error transfer function Hi — H m can be represented in 
state space as 



The cost J is then 

J{H m ) = ^tr latiCiQiC?} (5) 

»=i 1 ' 

where each Qi satisfies a Lyapunov equation 

AiQi + QiAJ + BiBj = 0 ( 6 ) 
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Appending these constraints to the cost with Lagrange multipliers Pi yields first order 
necessary conditions for a solution upon differentiation with respect to Qi, A m > > 

and C m - Each matrix and Qi has dimension (n; + n m )x(rii + n m ) and can be 
partitioned into n^xn^, riixn m , n m xrij and Timxrim blocks as 


1 

h, ' 

«o. 

m. 

II 

1 

1 

Q*12 

l Pin 

A, J 



Q*77 


(7) 


The necessary conditions are then Equations (6), and 


H- = PiAi + AJ Pi + C, T C. = 0 Vi 

dQ i 


dJ 


dB 


= y + Y, = 


o 


»» i=i 


dJ 

dC m 

dJ 


»=i 

r 


dA 


= C m Yj a iQi ? a — S — ® 

l=l »=1 


( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 


"» i=l 


The equation obtained from differentiation with respect to A m is of particular im- 
portance in simplifying and understanding the structure of the necessary conditions. 
In the case of a single model (r = 1), Equation (11) yields 


/nm = (-P 2 VP2l\(QllQnl ( 12 ) 

r g t 


A projection operator t = G^T = t 2 is then used to simplify the equations [12]. 

For the multi-model approximation case, from Equation (6) each Q% 33 satisfies an 


identical equation, 

AmQijj + QinAm + B m B m — 0 

(13) 

Hence the Q ij3 satisfy Q t33 

= Q 22 , i = l---r. Similarly, each P in satisfies 



P i „A m + A^P l „+C^C m = 0 

(14) 
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and hence = P 2 2 , * = 1 • • r. Furthermore, from Equations (13) and (14), Q 2 2 and 
P 22 are the controllability and observability grammians, respectively, of the system 
H m . With these simplifications Equation (11) can be written as 

r 

P 22 Q 22 + 51 Qi ia ) = 0 (15) 

1=1 

This immediately gives the following result: 


Proposition 2 Given a fixed order model that is optimal for Problem 1 , of order 

r 

Km > N — £ Hi, there exists a model of order N with the same cost. Hence with no 

i=i 

fixed order constraint , the optimal system for multi-model approximation has order 
no larger than N . 

Proof: Vi, rank < n t , hence rank j £ oiiP lai Q {l , j < N. So, from 

Equation (15), rank {P 22 Q 22 } < N. If n m > N, either P 22 or Q 22 or both must not 
be full rank, and thus the representation of H m must have states which are either 
uncontrollable or unobservable. (The maximum number of states which are both 
observable and controllable is N . ) Removing any uncontrollable or unobservable 
states yields a system with identical cost and at most N states. □ 

With the controllability and observability assumptions on the representation of 
P 7 2 Q 22 must be positive definite, and thus Equation (15) can be written 




(16) 


1=1 


r. 


GfT 


Define 


r 

G 




r i ••• r r 

[ ••• Gr 

T = G t T 


(17) 

(18) 


t± = In ’ - t 
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(19) 

( 20 ) 



Then r is again an oblique projection operator, that is r 2 = r. Note that in general, 
r is oblique rather than orthogonal, since it need not be symmetric. 

The following lemma from [12] is required for the statement of the main theorem. 

Lemma 3 Suppose Q, P € IR Wx * a™ positive semi-definite. Then QP is nonneg- 
ative semisimple (has non-negative eigenvalues). Furthermore, i/rank{<?P} - n m , 
then there exist G,T € IR w "‘ xW and positive semisimple M € W"**" such that 

qp = g t mt ( 21 ) 


Matrices G, I\ and M satisfying the conditions of the lemma will be referred to as a 


projective factorization of QP. 

It will be convenient to compile the state space information about all of the models 
into a single set of matrices (A, B , C a ), where 



Ai 0 • • • 0 


B 1 ' 

A = 

0 a 2 

B = 

b 3 


0 At 


Br 



QlCi 



(23) 


The subscript a on C a indicates that it depends on o^. 


Theorem 4 Suppose ( A m ,B m ,C m ) solves the optimal multi-model approximation 
problem (1). Then there exist positive semi-definite matrices Q, P € IR such 
that, for some projective factorization of QP , A m , B m , and C m are given by 

A m = TAG T 
B m = TB 
Cm = CaG T 
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(24) 

(25) 

(26) 


and such that the following conditions are satisfied: 


rank {£} = rank {/>} = rank {<?p} = n T 

AQ + QA t + BB t - t ± BB t t 7 = 0 
PA + A T P + C T C - rlClC a r x = 0 


(27) 

(28) 
(29) 


Proof: Define Q = &Q 22 G and P = T t P 22 T, and note that tQ = Q, and Pr = P. 
Pre- and post-multiplying the Lyapunov equations (13) and (14) for Q 22 and P 22 by 
either 7„ m = TG 7 or 7 nm = GY T yields the following equations: 


AQ + QA t + BB t J = 0 
PA + A T P + ClC a ] T = o 


(30) 

(31) 


The (1, 2) sub-blocks of the Lyapunov equations (6) and (8) yield identical equations. 
Equation (30) is equivalent to Equation (28) since Equation (28)=(30)-|-(30) r -(30)T, 
and Equation (30)=r(28). Similarly, Equations (31) and (29) are equivalent. Note 
that only two Lyapunov equations are required for the necessary conditions because 
the (1, 1) sub-blocks of both Equation (6) and Equation (8) are superfluous. 

Equations (25) and (26) follow directly from (9) and (10). Equation (24) for A, m 
is obtained from the (2,2) block of either (£ o,(Eq’n 8)&) or (£ o^Eq’n 6))), 
either of which yield that £ = 0. D 

Because the form of the equations is identical to that of the single model case, 
the discussion in [12] applies for this problem as well. As in [12], the form is a result 
of optimality, and not fixed beforehand. If (A m , B m ,C m ) satisfies the necessary con- 
ditions, so does (rA m 7 , - 1 ,rB m , CmT' 1 ) for an arbitrary nonsingular transformation 
matrix T. Further, there exists a similarity transformation which diagonalizes QP 
and r simultaneously. Representing r in terms of QP as in [12] leads to numerical 
algorithms for the optimal multi-model approximation problem. 
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Remark 5 In the “fall order ” case n m = N, then t = G = T=I Ni giving A* - A, 
B m = B, and C m = C a . Thus H m = £ c^H,. This is exactly the expected result; the 

best possible approximation is simply the weighted average of all the models. 

Remark 6 For a single model (r = 1), the equations clearly collapse to the equations 

of [12]. 


3 Multi-Model Estimation 


Consider the following problem. 

Problem 7 (Optimal Multi-Model Estimation Problem) Given a set ofr sys- 
tems H i = \ ...r, with state space representations 


H, = 


A 

B t 


c 


0 

C h 

C it 

m 

0 

Di 


. Ci, 

(si — -A*) 1 Bi + 

Di 


(32) 


and a set ofr numbers W € 1R, <* > 0, i = 1 . • r, find a single estimator of fired 
order n e , with state space representation 


H e = 


' A e 

— » 

C t 

0 


(33) 


that minimizes the weighted Tii model- estimation criterion, 

= < 34) 

1=1 

where Hi is partitioned into H h and H tt according to the two outputs. 

The estimation problem can be illustrated by the block diagram as shown in 
Figure 1. 

The following assumptions about the problem will be made: 


8 



Figure 1. Estimation problem for each system. 

(i) Each H t is stable, and each (>l t , C* a ) is detectable. 

00 E = 1 

*=i 

(iii) For clarity in understanding the form of the equations, the process and mea- 
surement noise for each model will be assumed to be uncorrelated, so B x Dj = 0. 
Without loss of generality, assume Bi = [£. o] and D { = [o £>•] 

(iv) Require thaUhe measurement noise have no singular directions common to all 
models, so £ a x D t Dj > 0. This is a generalization of the usual single model 
assumption of nonsingular measurement noise, DD T > 0. 

It is interesting to note that each D t need not have full row rank, hence the 
estimation problem for each individual model may be singular without the multi- 
model problem being singular. 

As in the multi-model approximation case, require (A e , B t , C ,) 6 71. 

The model estimation error transfer function H tl - H.H it can be represented in 
state space as 



The cost J can again be written in the form of Equation (5), 

W) = Etr|o.C,4(7 l T J 
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where each Q , satisfies a Lyapunov equation identical to Equation (6)). 

Necessary conditions can again be obtained using a Lagrange multiplier approach. 
The Lyapunov equations for P{ are identical to Equations (8). The equation obtained 
by differentiating with respect to A e is the same as Equation (11), and once again this 
will be the key equation for understanding the structure of the necessary conditions. 
For this problem, each Qi 3J satisfies 

A e Q i2i + Q t27 A T e + B € Ci 7 Qi x2 + Qi 7l C?Bj + B e Bj = 0 (37) 

Each Qi 22 now satisfies a distinct equation, and thus Qi 22 ^ Qnu * ^ b The critical 
observation for this problem, however, is that each Pi 22 still satisfies Equation (14). 
Thus it is still true that P t „ = P 22 , i = 1 . . .r. This is sufficient to obtain the 
elements of a projection operator from Equation (11), and to prove the following 
result, analogous to Proposition 2. 

Proposition 8 Given a fixed order model that is optimal for Problem 7, of order 

n m > N = 53 there exists a model of order N with the same cost. Hence with 
1=1 

no fixed order constraint, the optimal system for multi-model estimation has order no 
larger than N. 

Proof: As in the Multi-Model Approximation case, rank j £ OiAj , Qua } ^ N- 

From Equation (11), rank|p 22 E <*&„} < N. If n ra > N, either P 32 or £<*<?,„ 

or both must not be full rank. P 22 is the observability grammian of the system 

(AcSe,^), and thus is not full rank if and only if (A e ,C 9 ) is unobservable. Also, 

E <*iQin is not toll ranlc and onl y ^ £«) is not controllable. This result will 

i=l 

be proven in Proposition 13. Proposition 8 then follows in the same manner as the 
proof of Proposition 2. a 

Remark 9 The estimator must obtain all the information possible about the state 
from the output y. Since all state information from all the models has a finite di- 
mension N , there is an estimator state vector of dimension N that contains the most 
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information possible about the state vectors of the Hi. Any additional estimator states 
must be redundant. 


As noted earlier, P 2 2 is the observability grammian for H t and therefore must be 
positive definite. Proposition 13 proves that 53 must also be positive definite. 

»=i 

Hence for the multi-model estimation problem, Equation (11) can be written as 



,=1 T. — 

or 


(38) 


With G , r, and r defined as in Equations (17—19), r is again a projection operator, 
satisfying r 2 = r. 

In addition to the definitions of A and B, given in Equation (23), this problem 
requires C*> and C a j, defined analogously to C a , and 


BiBf 0 ••• 0 

0 B 2 BJ 

o B r Bj 


(39) 


V 2 = Y.^D x Dj (40) 

1=1 

Theorem 10 Suppose ( A e , B e , C e ) solves the optimal multi-model estimation problem 
(7). Then there exist positive semi-definite matrices Q, Q, P € JR NxN such that , for 
some projective factorization of QP, A e> B e , and C t are given by 


A e = rAG T -B e C a ,G T 
Be = TQClVf' 

Ce = C ai G T 

and such that the following conditions are satisfied: 


(41) 

(42) 

(43) 


AQ + QA t + V,~ QCl t V 2 ~'C ai Q + t^QCI^-'C^Qt^ = 0 (44) 
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AQ + QA r + QCl 7 V 2 ~ l C ai Q - r^QCl^C^Qr^ = 0 


(45) 


P(A - QCl Vf l C a , ) + ( A - QCl 7 V 2 - l C ai ) T P + Cj t C ai - t l C J C ai r x = 0 (46) 


Proof: The derivation of these equations is similar to that for the necessary con- 

ditions for the multi- model approximation problem. Define Q = G 7 G, 

P = T t P 33 T } and Q = diagja-^^} - Q. Substituting into the Lyapunov equa- 
tions defining Qi and Pi yields Equations (45) and (46) from both the (1,2) and 
(2,2) sub-blocks. The (1,1) sub-block of the Qi Lyapunov equation can be used 
to obtain Equation (44), and the (1,1) sub-block of the P, equation is superfluous. 
Equations (42) and (43) follow directly from the equations obtained by differentiating 
the augmented cost with respect to B e and C t . Equation (41) for A« is obtained in 
an analogous fashion to the approximation problem. D 

As in the multi-model approximation case, the necessary conditions obtained here 
are similar in form to those for the single model case [13]. Again, the necessary 
conditions hold for any non-singular state transformation of the estimator. Numerical 
algorithms developed for solving the equations in [13] can be applied to this problem 
as well. 


Remark 11 In the “full order ” case n e = N, then r = G = T = In, giving A* = 
A- B e C a]l B e = QClV 2 l , and C € = C ai . Only the Riccati equation for Q needs to 
be solved, and this has the same form as the Kalman filter equation. Because of the 
coupling of the multiple models in Q, the full order estimator is not simply a weighted 
average of the individual model estimators. 

Remark 12 For a single model (r = 1), the equations clearly collapse to the equa- 
tions of [12]. For r = 1 and n t = N, the equations collapse to the standard Kalman 

filter result. 

Finally, the proposition used in the proof of Proposition 8 needs to be proven. 
Proposition 13 Q = £ <*iQin » fall rank if and only if (A*, B e ) is controllable. 
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Proof: Q satisfies the Lyapunov equation 


(A. + B t C a3 G T )Q + Q{A e + B e C a ,G T ) T + B e V 3 Bj = 0 (47) 

This follows from summing Equations (37) and representing Qi lt in terms of G 
and Q. Q is therefore a controllability grammian, and is full rank if and only if 
( A e + BtCajG 7 , B e ) is controllable. This system is controllable if and only if ( A t , B e ) 
is controllable. □ 


4 Multi-Model Control 


A simple form of the necessary conditions for the multi-model control problem is sig- 
nificantly harder to obtain than for either of the two previously considered problems. 
A form of the equations similar to the single model case has not yet been obtained. 
The problem will be set up here, and the critical issues discussed. In particular, the 
question of controller order is investigated. 

Consider the following problem. 


Problem 14 (Optimal Multi-Model Control Problem) Given a set of r sys- 
tems H l} i = 1 . . .r, with state space representations 



Ai 

B n 

B» 

H, = 


0 

Bin 


_c l3 

A, t 

0 


(48) 


and a set of r numbers ati € IR, > 0, » = 1 . . . r, find a single compensator of fixed 
order n ci with state space representation 


H e = 


u 

B c ' 

Oc 

0 


that minimizes the weighted Tij model-control criterion, 


(49) 


j(H e ) = £> Usui* 

«=1 


(50) 
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Hi is partitioned into H ni , H , l3 , B,,, and Hi n according to the two inputs and two 
outputs. The closed loop transfer function Hi tu is obtained from the lower linear 
fractional transformation, Hi tv = F(Hi,H c ). 

The control problem can be illustrated by the block diagram as shown in Figure 2. 


w 


u I 


Hi 


H e 


Figure 2: Control problem for each system. 


The following assumptions about the problem will be made. 

(i) Each system Hi must satisfy (A^B,,) stabilizable and ( Ai,Ci 3 ) detectable. 

(ii) £ a, = 1 


»=i 


(iii) For a compensator H e to exist which gives finite cost J, the set of systems Hi 
must be simultaneously stabilizable. Conditions for simultaneous stabilization 
have been studied by Ghosh and Byrnes [6]. 

(iv) As in the estimation problem, assume uncorrelated process and measurement 

noise, so B^ DJ 31 = 0. Without loss of generality, again take B tl = [B^ 0] and 
Di tl = [0 • Further, require that the measurement noise have no singular 

directions common to all the models, so £ a.Bj,, > 0. 


«=i 


(v) The dual assumptions to (iv) will also be made. That is, Cf x Di lt = 0, Ci x = 


» Bi lt — 


D, 


, and £ o^Di Xi Dj > 0. Note that for any individual 


Ok 
0 

Hi, the control weighting Di l3 Df may be singular. 


i=i 





The optimization will be restricted to the (non-empty) set of simultaneously sta- 
bilizing compensators H c , with controllable and observable realizations. 

The closed loop transfer function can be represented in state space as 



The cost J can again be written in the form of Equation (5), 

W) = I>{<*C,Q.6 T } (53) 

t=l 1 J 

and again, each Qi satisfies a Lyapunov equation identical to Equation (6). 

Necessary conditions can again be obtained using a Lagrange multiplier approach. 
The Lyapunov equations for P x are identical to Equations (8)* Once again, the equa- 
tion obtained by differentiating with respect to A c is the same as Equation (11). 
However, for the control problem, there is a crucial difference. Each Qi 27 and Pm 
satisfy, respectively, 

+ Q,„A* + B,C„Q n , + Q in C?Bj + B'BJ = 0 (54) 

ft,. A: + A r c p.„ + Pi., B h C c + CfBjPi,, + CjC t = 0 (55) 

Thus for this problem, every Q tJ3 and every Pm is different, that is Q til ^ Q^, i ^ j, 
and P Xn P n3 , t ^ j. As a result, Equation (11) is difficult to factor, and this also 
has serious implications on the order of the compensator. 

Proposition 15 There is no a priori bound on the order of a compensator which is 
optimal for Problem 14- 

Proof: Ghosh and Byrnes [7] give an example of two second order systems, param- 

eterized by A, which require an arbitrarily high order compensator for simultaneous 
stabilization as A tends to some limit. Since any optimal compensator must be si- 
multaneously stabilizing, it also may be of arbitrarily high order. D 
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Remark 16 The result in Proposition 15 has been shown before; the purpose of re- 
stating it here is to illustrate how the result manifests itself in the present context. 

r . 

For all three of the problems investigated in this paper, 53 mos ^ 

«=i 

rank N. Equation (11) then yields that £ OiP t3i Qi 33 has rank less than or equal to 

t=l 

N. For controllable and observable systems H m , H e , and H c , each term in this last 
sum has rank rim, n e , or n c . In the approximation case, this sum can be factored as 

r 

P 22 Q 711 and i n the estimation case, it can be factored as P 22 E ^ Q*n- Sylvester’s 
inequality [15] can then be used to show that this second sum has rank equal to n m 
or n e . From this, the conclusion that n m < N, and n e < N follows. In the multi- 
model control problem, the sum £) a,P, 33 Q, 33 may have maximum rank N while the 

t = l 

individual terms in the sum can have larger rank n e . That is, the optimal compensator 
may be both observable and controllable for arbitrarily large order n c . 

Theorem 17 Suppose ( A c ,B c ,C c ) solves the optimal multi-model control problem 
(14)- Then there exist positive semi-definitt matrices Qi, P 3 6 ]R( n ' +n ') x ( n * +n *) such 
that Ac, B c , and C c are the solutions of 

j2ai(P, i3 A c Q l]3 + P i3l AiQ il3 + Pi 31 B c Ci 3 Qi l3 + P i3l B i3 C c Q i33 ) = 0 ( 56 ) 

»=1 

'£a i (p i „B,D iu Dl l +(P i „Qi. l +Pi„Qi, x )C^) = 0 ( 57 ) 

«=1 

(Dj, A.,C«<3 i„ + + K.Q.J) = 0 (58) 

»=1 

where Qi and P 3 satisfy Equations (6) and (8) respectively, with the appropriate par- 
titioning given by Equation (7). 

Proof: Equations (57) and (58) are the necessary conditions obtained directly 

from differentiating the augmented cost with respect to B e and C e • Equation (56) for 
Ac is obtained in the same manner as for the approximation and estimation problems. 
□ 
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Remark 18 Equation (56) can be solved for A* using Kronecker algebra [16]; 
vec {.A e } = x 

£ a * ( (Qu, ® A S1 ) vec {A,} + vec {P,„ B e C ia Q il7 + P ijt C c <?,„ }) (59) 

1=1 

Note that the inverse in Equation (59) exists. To see this, note that each Pj„ and Qi 17 
are positive definite, and their Kronecker product is therefore positive definite [16]. 
The sum of these products is therefore nonsingular. 

Remark 19 If D %ai = piDa, i = 1 . . .r (which may not be an unreasonable assump- 
tion,) then Equation (57) can be written as 

B c = (£<**„) + P,„Q.»)cz) Vf' = 0 (60) 

where 



In general, B c can be solved using Kronecker algebra. Similar comments apply to the 
calculation of C c . 

5 Conclusions 

The simultaneous optimal approximation, estimation and control problems for mul- 
tiple models has been investigated. In each problem, the order of the system to be 
found is fixed, and the necessary conditions that an optimal solution must satisfy are 
found. For both the approximation and estimation cases, the optimal model can be 
written as an optimal projection of a “full order” model with order N = £ n i- There 

»=i 

is no improvement in the optimal cost that can be obtained by using a model with 
order larger than N . In the control case, there is no such a priori bound in terms of 
the individual model orders n, that can be placed on the optimal compensator order. 
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ABSTRACT 

The exact Linear Quadratic ,^ 0 ^ 

the convolution of spatially distributed ^““^^e funct.ons exist that 
continuous state functions. Fo r s , nK . tu ' re a t any given point in time. The 
completely describe the state of the convolution! of one of these state 

upon the ease with which it can be measured. 

This paper discusses the estimation of exact 

rate feedback kernels from to equivalent 

element structural models. Th kernels These curvature kernels are 

curvature and curvature rate feedback ker . q feedback to account for 

augmented with single point displacemei -.tvature rate state functions can be 
rigid body motions. The curvature and ^tme rate sta sens0 rs 

measured using a growing class of sensors “0 convolution of all structural 

Curvature' sites" wXThe Tern's 

an^usin^ are Taveraghig sensors enables the implementation of full state 
feedback for infinite order structural systems. 
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Introduction 


Rationale 

Structures are infinite order systems. To obtain a mathematically exact 
structural model requires the use of a set of partial differential equations 
subjected to the appropriate boundary conditions. However, in practice it is 
difficult or impossible to find the exact closed-form Linear Quadratic 
Regulator^ (LQR) solution for most structures. Therefore, structures are often 
modelled by discretization of the structure. This is even true for some very 
simple structures. The most common method of discretization is finite 
elements When the stnicture is discretized, the order of the model is reduced. 
Instead of being modelled by an infinite order system, the structure is now 
modelled with a finite number of degrees-of-freedom. The result is a matrix 
ordinary differential equation which will approximate the temporal and 
spatial behavior of the structure. 

Given the possibility to model a structure as an infinite order system or as a 
discrete finite dimensional model, it is prudent to define the terminology used 
in this paper. A state function corresponds to a motion variable which is a 
continuous function of both space and time. Discrete states or degrees-of- 
freedom correspond to point motion variables, which are functions only of 
time, at a finite number or locations throughout the structure (Fig. 1) The 
feedback of spatially discrete structural states involves feedback gains 

w ereas the feedback of a spatially continuous state function involves feedback 

Kernels. 

For control design Linear Quadratic Regulator (LQR) methods exist that can 
be used to formulate optimal structural control solutions for these matrix 
ordinary differential descriptions. Given that model truncation is one of the 
major contributors to the control spillover problem, it is desirable to include as 
many degrees-of-freedom as possible in the control model. This is a costly 
approach, both in terms of money and in implementation since the increased 
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"^^Tto^vd^ptysmte* feed^ i^CTeasc^^mberTfstat* 

measurements 10 s f ofution k ba™”n a^nfiXU' ordefmTdXthat 

lmplementoti^ of a Qlt ^1^ feedback kemel(s) with a spatially continuous 
state°fmcti 0 n(s?wovdd completely avoid the model truncation, spatial aliasing 
and cost of implementation problems. 

This approach contradictstw^co^ 

^These 

meMu^mtn“rt^kally correspond Sacement 

state variables to measure. 

It is also important to realize that the Wh4h»*»' ^"‘el m 
desired to accommodate measurements other than .the state splacement 

the model. Such a transformation can allow the ^e of not only P 
or rotation but also curvature « m^ from 

controUer. This paper discusses the estimation^eiia ct^^ rf kemeta to 

finite dimensional control ^liitions Posing the feedback in terms of 

accommodate the measurement of c \ r sensors known as area 

curvature allows the^J. < ^df the spatially continuous 
measurement of the curvature required by the infimte order control er. 

Implementation issues associated with these sensors are ^disci^in 
ord P er to demonstrate one technique for realizing th f^ l state feedback 

kernels. In this sense, the continuous kernel ' amiability of at least 
^^^^Sr^«s^n V more tLn iust a 
mathematical exercise. 


Background 


form 


Mx(t) + Cx(t) + Kx(t) = f(t) 


( 1 ) 


where M, C and If are the mass, dar nping 
The vectors * and f contain discrete 

respectively. This system can be placed in first order, sta p 

V n / a \ { O « ^ 


ifn - At( t\ 4- Rtl(t) 


where 
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(2b) 


x 

Ul 
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The Linear Quadratic Regulator minimizes a cost 


J = ±J{z T Qz + u T Ru}dt 

0 

for this system by formulating a feedback gain matrix G such that 

u(t) = - R~ I B T Pz(t ) = -Gz(t) 

where P is the solution to the steady-state matrix Riccati equation 

PA + A t P + Q - PBR^BpP = 0 


(3) 

(4) 

(5) 


The feedback form in Eq. 4 consists of multiplying the feedback gains contained 
in G by the state vector in z, whose entries correspond to the temporal motions 
of spatially discrete points throughout the structure. The resultant products 
are summed to arrive at the appropriate control actions which are placed in 
the vector u. This feedback architecture is simply an artifact of the need to use 
a finite dimensional (reduced order) structural model to implement the control 
design. 


In actuality, structures do not undergo motion only at discrete points 
corresponding to the model's nodes but also deform continuously throughout 
the region between nodes (See Fig. 1). The exact motion of the structure is 
described by state functions which are continuously distributed along the 
length of the structure. Therefore, for infinite order structural systems, the 
mathematically exact control inputs are not the sum of products of discrete 
gains with discrete motions but the general form of the control is the spatial 
convolution of the state function with a feedback kernel. 


In order to demonstrate the concept of using infinite order structural models 
for control, a simple structural beam can be used as an example. The partial 
differential equation description of a uniform beam is 


r , T d 4 v(x,t) ^ A d 2 v(x,t) r , 

EI — T~4 — T~2 — = f(x,t) 

dx* dt 1 ( 6 ) 

This description can be placed into state-space, spatial operator form 2 


^-z{x y t) = aix)z{x,t) + b(x)u(x,t) 
at 


(7) 


by choosing the state functions as those which determine the potential (strain) 
and kinetic energy in the beam (curvature and transverse velocity) 
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The parameter E ia the modulus of elasticity, / is the area moment of inertia, 
A is the cross-sectional area and p is the volumetric mass density. 

The cost is defined by 


oo oo 

J = - { J(< qz,z > + < ru,u >)dxdt 


2 O-o 


(9) 


where the matrices q and r are matrix operators penalizing the state and 
control functions, respectively Note that , and r are not constant .tat. » 
sDatial onerators and that the inner integral indicates that the beam is 
assumedCbe of infinite extent. An infinite extent beam was chose n to 
facilitate the acquisition of a closed-form, exact solution. However, the 
operator form for a finite extent beam can also be posed, although the solution 

P ^ffi re s=tettu h nT?mm the solution ps to the fimctional 
Riccati equation 

0 — pox + o pz + qz — pbr ^b pz V* (10) 

where the symbol signifies the adjoint operator. The feedback is the spatial 
convolution of a kernel matrix k with the state function z 


oo 

u(x,t) = -r~ 1 b*( pz/ x,t) = - \>c{x-w)z(w,t)dw 


(ID 


where x corresponds to the location on the structure general 

and w indicates where states are measured. Equation 11 is the genera 
solution because it represents the control action at any location %^ b ack is 
the state functions along the entire extent of the structure. This teedbacK is 
analogous “that in Eq. 4 in the sense that it represents the contmuous sum of 
gains times the states of the structure. 

The implementation of these continuously distnbuted feedback kernels 
requires^ the use of a continuously distributed sensor. Several researchers 
have demonstrated the use of continuously distributed cun^atu re sensors and 
actuators C K. Lee 3 - 4 , S.E. Miller 5 , S. Collins 6 and D. Miller have worked on 
the deveiopment of area averaging sensors. These authors use spatially 
distributed sensors to achieve certain measurement charactensucs^ CX 
Lee 3 .* and S. Collins® used sensors shaped as particular mode shapes to obtain 
a measurement of the generalized coordinate of that mode. S_ Collins® and a 
Miller 7 developed sensors which roll off with frequency without g 

rdiatelag 6 It will be shown in this paper that area averaging sensors can be 
used to implement the feedback solution to a partial differentia equa 
description of a finite extent structure 
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Approach 

An over optimistic goal for this research would have been to attempt to solve 
the infinite dimensional structural control problem. This goal is not realistic 
because first it would require the exact partial differentia] equations (PDE) and 
boundary conditions (BC's) that describe the dynamics of the structure and 
second it would be impossible in most cases to find the LQR solution even if an 
exact model existed, de Luis 2 , for example used an infinite extent beam in 
solving the infinite dimensional control problem in order to find a closed-form, 
exact control solution. The same infinite dimensional control problem can 
also be posed for a finite extent beam 8 . However, this problem is much more 
difficult to solve due to the existence of boundary conditions. 

A more realistic approach is to model the structures with the more fami liar 
finite reduced order models (Eqs. 2 through 5) and to hope that by increasing 
the fidelity (number of degrees-of-freedom) of the model, the continuous 
feedback kernel can be inferred from the distribution of the discrete gains. 

The following section discusses the derivation and implementation of 
continuously distributed feedback kernels. Several important steps are 
involved in this derivation. First, spatially discrete displacement and rotation 
gains derived using standard matrix Riccati techniques on finite element 
structural models must be converted into spatially continuous feedback 
kernels. Second, these kernels must be transformed into feedback of 
distributed curvature to facilitate implementation using area averaging 
sensors. An alternative approach, also discussed in the next section, is to first 
convert discrete displacement and rotation gains into discrete curvature gains 
and then to convert these gains into a spatially continuous curvature feedback 
kernel. Numerical examples are interspersed with these formulations to 
demonstrate these techniques. After the section on feedback kernel 
formulation, a discussion of general control issues of interest is presented 
along with an additional numerical examples. 


Reference Example 

Throughout the rest of this paper, these techniques are formulated for the 
cantilevered beam of length L shown in Figure 2. A control moment is applied 
to a point on the beam 1/10^ of the distance from the clamped root to the free 
tip. This moment actuator is used to represent an equivalent piezoelectric 
actuator at the cantilevered end. de Luis et afi demonstrate that one valid way 
of modelling the influence of a piezoelectric curvature actuator is to derive 
equivalent moments at the two ends of the actuator, which are of equal 
magnitude but of opposite sign. In this problem, if it is ass um ed that the 
piezoelectric actuator runs from the root, the companion moment at the 
clamped end of the beam does not enter the problem and is therefore not 
shown. The pertinent parameters of the problem are listed in Table 1. The 
performance metric is the transverse displacement of the tip of the beam (vtip). 
The entry in the state penalty matrix (Q) corresponding to this displacement is 
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assumed to be unity. This state penalty in equation form, from Eq. 3, 

fM O 


* T Q z = v li 



Feedback kernel formulation 

This formulation process is shown in Figure 3. Theupperleft 
and curvature gains and kernels. 

agag£ st& 
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action. This path is discussed in detail in the rest of this section. 

Table 1. Parameters for cantilevered b eam example. 

Bending stiffness ' \ 1.0 Nm 


Mass per unit length 


Length 


Actuator location 


Control effort penalty 
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and rotations gains 
from classical FEM 
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Use finite element 
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interpolation 
functions 


Discrete state 
transformation 


(Continuous 


Continuous displacement 
kernel function 


lb I Integration 


Discrete 
curvature gains 


Continuous curvature 
kernel function 


Use finite element 

curvature \ 

interpolation 

functions 

Figures. Matrix of feedback options. 


Following path two, the second step involves transforming the discrete 
displacement and rotation gains into discrete curvature gains (2a). The third 
step then involves using these gains to find the continuous curvature kernel 
(2b). This path is used, in the following discussion, as a check of the first path 
since both paths should yield approximately the same curvature kernel. 


Evaluation of the Discrete Displacement and Rotation Gains 

The first step in evaluating the discrete displacement and rotation gains is to 
develop a finite element model of the cantilevered beam. The mass and 
stiffness finite elements that are used in the following analysis are 

'156 221 54 -23/1 [22 61 -12 61' 

_ pAl 22/ *l 2 131 -3l 2 _ EI 61 4l 2 -61 2l 2 

mele 420 54 131 156 -221 ~ ~F -12 -61 12 -61 

-131 -3l 2 -221 4l 2 J [$/ 2l 2 -61 4l 2 J (13) 
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with the corresponding finite element nodal degrees-of-freedom 

", "i+l Vl] (14) 

where l is the element length and i. .equal* 

utsls^dXX mod“l"s undamped. The entry in the state penalty 
matrix Q corresponding to this displacement is set equal to one. 

Using a ten element model of tfie beam the gains ob£n e d from ^arfth" 

window shows the rotation gams. N « t ' c V^ s n “ n ^e upper and lower 
"X^ ^s^ernemrnte 6 and rotation rate gains, respectively. 
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These displacement and displacement rate gain distributions indicate the 
shape that the respective continuous kernels will have, but not the 
magnitudes. This is only an approximate indication of the continuous shapes 
since the displacement kernel combines the information from both the discrete 
displacement and rotation gains. In other words, a single kernel contains all 
of the gain information displayed in a single column of Figure 4. 


Evaluation of the Spatially Continuous Feedback Kernel 


The next objective in the analysis is to find the spatially continuous feedback 
kernel from the spatially discrete gains evaluated in the previous section (Path 
a in Fig 3). To this end, the beam finite element displacement and rotation 
gains will be used to derive the continuous displacement feedback kernel 
which convolves with the displacement state function. Since the reference 
example has a point actuator, the feedback convolution in Eq. 11 degenerates to 
the integral of a kernel times the state function. It is also convenient to use a 
kernel that is defined over the length of the beam, rather than having the 
kernel be defined, as in Eq. 11, in terms of the actuator location (x a ). Using 
this kernel transformation, the feedback is given by 


u{t) = - J K{w)z{w,t)dw = - 


r 

\d 2 v 1 

[ K C k Dr](w> 

dw 2 

dv 

J 

. dt . 


\(w y t)dw 


^2 - . 

kq{w)— 2 ^ w ^)dw- [ kdr(w)— (w,t)dw 
dw J dt 


L ^ 

^ r dt) 

= -iK D (w)v(w,t)dw- K DR (w)—(w,t)dw 
o J at 


-U[)(t) + U[) R (t) 


(15) 


Note that the state functions shown in Eq. 8 include the curvature of the beam. 
Eq. 15 shows part of the feedback to be the integral of curvature times a 
curvature kernel. Alternatively, this can equivalently be expressed as the 
intep-al of the displacement state functions times a displacement kernel. This 
displacement kernel is derived in the next paragraph. 


The integration over the entire length of the beam can be divided into the sum 
of integrals over segments of the structure corresponding to the finite element 
domains as shown in Fig. 5. 
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The first half of Eq. 15 then becomes 

L 

u D (t) = -\K D (w)v(w,t)dw 

0 

Wi+1 

=...- \ K[){w)v(w t t)dw 

w t 


u, i+2 

1 k d( w ) v ( W,t)dw~... 

Ui+l 



=...UD t + u D l + i + - (16) 

The element interpolation function description of the displacement anywhere 
within the element located between Wi and Wi+j 


v(£,t) 
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for 0 < % < l 

( 17 ) 
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can be substituted into each of these element convolutions. Then, the control 
action associated with an element is a function of that element's nodal 
degrees-of-freedom 
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dw where t; =w-Wj 


(18) 


If the form for kq were known, then the integral in Eq. 18 could be evaluated to 
find the the gains for the nodal degrees-of-freedom. Conversely, in this case 
these gains are known from the solution to the matrix Riccati equation and 

instead it is the form of the kernel kd that is being sought. To estimate this 
kernel, a form for the kernel, containing unknown parameters, can be selected 
so that the spatial integral in Eq. 18 can be evaluated. Then, these parameters 
can be found by equating the elements of this integration to the discrete gains. 
A cubic form for the kernel is chosen 

KQ.(w) = Oi(w-Wi) 3 J rb i {w-w i ) 2 +Ci(w-Wi) + di forWi<,W<,Wi + i 

Given the polynomial order (cubic) assumed for the four degree-of- freedom 
finite element, a cubic polynomial for the internal curvature distribution is the 
highest order polynomial for which the unknown coefficients can be uniquely 
found. 


If the form in Eq. 19 is inserted into Eq. 18, and the integral is evaluated, the 
result will be the contribution that the continuous kernel across that element 
makes to the gains associated with that element's nodal degrees-of-freedom. 
In other words, at one of the element's nodes, Eq. 18 yields partial gains for the 
nodal degrees-of-freedom which, if summed with the gain contributions from 
the adjacent element, will give the total gains associated with that node's 
degrees-of-freedom. Thus, the gain contributions from the elements 
neighboring a shared node can be used to find the total displacement and 
rotation gains associated with that shared node 



(20a) 
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(20b) 


where g v . and g • are that node's displacement and rotation gains, 
respectively, as shown in Fig. 4. 

These two relations give two conditions for finding for the elemental kernel 
coefficients ai, bi, ci and di. Two more conditions are required in order to 
ensure a unique solution. These two additional conditions are found y 
requiring continuity of the kernel magnitude and slope at a shared no e. 
These are found by using Eq. 19 to evaluate the magnitude and slope at the 
right end of the i th kernel and equating that to the magnitude and slope of the 
( i+1 )th kernel at its left end yielding 

a t l 3 +2bil 2 +c l l + d l -d i + I =0 (21) 


3a t l 2 + 26j/ + Cj -c; + ; =0 (22) 


These four conditions can be expressed in matrix form as 
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where the first two rows are found by evaluating the integrals in Eqs. 20a and 
20b A global matrix can be assembled, using Eq. 23 as the sub-matrices, o 
yield a linear equation relating the coefficients of the cubic-fitted kernel 
functions to the discrete gains 

Tc = g (24) 
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( 25 ) 


The desired coefficients are then given by 

c = T J g 


The approximate shape of the spatially continuous displacement feedback 
kernel can be calculated by evaluating this piecewise cubic kernel along the 
length of the structure. This evaluation is made by using the coefficients in c 
which are appropriate for the given segment of the structure within which the 
kernel is being evaluated. 

Using the discrete gains of the ten element finite element model (shown in Fig. 
4) to evaluate the coefficients in Eq. 19 of the piecewise cubic displacement and 
displacement rate kernels, the functions in Figure 6 are found. These 
functions are the piecewise cubic kernels combined into a single curve. 

Notice the erratic shape of the displacement kernel. This erratic shape may 
correspond to some weighted sum of mode shapes. Given that the tip 
displacement (performance metric) can be represented as a sum of 
displacement mode shapes, and that the applied moment (control input) can be 
represented by the sum of curvature mode shapes, the shapes in Figure 6 
could correspond to some combination of the displacement and/or curvature 
mode shapes. In other words, these shapes may correspond to some type of 
mode shape 'feedthrough' from the control input to the performance metric. 



Figure 6. S patially con tinuo us feedback kernels as a function of location along the beam for 
controlling tip displacement. Hie individual windows show the kernels for (a) displacement | 

and (b) d isplacement rate. 1 
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Equivalent Feedback Using Alternative State Functions 

rpi f fk architecture using the kernel derived in the previous section has 

Thl fn^ fn E^16 Thi8 involves the spatial integral of the kernel times the 
the form in Eq* lo- However the continuously distributed 

displacement state function. H.o e , . ■ . Therefore it 

displacement state function may not be a measurable q^nt.ty. TOerefo^ .t 
may be convenient to express the same control in terms of another, more 

measurable state function. 

Extensive work in the area of area averaging sensors 4 - 5 -®’ 7 ’ 8 has shown that 
continuously distributed measurements of curvature-induced strain can be 
made 1 using polyvi ny lide ne flounde (PVDF). Therefore the displacement 
feedback kernel of Eq. 16 must be transformed into equivalent feedback of the 

curvature state function. 

Integration by parts can be used to transform the feedback form m Eq. 16 into 
equivalent feedback of the curvature state function 

rotation and displacement, in order to retain rigid body control, this 
transformation is given by 

L 

u D (t)= \K D (w)v(w,t)dw 


= v(0,t) \ k d ( w)dw + 1 \ kd( y)d)dw 

o Ow 

LwL d 2 v ( w ) 


\j\K D (t)dt dy°-ff-dw 


ooy 


(26) 


While the point measurements of displacement and rotation ^ st ,^ e 

order to retain rigid body measurement, the actual location on the ^ructure 

where these measurements are made >s arb. rary. The « 

"Sit" :^y^»™oiX^,ng n th b e e il™n£ ‘strain ipropnateiy. 
For example, 


du(L,t) dufOdl 


dx 


dx 


d^v(_ wd) 

dw 2 


dw 


(27) 
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Substituting this translation of the rotation measurement into Eq. 26 gives the 
equivalent feedback as 


?h)( 0 t) 

uj)( t) = v( 0,t) S k d ( w)dw + — — 1 — / / k d ( y)dydw + 

O ow 

LLL d^l)( w t ) 

\]\K D (x)dxdy - v -^f+ dw 

owy dw (2g) 


Notice that in this equation, the first two terms, representing the feedback 
gains associated with point displacement and rotation measurements, can be 
evaluated directly from the displacement feedback kernel. The outer integral 
in the third term corresponds to the integration of the distributed curvature 
kernel with the curvature state function. The inner two integrals evaluate the 
curvature kernel from the displacement kernel. This curvature kernel is 
given by 


LL 

k c (w)= j j< D (x)dxdy 

wy 


(29) 


The boundary conditions in the reference example were conveniently chosen to 
exclude rigid body motion thereby eliminating the need for any point 
displacement or rotation measurements. The motion of the structure is 
completely describable by the curvature state function because 


v( 0,t) 


dv( 0,t) 
dx 


(30) 


Substituting Eq. 30 into Eq. 28 yields 

u d( O-iii tjdxdy - ^^- dw 

owy dw 

as the feedback law in terms of the displacement kernel. To calculate the 
shape of the continuous curvature kernel, Eq. 29 is employed. Equation 31 can 
also be used to evaluate the curvature rate kernel if the displacement rate 

kernel (k dr ) is used in place of the displacement kernel (»r D ). 


Figure 7 depicts the resulting curvature and curvature rate kernels for the ten 
node finite element model (Figure 4 and 6). Notice that, unlike the 
displacement kernel, the curvature kernel is smoother. This is predominantly 
due to the smoothing process inherent in the double integration of Eq. 31. Also 
notice that the magnitude of the kernels in Fig. 7 are largest where the 
cantilevered beam tends to exhibit the largest curvature: the root. In Fig. 6, 
the magnitude of the displacement kernel is not the largest where the beam 
tends to exhibit maximum displacement; namely the tip. 
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Figure 7. Spatially continuous feedback kernels as a function of location along the beam for 
controlling tip displacement The individual windows show the kernels for (a) curvature and 
I (b) curvature rate, 


Although not shown, for this reference example, increasing control authority 
by decreasing the control effort penalty (R) does not change the kernel shapes 
Instead, it changes the absolute and relative magnitudes of the kernels A 
change in the shape of the kernel will be achieved by a change in the spatial 
nature of the problem such as moving the actuator or selecting a different 
performance metric. This observation is supported by an additional example 
presented later in the paper. Actual implementation of these sensors is the 
topic of a follow-on paper. 

The results in Fig. 7 correspond to the objective represented by the lower right 
box in Figure 3. The next step would be to implement these two kernel shapes 
using area averaging sensors. The details of this process wil e scusse in 
the section on implementation issues. Prior to that, the next section discusses 
the alternate path in Fig. 3; namely path number 2. 
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Equivalent Feedback Gains using Curvature States 


An alternative procedure to evaluating the continuous curvature kernel is to 
first derive the discrete curvature gains from the discrete displacement and 
rotation gains, as shown by path 2a in Fig. 3. This can be done in two ways. 
The first involves using the transformation matrix given by de Luis et a/ 2 



2 

' l 
4 

l 


' Vi ' 
Vi 

Vi+1 
Pi* l . 


(32) 


This elemental sub-matrix can be assembled into a global state transformation 
matrix. The number of degrees-of-freedom are not reduced by this 
transformation because now there exist two independent curvatures at each 
node. Remember, curvature is not constrained to be continuous in the beam 
finite element formulation because applied point moments can induce 
discontinuous curvature. Originally, displacement and rotation were the two 
nodal DOFs. Now, a node has two curvatures, one associated with the left and 
one with the right-hand element. 


The ’o’ symbols in Figure 8 indicate the net curvature gains at each node as 
derived using this transformation. The net curvature gain at a particular 
node is found by summing the individual curvature gains at that node. This 
procedure is justified at nodes where external moments are not applied 
because the two curvature gains correspond to the feedback of curvature 
measurements acquired an infinitesimal distance to either side of the node. 
Without an externally applied moment, it can be assumed that these curvature 
measurements are identical and therefore the net gain is the sum of the two 
gains. 


The second approach to deriving discrete curvature gains is to integrate the 
displacement and rotation gain vectors to get the curvature vector. Unlike Eq. 
29, this integration process involves both the displacement and rotation gains. 
This integral can be approximately evaluated element by element by summing 
the products of the gains with the element width. Other standard numerical 
integration techniques can also be used. 

The symbols in Fig. 8 represent the curvature and curvature rate gains 
found using this integration approach. 
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Figure 8 . Discrete gains at nodal locations along the beam for controlling tip displacement. 
The individual windows show gains for (a) curvature and Ob) curvature rate. Tne 
indica tes gain* found through integration while the 'o' indicates gains found through 

transformation. 


Notice the good agreement between the curvature and curvature rate gain 
found using the transformation and integration techniques. The agreemen 
may seem to improve near the tip of the beam but when calculated it was found 
that the relative error is more less constant along the beam. 

The final step in Figure 3 (2b) involves calculating the curvature and curvature 
rate kernels from the discrete curvature and curvature rate gains. Kather 
than using the technique in Eqs. 15 through 25, it can be observed that each of 
the discrete gains at a node roughly represents the area under the continuo 
kernel for the region between the midpoints of that node s f 

elements Therefore, if the gain is divided by the length an elem f nt ; \ 
result should be approximately equal to the magnitude of the kernel at the 

nodal location. 

Figure 9, when compared with Fig. 7, shows that this is the case. 
Furthermore, Fig. 9 shows the overlay of the gains divided by respecti e 
element lengths for different fidelity models. This reveals that the magnitude 
of the kernel is captured quite well at nodal locations for rather coarse mo 
for this simple reference example. This is an important result since 
practice it would be generally impossible to find the exact feedback kemel fr 
a continuous model. However. Fig. 9 illustrates that as the order of the model 
is increased the kernel shape is asymptotically approaching some shape. It is 
this shape that represents the infinite order solution and that must be 
implemented. 
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Figure 9. Discrete gains divided by element length at nodal loca ti o n s along the beam for 
i. controlling tip displacement Gains for 5, 6, 7, 8, 9 f 10 and 11 finite elements are overlaid. The 
| individual windows show g ains for (a) curvature and (b) curvature rate. 


Implementation Issues 

The possibility of implementing infinite order structural controllers is made 
possible by the existence of area averaging sensors such as those described in 
References 4,5,7 and 8. Once the curvature kernel is obtained, it is a simple 
calculation to alter the kernel for equivalent feedback of curvature-induced 
surface strain. This simply requires knowledge of the distance of the surface 
mounted sensor from the centroidal axis in the structure. Once this kernel is 
found, the sensor can be built. 

Polyvinylidene flouride (PVDF) 9 is suggested for this sensor for several 
reasons. First, PVDF is a strain sensitive material which can be continuously 
distributed along the surface of a structure and whose accumulated charge on 
a surface electrode equals the integral over the length of the PVDF of the 
electrode width times the surface strain in the structure. Second, PVDF has 
an elasticity which is relatively small when compared to the elasticity of 
conventional structural materials. This allows the sensor to be rather non- 
intrusive into the dynamics of the structure. Third, the shape of the electrode 
can be easily altered to equal that of the kernel while leaving the actual PVDF 
material uniformly distributed. This achieves the strain sensitivity 
appropriate for implementing the kernel while keeping the small dynamic 
influence that the PVDF does exert on the structure uniformly distributed. In 
addition, removal of electrode from near the edge of the PVDF greatly reduces 
the possibility of the sensor shorting its bottom and top surface electrodes. A 
fourth and final reason for using PVDF is its high strain sensitivity which 
provides an excellent signal to noise ratio for control purposes. 

One drawback of implementing the feedback kernel through the shaping of the 
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electrode is that once the electrode is shaped and the material is mounted on 
the structure, the kernel is effectively fixed and cannot be altered. Feedback 
ea ins which reside in a computer can be readily altered if alteration is 
reauired However, C. K. Lee in Reference 4 has developed a method which 
could be used to circumvent this inflexibility in the gains. He uses an area 

averaging sensor whose electrode is segmented into n ^ er ^ !2’^[® 8 ^ 
the voltales on these squares are summed as appropriate for a particular gain 
b£ C distribution needs to be altered, the manner m which 

these voltages are summed can be changed. 

Throughout the discussion of full state feedback for infinite order systems, 
there ^was an implicit assumption that high frequency dynamics in the 
structure consisted solely of additional modes whic wou ^ 
modelled given the use of a sufficient number of finite elements. Howevertms 
^seldom if ever, the case in actual structures. Often, torsion or out-of-plane 
bening modls exist irrespective of whether only in-plane bending was 
modelled. These dynamics may feed through to the output of the sensor. 
Therefore the spatial wavenumber filtering concepts presented m Reference 8 
could bemused to roll off, without phase lag, the frequency response of the 
spatially continuous sensor. 

Ficnn-P 10 illustrates the way in which a PVDF area averaging sensor was 
Fomented in Reference 8. The electrode is shaped as a decreasing 
exponential in two directions. Note that the sensor may h^to be 
if thp PVDF sheet is not as long as the kernel. Given that PVDr is V° . 
maUnTa negatiie part of the kernel can be implemented by e.ther (hpp.ng 
that segment of the PVDF or reversing leads (see Reference 8). 

For the reference example discussed thus far, two PVDF electrodes could be 

shaped: on“as sho'wn in Figs. 7a and 7b. B -f^^ t ^kemel and 
oithpr side of the cantilevered beam, one sensor for the curvatu 
oneltor thecurvature raU kernel, the two sensor voltages can be summed 
appropriately and used to drive the control moment. 

The unique feature of this technique is that the processes 0 f” ult 'P 1 ^"8 th * 

pains times the curvature measurements and accumulating these products is 
fer“rmed by t^e"ensor. Th.s feature significantly reduces the control 
implementation effort associated with numerous point sensors. 


Segment \ 


Segment 2 


Segment 3 


Copper Lape 


Positive 

lead 


Etched electrode 


Copper upe _ 
(below PVDF) 


Negative 

lead 


Figure 10. Illustration of the implementation of an area averaging sensor using 1 

in an exponential fashion. 
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Issues qggnrintpd with controllers based on classical beam finite 
elements 


The previous section has shown how PVDF sensors can be used to implement 
infinite order controllers. It was also shown in Figure 9, that finite elements 
models can be used to predict the shape of the infinite order feedback kernel. 
The hope is that by progressively increasing the order (accuracy or fidelity) of 
the finite element model, the shape of the feedback kernel will approach some 
asymptotic shape. It is this shape that represents the infinite order feedback 
kernel and that must be implemented with PVDF. 


Classical finite elements are the obvious elements to be used in such a model 
refinement process. This study has identified two implementation problems 
that are uniquely associated with these classical beam finite elements. The 
first is that the stiffness matrix obtained with these classical beam elements 
becomes ill-conditioned as the element size decreases. Decreasing the element 
size is typically associated with increasing model fidelity. This is illustrated by 
looking at the conditioning number of the stiffness matrix of a cantilevered 
beam obtained by using the following classical beam finite element: 
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(33) 


The conditioning number for a matrix is the ratio obtained by dividing the 
largest eigenvalue by the smallest eigenvalue of the matrix. The higher the 
conditioning number of a matrix, the more ill-conditioned the matrix is and 
the more likely that matrix will be susceptible to computer round-off errors. It 
can be shown 10 that the conditioning number is proportional to: 


Cond « 7) 

l 2 


(34) 


Thus as more elements are used and the element length (/) decreases, the 
matrix becomes ill-conditioned and results from the LQR routine will become 
less reliable. 


A second problem associated with classical finite element models is a problem 
of non-uniqueness. From finite difference theory it is known that rotation can 
be estimated from discretized displacements as: 




(35) 


In Eq. 35, v t is the nodal rotation and Uj’s is the nodal deflection. The 
truncation error, which is of order ft, will decrease as the element size ( l ) 


22 



decreases, indicating that the finite element nodal rotations iv t ) can be 
expressed as linear combinations of the nodal displacements (i/j) with 
incrcssin^ accuracy. The manifestation of this problem lies in the 
interpretation of the feedback gains calculated by the LQR algorithm. Given 
that the nodal rotations may become linearly dependent on the nodal 
displacements (or vice versa), the gains obtained by the LQR algorithm may 
yield an optimal solution but the displacement and rotation gains may not be 

Thase^two problems are investigated by comparing the results of two 
discretized models used to solve the reference cantilevered beam example. The 
first model is the classical finite element beam model, while the second is a 
second order accurate finite difference model. In the finite difference mode, 
the stiffness term in the governing differential equation (Eq. 6) is approximated 


by: 


»m 

EIvi 



- 6^-1 


+ 4 Vj — 6t?i_i -t-^i-2 
l A 



(36) 


The effects of ill-conditioning and non-uniqueness are investigated by 
comparing the results of models in which the fidelity of the model is increased 
by increasing the number of nodes. Both these models should exhibit the ill- 
conditioning problem since the finite difference model also has a conditioning 
number that will increase (deteriorate) as the element size decreases since the 
conditioning number is approximately 1 HI 2 ) 10 . The finite difference model, 
however, should not exhibit the non-uniqueness problem associated with the 
finite element model. These conclusions are supported by the results ot the 
investigation. Although not shown, both the models exhibit ill-conditioning 
problems and the Riccati solver failed to yield a solution for a model with 40 
nodes (or 80 degrees-of- freedom) for the finite element model and 80 degrees-ot- 
freedom for the finite difference model. However, the finite element model 
may exhibit the non-uniqueness problem as the fidelity of the model is 
increased. In Fig. 11 the distribution of curvature gains becomes erratic as the 
number of nodes are increased above 10. The finite difference model on the 
other hand, as shown in Fig. 12 does not exhibit this behavior Even with these 
erratic gains, the closed loop finite element models are stable with identical 
closed loop poles for the first five modes. This observation leads to the 
conclusion that this behavior may be due to the non-uniqueness problem 
associated with these elements. 

Note that the slow convergence to the “infinite” shape of the finite difference 
model is due to the method in which the point moment is applied to the 
structure. An applied point moment is achieved by applying appropriate 
forces to nodes neighboring the node to which the moment must be applied. 
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location on beam 


location on beam 


(a) (b) 

Figure 12 Convergence of finite difference model with increasing model fidelity. 
(’* '■ 10, V » 20, ’o' » 40 and 'x' ■ 60 Degrees -of-freedom) 
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an additional numerical example 

An additional numerical example involves the control ° f 
transverse displacement between the tip and the middle of the beam. This 
state penalty has the form 


Q{ v tip - ^middle) = <l[ v tip ~ 2 v tip v middle + v middle ) 


(37) 


For this example, the scalar q is unity. 

Figure 13 shows the discrete gains. Again, the displacement and rotation 
gains are rather erratic. However, the curvature and curvature rate grnns a 
smooth Figure 14 shows the continuous feedback kernels. While the 
curvature rate kernel has a shape similar to that in the previous example, the 
curvature kernel now undergoes a change in sign. All the curves seem to 
have an inflection point near the midpoint of the beam U -O.b). 


CONCLUSIONS 

A technique has been presented for inferring the exact, sp^i^ly cctnliiriuous 
LQR feedback solution to the control of structures from the discrete feedbac 
Sns derived using finite dimensional structural descriptions. These 
feedback kernels possess several unique attributes. First, it has been s ow 

that feedback of the state functions can be transformed t0 'i'Xl'eedback^be'^ 
other stole functions. This aids in implementation because the feedback can be 
derived in terms of the state function that is most easily measured. Area 
averaging sensors provide one means for implementing these spatially 
continuous feedback kernels. Second, these continuous 

SDatial aliasing. Spatial aliasing is one of the primary causes of spillover in 
structural control. Third, all of the feedback computation can be effectively 
performed by an area averaging sensor. 

The research presented in this paper must be seen as the first step in an 
attemDt to formulate and implement full state feedback for infinite orde 
structural systems. Several issues must be resolved before this approach can 
be considered a viable alternative to reduced order controllers ^ForexampK 
the accuracy with which the area averaging sensors must match the desired 
kernel must be investigated. Robustness of this control app ™ a ^ “ ust T Jj 
determined and the theory must be demonstrated in the laboratory. The 

researchers are presently working on these topics and plan ta ^f^r and 
infinite order controller on a cantilevered beam using the actuator a 
performance metric presented in the reference example. 
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curvature gain rotation gain 



location on beam location on beam 

(c) (d) 



(e) (f) 


Figure 13. Discrete gains at nodal positions along the structure for controlling relative 
displacement between the tip and midpoint. The individual windows show gains for (a) 
displacement, (b) displacement rate, (cl rotation, (d) rotation rate, (e) curvature and (f) 

curvature rate. 
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on a structure in order to modify its dynamc behavior to 
moot it. parformanc. requirements. Unfortunately 
active control introduces the possibility of «nttng the 
structure in an unstable manner making it mtwd tha 
either confidence in the prediction of on-orbit behainor 
be improved or the types of tests required for 
qualification be identified. 

Before proceeding with a discussion of the 
experimental approach to developing 
procedures, it is necessary to present the rationale that 
lead to MACE. After all, conducting experiments on- 
orbit, even thoee which are performed on the b O 
middeck, is technically risky, expensive, rec l u >"* 
extensive planning, and produces less data than w 
be obtained in a comparable ground experiment, lhe 
program must clearly exploit the unique aspects of the 
on-orbit environment in order to justify its conduct. 

The objective of this paper is to portray the rationale 
for conducting this type of flight experiment and to pose 
the scientific questions to be addressed through this 
research. Additionally, the test article will be esen d, 
along with the ground and on-orbit experiment support 
equipment. This paper concludes with a discussion of 
planned on-orbit activities. 

Objectives and Rationale 
" The goal of MACE is to develop a well verified set of 
CST tools that will allow designers to either be able to 
predict on-orbit behavior or allow sufficient versatility in 
the design to allow identification and tuning of the 
structure on-orbit. A number of different options exist 
for deriving this set of tools. The first and least expensive 
is to rely on analysis for the design and qualification of 
spacecraft which incorporate CST. Unfortunately, this 
approach is far less than satisfactory TTie scientific 
literature is riddled with examples of both closed and 
open-loop experiments whose performance varied 
greatly from that predicted by state-of-the-art analytical 
methods. The reasons behind this are varied. Often the 
structural or sensor/actuator characteristic which 
contributes to this performance degradation is not the 
next detail that would have been included m the 
analytical model. Its existence is usually not predicted 
but instead is discovered through experimentation. This 
experience illustrates that analysis alone is not sufficient. 

The question that next arises is what sort of testing 
needs to be performed, along with analysis, in order to 
develop an effective and efficient spacecraft 
qualification procedure. Four different options exist. 
Li sted in ascending order from lowest to highest cost and 
complexity, they are: ground-based open-loop 

experiments, ground-based closed-loop experiments on- 
orbit open-loop experiments and on-orbit closed-loop 
experiment*. 

Ground-based open-loop testing is the simplest type 
of experimental program that can be earned out to 
verify the validity of analytical models. It is an 
absolutely necessary step, since the quantities that are 
most required for closed-loop control design are exactly 
those which are hard to predict analytically For 
example, structural modal frequencies can be predicted 
using numerical methods with a relatively hign degree 
of accuracy. Conversely, modal damping values are 
extremely hard to predict analytically on large comp. ex 
structures where many energy dissipation mecranisms 
are present. Unfortunately, closed-loop contro. ers .or 


structures usually require accurate knowledge of the 
modal damping because damping determines stability 
margins and therefore performance. This problem is 
exacerbated in structures that are lightly damped, such 
asLSS. 

It is easily concluded, therefore, that ground-based 
open-loop testing is essential to quantify the accuracy of 
analytical models. However, these tests by themselves 
are not sufficient to validate the appropriateness of an 
analytical model or the performance of a closed-loop 
system. Skelton 3 has demonstrated that no measures of 
accuracy of the open-loop model are sufficient to 
guarantee stability of a closed-loop system at arbitrarily 
high gain. This implies that the acquisition of the open- 
loop model can never be sufficient to predict closed-loop 
performance. Therefore, ground-based closed-loop 
testing is absolutely necessary for the successful 

application of CST to reidistic structures. 

Since CST structures will be used in the space 
environment, it is important to investigate whether thoee 
characteristics that are present on-orbit and cannot be 
adequately simulated on earth affect the open and 
closed-loop teste. In Table 1 various vehicle parameters 
are listed along with four significant differences that 
occur between on orbit and ground-bawd tests. The 
table indicates that these differences do affect the vehicle 
parameters. 

T*btal TWv«ite«**wWnLU*~"‘<“ dd ^" lt ?“”■•*•" *** 
j-** dlfftr b«rw t « oiHOttolt and yroondt fa- _ 

Aero/ Suspension Gravity Thermal/ 
Acoustic Radiation 


Stiffness 

Damping 

Mass 

Forcing 

Kinematics 


no 

yes 

yes 

yes 

no 


yes 

yes 

yea 

yes 

yes 


yes 

yes 

no 

no 

yes 


yes 

yes 

no 

no 

no 


The important issue is whether the differences in 
Table 1 cause regular or singular perturbations to 
problem. A regular perturbation is one whow affect on 
the vehicle parameter disappears as the P««urbaaon is 
allowed to approach zero. This is m cont ”»* ‘ 
singular perturbation whose presence substantially 
modifies the vehicle parameter even as the P« rt ^ rb *° 
approaches zero. If the perturbations are 
they can be modelled and the results from the ground- 
based tests can be more easily used to predict on-orbit 
behavior. However, they may still have a very 
substantial, although predictable effect on the stomctural 
parameter. For example, small changes in the 
often lead to large changes in the modal °I * 

the mod* shapes, two quantities that have a ^ect effect 
on closed-loop stability and actuator and 
performance. Therefore, if the plant is highly 
regular perturbations due to influences listed in Tabla 1, 
it is probably necessary to conduct open-loop on *? rb ‘ t 
testing. If the perturbations are singular, it is essentisj to 
conduct open-loop testing on-orbit in order to identify 
and adjust for these perturbations. ^ 

The only issue that now remains to be addressedis 
whether on-orbit closed-loop testing is still required. The 
answer to this question depends on whether any singular 
perturbations are identified during the on-orbit open- 
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loop •xperimenU, or whether any regular perturbations 
cause significant unpredictable changes in the plant. If 
the answer to either of these questions is yes", then on- 
orbit closed-loop testing is essential. 

A preliminary analysis does not reveal any singular 
perturbations arising from one of the four sources shown 
in Table 1. Non-convective potential aeroacoustic 
equations do not give rise to singularities, nor do 
conservative fields such as gravity. So long as 
suspension devices are passive or collocated active, they 
do not introduce singularities. Since the 
thermal/radiation terms only affect otherwise symmetric 
stiffness and damping parameters, they also do not give 
rise to singular perturbations. 

However, a situation in which a regular 
perturbation can have significant effect on the closed- 
loop performance of the structure can be easily 
imagined. The stiffness added by a suspension system, 
even if small, can change the modal structure. 
Additionally, for an articulated test article, a suspension 
system could introduce an unexpected kinematic 
constraint. Gravity can change preload on a joint, and 
hence damping. Gravity will also cause otherwise 
straight members to curve, causing significant changes 
in the modal structure, such as nonplanar coupling of 
modes. Therefore, while no singular perturbations have 
been identified, there are a number of regular 
perturbations which can cause significant changes in the 
plant that could result in control performance being 
degraded. 

Therefore, the conclusion that is reached is that 
ground-based open and closed-loop testing is not 
sufficient for the verification of CST technology. At a 
minimum, on-orbit open-loop testing would need to be 
conducted to test for the presence of any singular 
perturbations, or any significant regular perturbations. 

If these perturbations are found to exist, then on-orbit 
closed-loop testing becomes essential as argued by 
Skelton. If they are not present, then the closed-loop tests 
might still be needed if a suitable ground-based 
performance metric or disturbance environment is 
unobtainable, or, more likely, if the additional cost of 
conducting the closed-loop experiments were 
incremental. 

Having demonstrated the likely necessity of on-orbit 
closed-loop testing, a test article on which to perform the 
experiments must now be selected. A survey of proposed 
future spacecraft was undertaken and an evaluation was 
made on which type of spacecraft exhibit the most 
requirements for CST and which were most limited by 
earth-bound testing.* 5 Some of the spacecraft types that 
were considered included two point alignment occulting 
instruments, multipoint alignment interferometric 
devices, shape control of reflective surfaces, flexible 
manipulators, and multipayload platforms. This latter 
type was selected because the large angle motions of the 
payloads stress state-of-the-art suspension devices and 
because of its applicability to missions of near term 
interest. 

Proposed missions which will use this type of 
spacecraft include low and geosynchronous platforms in 
the Mission to Planet Earth, the evolutionary 
International Space Station, and the planetary orbicing 
platforms of the Exploration Initiative. As these 
platforms become larger and more complex, the 
propensity for individual on-board controllers to interact 


with each other and with the bus attitude control system 
will grow. This propensity is exacerbated by increasing 
payload mass fraction associated with larger 
instruments and robotic devices, decreasing structural 
bus stiffness associated with larger platforms, increasing 
authority of the controllers associated with tighter 
pointing and positioning requirements, and the 
increasing need to reject disturbances which originate at 
other payloads. This rationale makes clear the need to 
develop a well verified set of CST tools. This 
development must include: 

development of a comprehensive analytical 
CST framework for the design and analysis of controlled 
multibody platforms. This analysis begins with an 
understanding of how flexibility influences the pointing 
and tracking performance of multibody platforms, and 
must be able to include the influence* of suspension and 
gravity for use in correlating with ground test results, 
and to exclude the influence of suspension and gravity 
for use in predicting on-orbit results. 

2. The validation of the analytical framework by 
comparison with a set of ground baaed experiments with 
a test article which incorporates the essential physical 
characteristics of a multibody platform. This test will, of 
necessity, include the influence of gravity and 
suspension, and will be typical of the preflight ground 
testing of an actual platform. 

3. The validation of the analytical framework by 
comparison with a set of on-orbit zero gravity 
experiments which eliminate the influence of gravity 
and suspension. 

The specific criteria which will determine 
experiment success of MACE are the identification of the 
regular (and, if they exist, singular) perturbations in the 
dynamics which occur as a result of the change from one 
to zero gravity, and the production of the data for the 
final validation of the analytical framework. The 
ultimate result of MACE will be a well verified modelling * 
capability for the controlled structures design and 
qualification of future multibody platforms, and a 
detailed understanding of the parametric tendencies in 
vehicle dynamics, geometry and performance 
requirements, which cause the zero gravity closed-loop 
behavior to differ from the one gravity results. This 
capability can be exploited by future spacecraft 
designers to either obtain confidence in the on-orbit 
performance of their CST spacecraft before they are 
deployed, or to design enough versatility into the 
spacecraft in order to accommodate any unexpected 
deviation between ground and on-orbit behavior. 

Experimental Approach 

The fulfillment of the basic objective of the MODE 2 
program requires two steps. First, the research must 
validate the analytical framework for the design and 
analysis of controlled multibody platforms by 
comparison with a set of ground based experiments on a 
test article which incorporates the essential physical 
characteristics of envisioned multibody platforms. 
Second, the research must also validate the analytical 
framework by comparison with a set of zero gravity 
experiments with a test article similar to that used in the 
ground tests. These objectives necessitate two aspects of 
the experimental approach: the capture of the essential 
physical characteristics of multibody platforms in the 
design of the MACE test article, and the performance of 
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meaningful tests which validate the analytical 
framework through a coherent on-orbit and ground test 
program. 

Capturing the Es— ntial Physics 

To arrive at the essential physical characteristics of 
multibody platforms, one must consider the vehicle 
architecture of the missions which are envisioned by the 
international space community. 6 In such platforms, the 
payloads and articulating appendages each have 
pointing or positioning requirements, and corresponding 
attitude sensors, pointing gimbals and control systems. 
The spacecraft structural bus is flexible and has its own 
attitude control system. The simulation of this vehicle 
architecture, in its associated operational environment, 
necessitates a test article with the following attributes: 

• a test article designed with the appropriate multiple 
scaling laws to allow it to fit in the middeck, yet 
preserve the essential performance requirements of a 
full scale test article, 

• the incorporation of at least two gim balling payloads 
to enable the implementation of multiple interacting 
control systems with independent objectives, 

• the incorporation of two rigid payloads, 
representative of compact but high mass fraction 
devices, and a flexible appendage, interchangeable 
with one payload, representative of an articulating 
appendage such as a robotic servicer, 

• a sufficiently flexible structural bus such that flexible 
resonances lie within the controller bandwidth, 

• a sufficiently flexible structural bus which, when 
suspended even from state-of-the-art suspension 
devices, exhibits a degree of suspension coupling, 
gravity stiffening and droop, 

• a sufficiently low structural damping so that the test 
article is representative of structures incorporating 
typical aerospace materials, 

• and a sufficiently complex geometry so that the test 
article undergoes full 3-D kinematic and coupled 
flexible motion further stressing state-of-the-art 
suspension systems. 

In order to develop the appropriately refined CST 
tools, representative test objectives with appropriate 
disturbances and performance metrics must be used. 7 
The tests that will be carried out as part of MACE includs 
pointing and tracking of single and multiple payloads. 
For each experiment run, performance will be measured 
in the presence of random broadband disturbances, 
which originate on the structural bus, and narrowband 
disturbances due to the planar and non-planar slewing of 
a second payload. 

The performance metrics of all the closed-loop tests 
will be derived from inertial angular rate data obtained 
from bi-axis gyroscope packages mounted on the 
payloads. Specifically, the performance metncs for the 
various tests are stability (i.e, RMS 2-axis angular 
position about pointing line of sight or tracking reference 
profile), jitter (i.e., RMS 2-axis angular rate about 
pointing line of sight or tracking reference profile , slew 
response time (i.e., time required to complete 
maneuver) and percent degradation of stability and 
jitter from single payload performance u e , 
quantification of multiple interacting control 
performance). 


Different types of controllers, both linear and 
nonlinear, will be implemented on the MACE test article 
depending on the performance objective and payload 
amplitude. Three families of controllers will be used 
during the on-orbit test. One family will be identical to 
those used in the ground test. This family will explicitly 
identify the differences in one-gravity and zero-gravity 
performance. The second family will be those which 
analytically corrected beforehand for the absence of 
suspension and gravity effects. This family will 
explicitly verify the ability to model the known 
differences between ground and flight and identify the 
importance of unexpected perturbations. The third 
family will be based upon on-orbit identification of the 
test article. Between these three families, the objectives 
of MACE will be met. 

Validation of the Analytical Framework 

Given a test article which captures the essential 
physical characteristics of the generic class of multibody 
platforms, a test program which validates the analytical 
CST tools must be formulated. Such a program must 
incorporate both ground-based and zero-gravity testing. 

Based upon SERC's previous experience in 
laboratory active structural control experiments, it was 
concluded that a challenging yet realistic goal for MACE 
would be to attempt to improve closed-loop 
pointing/tracking performance by 40 dB over its open- 
loop value (Pig. 2). Independent of the absolute level of 
performance, this level of performance improvement 
will demonstrate the effectiveness of the controlled 
structures technology. 
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Both the ground testing and on-orbit testing will 
begin by measuring the open-loop performance. Then 
the authority of the controller will be increased, and 
closed-loop performance in the presence of scaled 
disturbances will be measured. By comparing closed- 
loop performance as a function of control authority 
between ground and on-orbit testing, regular (and if they 
exist, singular) perturbations in the dynamics which 
occur as a result of the change from one to zero gravity 
will be identified. To extract maximum benefit from the 
on-orbit data, it is desirable for these perturbations to 
begin to manifest themselves at the level of control 
authority which achieves half of the performance in the 
1-g environment ( Le. t at 20 dB). In this way, there is a 
senes of tests (i.e., 0 to 20 dB) where ground and orbital 
results should be similar, and a series of tests (i.t, 20 to 40 
dB) where significant deviation might be expected. 
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Singular perturbations could causa significant deviations 
throughout the 0 to 40 dB range. 

This experimental approach is formulated to study 
the levels of control authority where the gravity 
perturbations become important (i.e., the transition 
regime ). Testing only at levels below this transition 
regime does not justify an on-orbit experiment. Testing 
only at levels above this transition regime may not yield 
meaningful data. Valuable information can only be 
uncovered by testing at levels which span the transition 
regime because these tests gradually reveal the 
fundamental ways in which the pertinent gravity 
dependent phenomena perturb the control problem. 

Thus the MACE test article and associated tests are 
representative of an important class of future NASA, 
ESA, and NASDA missions, and they are designed to 
exhibit gravity dependent characteristics which become 
important to closed-loop performance as control 
authority is increased. By its design, the program 
exhibits mission applicability, technical relevancy and a 
fundamental exploitation of the environment unique to 
the STS system. 

Pointing on a Flexible Structural Bus 

A preliminary analysis of the linear pointing 
problem is presented to illustrate the research approach. 
In this section, performance degradation due to 
unmodelled flexibility will be investigated. There are 
two fundamental questions that need to be answered for 
the problem of pointing while mounted on a flexible 
structural bus. They are: 

1) How does unmodelled flexibility degrade payload 
pointing performance ? and 

2) How art controller s designed and implemented 
on a modelled flexible bus ? 

The first identifies the problem and the second identifies 
the solution. The actual control analysis tasks that will be 
used as this research progresses are: 

Task 1. Design a controller assuming the structural bus 
is rigid. 

Task 2. Evaluate the performance of this controller on 
an evaluation model which incorporates 
flexibility in the structural bus. 

Task 3. Use a flexible model to design the active 
controller using existing pointing and tracking 
hardware. 

Task 4. Allow the flexible model controller to use 
additional sensors which measure flexible 
motion of the bus. 

Task 5. Allow the flexible model controller to also use 
actuators to control this flexible motion. 


of the payload is inferred from the inertial attitude of the 
structural bus at the IP and a measure of the relative 
angle at the gimbal. Now, flexibility lies between the 
payload and the inertial measurement If the structural 
bus were rigid the performance using the centralized and 
localized configurations would be equivalent. In the 
centralized configuration, however, flexibility in the 
structural bus can introduce an additional angle between 
the IP and the end of the structural bus where the gimbal 
is located. Left unmeasured, this flexibility induced angle 
can degrade pointing performance. 



In this paper, only typical eection analyses* will be 
dealt with to investigate the manner in which structural 
bus flexibility degrades payload pointing performance 
(tasks 1 and 2). The typical section models employ 
lumped masses and inertias to capture the fundamental 
physics embedded in the linear pointing problem. 
Ultimately, these various control design and analysis 
tasks will be performed on models of increasing 
complexity. 

There are two basic classes of rigid payloads: center 
of gravity (CG) mounted payloads and non-CG mounted 
payloads. As will be shown, CG mounted payloads 
exhibit certain desirable characteristics which make their 
control significantly easier. 

The simplest model which captures the 
fundamentals of CG mounted payload pointing is the 
two inertia model shown in FHg. 4. The inertia J , 
represents a sfructural bus on which an attitude control 
torque r x is applied. The inertia J 2 represents the pointed 
payload with the torque r representing the gimbal torque 
between the payload and the structural bus. The two 
angle coordinates and are the inertial rotations of 
the structural bus and payload, respectively. This model 
is used as the rigid control design model. 

In the Linear Quadratic Regulator (LQR) 
formulation the inertial angle of the payload can be 
penalized to improve payload pointing stability as 

J = ^j(x r Qx + u r Ru)<fc 

0 ( 1 ) 


The first two tasks address the first question. The 
control algorithm derived using the rigid design model in 
task 1 will be applied, in task 2, to a flexible evaluation 
model using two different sensor configurations referred 
to as localised and centralised , which are depicted m Pig 
3. 

In the localized configuration, the inertial atncude of 
the payload is measured directly by an inertial platform 
(IP). In the centralized configuration, the inertial attitude 


f d x ) [oooo' 
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( 2 ) 


where J is the cost, x is the state vector, Q is the state 
penalty matrix, u is the control input vector, and R is the 
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control effort penalty matrix. The feedback solution to 
the steady-state Riccati equation gives 




Flgar»4 


Notice that this control only feeds the inertial payload 
angle and angular rate to the gimbal torque. No attitude 
control or measurement of the structural bus attitude are 
required. The control stiffens and damps the payload 
motion with respect to a particular orientation in inertial 
space by using the structural bus as a reaction inertia. 
The closed-loop eigenvalues are 


* = 0 , 0 , 



14) 


As might be expected, the pointing mode is in a 
Butterworth pattern with damping equal to '0.71‘fc ot 
critical. 

The closed-loop variance of the payload inertial 
angle about its nominal line-of-sight can be calculate 
assuming a steady-state additive white noise 
disturbance. This disturbance is assumed to be present 
either at the attitude control location or at the payload 
gimbal. Other work has looked at stability bounds 
associated with unmodelled flexibility. 

The variance is found by solving the closed-loop 
Lyapunov equation relating the driving noise covanance 
matrix V to the state covariance matrix X. 

XA^+A^X = -V {5) 

where A,, is the closed-loop state dynamics matrix of the 
plant. Tne variance of the payload inertial angle is 



Notice that the variance is only a function of the 
additive gimbal torque noise ( If there i s no gi mbal 
torque noise, the variance is zero. The attitude control 
noise does not disturb the payload because the motion of 
the payload is decoupled from the motion of t..e 
structural bus. The cost is proportional to me g-.m oal 
torque noise and decreases with increasing payload 
inertia and increasing control authority ( v -jj). 


Having derived the controller using the design 
model, it is now possible to investigate how unmodelled 
flexibility degrades the pointing performance by 
impinging the control law (Eq. 3) upon a flexible 
evaluation model (Pig 5). 



ngoriS 

Assuming that 6L can be measured directly (the localized 
configuration), the cloeed-loop eigenvalues are given by 


*= 0 , 0 , 



(7) 


Notice that the rigid body mode is unaffected since the 
attitude control torque is not used The poles associated 
with the pointing mode are equivalent to the pol«« forth* 
system without flexibility (Eq. 4). The remaining poles 
are identical to the flexible mode poles of the open-loop 
system. 

Control spillover exists because the gimbal torque 
disturbs the structure. However, there is no observation 
spillover because there is no measurement of any morion 
associated with the mismodelled structure The 
measurement of the payload inertial motion is 
reconstructed exactly and therefore eliminates spillover. 
The closed-loop variance of the payload angle » iden ti cal 
to that in Eq. 6. Therefore, flexibility does not degrade 
the pointing performance when local inertial 
measurements are fed back to a CG mounted payload. 

In the centralized configuration, theinertial angle of 
the payload equals the inertial angle of «/,, (»i2> P lu * 
gimbal angle (fi c ). However, the inerri^angle ofthe 
structural bus is assumed to be measured at the attitude 
control location on J lv Therefore, the flexibility induced 
rotation 0 12 -0 U is not measured. 

The closed-loop variance of the payload s inertial 
angle is shown in Pig. 6a (for gimbal noise) and Fig^ 6b 
(for attitude control noise). The horizontal axis 
represents the ratio v//» as the cost of the control (£> i* 
decreased. The solid line in Fig. 6a is the van.nc. from 
Eq. 6, for the rigid design model subject to pmbal noise. 
The dashed line represents the variance associated with 
the flexible evaluation model. Notice that feedback from 
inertial measurements at the attitude control locarion to 
the gimbal, across the flexibility, couples the flenble 
motion to the payload angle causing performance 
degradation which increases with increasing contro 
authority (v//J). __ , . 

In the case of attitude control noise (Pig. 6b), only the 
variance associated with the evaluation model is shown 
because the variance associated with the rigid design 
model (Eq. 6) is zero. This variance is now nonzero 
because the centralized configuration fails to account fo 
the flexibility induced angle between the inertial 
platform and the location where the gimbal is attached 
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^U’^12^ flexibility-induced angle corrupts the 

estimate of the payload inerdal angle. Since the payload 
attempts to track this estimate, this error causes a 
degradation in pointing performance. 

Multibody platforms can also have non-CG 
mounted payload* attached to the structural bus. The 
non-CG mount couples rotation (0,) of the structural bus 
with rotation 6^ of the payload. Tne rigid control design 
model is shown in Fig. 7. 


k _ (m 1 +m 1 m 2 (g 1 +g 2 ) 

~ ( 10 ) 

den = (m i ^m 2 )J 1 t/ 2 + (Ji£$ + «/2 £ l) #n l m 2 (H) 

Notice that while both the attitude control and gimbal 
actuators are used, only the inertia] states of the payload 
are measured. The closed-loop poles are given by 



Rgnrc7 Bifid dHtgnnodal for oooOG novated pvytoed 

Penalizing the inertial angle of the payload gives the 
feedback as 


{:-H 


where 


° a* 52 ^ 
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V - “*32 + 4*33 
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= -Gx 


( 8 ) 


■ 9 1 


™ den 


s = 0, 0, -ift5*(-l±i) 

V2 (12) 

Again, the Butterworth pattern exists. The control now 
requires feedback to the structural bus' attitude control 
torque since angular motion of the structural bus and 
payload are coupled in open-loop. 

This control can be impinged upon a flexible 
evaluation model such as the one shown in Pig. 8. 



FlfuraS Fkdbit tvihadon wodd for nao-OQ mninarf pmykmd. 

Flexible motion of the structural bus, caused by gimbal 
and attitude control torque noise, perturbs the angle of 
the payioacL This results in both control and observation 

spillover. 

Impinging the feedback in Eq. 8 on the evaluation 
model in Fig. 8, using the localized configuration, gives 
the results shown in Figs. 9a and 9b. The overlaid solid 
and dashed lines in Fig. 9a show that the level to which 
gimbal noise disturbs the payload angle barely changes 
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between ths design and ©valuation model*. The solid 
curve in Fig. 9b show* the variance of the payload angle 
associated with the deeign model (Fig. 7) in the presence 
of attitude control torque noise. Notice that since 
structural bus rotation couples with payload rotation, 
attitude control noise now disturbs the payload in the 
design model. The dashed line in Fig. 9b shows the 
variance associated with the evaluation model. 
Excitation of the flexible motion couples with payload 
rotation to cause performance degradation, even though 
a localized configuration is used. The evaluation model 
is more susceptible to performance degradation as a 
function of control authority when the noise is 
introduced at the attitude control location than when it is 
introduced at the gimbal. This is because the unmodelled 
flexibility lies between the disturbance and the payload 
thereby frequency shaping the disturbance on the 
payload in the former case, while the disturbance is 
impinged directly upon the payload in the latter. 

Figures 9c and 9d show the variance caused by the 
two different noise sources for the non*CG mounted 
systems when a centralized configuration is used. In 
both figures, the solid curves represent the variance 
associated with the design model. The dashed curves are 
the variances of the evaluation model. Notice in Fig. 9c 


that the variance associated with gimbal noise deviates 
from that for the design model at high levels of control 
authority. This was not the case for the localized 
configuration (Fig. 9a). For the case of attitude control 
noise (Fig. 9d), deviation again occurs between the 
variance of the design and evaluation models. Note, 
however, that for either noise source the variance 
eventually increases with increasing control authority 
and that the level of control authority which minimizes 
the variance dspsnds on which noise source exists. 

The above analysis has served to illustrate the 
degradation in performance that can occur when 
controllers designed using rigid models are applied to 
flexible spacecraft. The open-loop coupling of the 
unmodelled flexibility to the payload angle makes the 
non-CG systems more susceptible to performance 
degradation than the CG system. Centralized 
configurations exhibit mors deviation from the expected 
rigid body performance than localized configurations 
because the feedback paths are closed across the 
flexibility thereby coupling the unmodelled flexibility to 
payload motion. However, centralized configurations 
are programmatically advantageous because the various 
payloads share an expensive common resource, the IP. 
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Effect of Gravity on the Pointing and Tracking 
Probiair 

Multibody platforms wars chosen as the reference 
mission configuration not only because they characterize 
many proposed missions but also because they are 
arguably the most susceptible to gravity influences. The 
essence of the on-orbit phase of the MACE program is to 
identify and characterize these influences. To this end, a 
set of sample problems was selected each of which 
captures a different type of gravity perturbation. The 
objective of this line of research is to analytically predict 
the manner and degree to which these influences perturb 
the closed-loop control problem. 

Gravity will cause changes between dynamics 
measured on the ground and on-orbit. These 
perturbations can be grouped in two broad categories: 
those resulting directly from the presence of the gravity 
field, and those which are a result of the mechanical 
suspension system required for 1-g tests. These are 
illustrated in Fig. 10. The first category includes: modal 
coupling which occurs due to the static sag of a 
structural member, gravity stiffening (in tension) or 
destiflfening (in compression) of structures along the 
gravity vector, and dynamic buckling which occurs 
when the structural members deform transversely to the 
gravity vector. The second category of problems 
includes: added stiffness and mass of the suspension 
system, added damping of the suspension system, and 
modal coupling of the suspension dynamics with the test 
article. All of these influences result in perturbations of 
the system frequencies, damping and mode shapes 
which can fundamentally alter the stability and 
performance of a controller, and must be taken into 
account in design. 

Ground-Based Engineering Model Testbed 

The initial configuration of the mace test article is 
shown in Fig. 11. It consists of a segmented straight 
tubular bus with a two axis pointing/tracking payload at 
each end. An active, strain -inducing segment is located 
along the bus. The MACE test article will have a closely 
coupled set of flexible modes with a fundamental 
bending frequency below 2 Hz. This is done through the 
choice of material (Lex an) and geometry of the bus. 


A segmented design of tubular members connected 
by universal joints was chosen as the bus structure for a 
number of reasons. First, it provides an evolutionary test 
article since it is straightforward to modify its geometry 
to represent more complex structures. It is also possible 
to add and change the locations of passive and active 
members. These include piezoelectric members and 
members with a high level of passive damping. Discrete 
devices such as torque wheels, accelerometers and proof 
mass actuators can be attached at the joints. 

The overall length of the test article is approximately 
1.5 m. The MACE engineering model (EM) node 
provides for attachment of the members through the 
MACE joint and provides a standard hole pattern for 
attachment of the payloads, inertial platforms and other 
instrumentation. Each member is .4 m in length and 25.4 
mm in diameter. Four members are used in the MACE 
initial configuration. 

Two types of payloads are currently envisioned: 

• Pointing /tracking. These payloads are mounted to 
the bus through a two axis motorized gimbal 
mount. The payloads are rigid, and capable of 
120° motion in two axes. 

• Flexible appendage. This payload consists of a 
flexible, instrumented boom mounted on a two 
axis motorized gimbal. The gimbal is capable of 
120° motion in two axes, and the fundamental 
frequency of the flexible boom is less than the 
fundamental frequency of the bus structure (<2 
Hz). 

The DC torque actuators in the gimbals will be used 
to align the payloads or to sweep them through a pre- 
determined tracking profile. Rate gyroscopes located on 
the rigid pointing/tracking payloads and the flexible 
appendage will provide a measure of the inertial angular 
rate of the payloads for feedback and performance 
measure. The gimbal motors will have integrated 
encoders. 

In addition to the sensors and actuators located on 
the payloads, the following sensors and actuators will 
also be used: 

• Torque Wheels . A set of three torque wheels is 
situated at the center node of the structural bus. 
The purpose of these torque wheels is to provide 
both three axis attitude control and structural 
control. 
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.. Active Member. The MACE active member 
consist* of a square Lexan rod with piezoelectric 
ceramics mounted on the sides. It will be capable 
of bending about two axes. The member will be 
instrumented with surface bonded strain gauges. 

• Rate Gyroscopes. A set of three rate gyroscopes 
will be collocated with the torque wheels forming 
an inertial attitude control platform. 

Additional sensors such as strain gauges, 
accelerometers, etc. can be placed along the test article as 
required by the various control algorithms. 

Given the recognized need to perform closed-loop 
ground-based tests, the question arises as to how does 
on* best approximate the boundary conditions of space^ 
Required is a system which will support the payload 
weight while having a minimal impact on the test article 
dynamics. A zero spring rate pneumatic/electnc 
suspension device from CSA Engineering Inc. of Palo 
Alto, California will be used to support the test article in 
1-g. The suspension system will have a 63.5 mm 
marimum vertical stroke, a maximum payload of 17.4 
kg, and will us* displacement and acceleration feedback. 


Flight Testbed 

The MACE flight testbed consists of (1) the 
Experiment Support Module (ESM), which contains all 
experiment electronics in one standard middeck locker, 
and (2) the MACE test article which is stowed in a second 
middeck locker (Fig. 12). 10 The primary difference 
between the ground-based EM and the flight testbed will 
be the manner in which the various active components of 
the test article will be connected to each other and to the 
ESM. Electrical connections along the bus will be 
accomplished by modifying the EM joint to provide 
simultaneous electrical and mechanical connection*. 
This will be accomplished by inserting a multipin 
electrical connector inside the joint. Wiring will run 
inside the hollow Lexan members. Finally, the test arnde 
will be connected to the ESM through a single umbilical 
which will also attach to a test article node. This greatly 
simplifies on-orbit assembly time thereby maximizing 
testing time. 


Experiment Support Module (ESM) 

Much of the MODE 2 ESM will be identical to the 
MODE 1 ESM, utilizing many similar or identical 
components. These will include the ESM support frame, 
data storage device, analog drcuit card cage, and the 


majority of the computer system. Modifications will 
include the addition of a real time high speed control 
computer, and downlink/uplink capability. All MAC h 
data acquisition, storage, signal processing and signal 
generation will be performed by Payload Systems 
SensorNet Experiment Computer. 

The purpose of the downlink/uplink is to allow on- 
orbit identification, downlink of identified parameters 
and uplink of new control algorithms in the event that 
unexpected behavior occurs. Downlink w»U be 
accomplished through data interleaving on the STS 
video channel. Uplink will be accomplished through the 

O'T'O nn/i flranViic* SvfitAIQ (TAGS). 


Required Resources 

MACE resource requirements are summarized in 
Table 2 below. 


Table Z 


MACE Roouia R«qq truncate Soaanury Table 


ESM 

Weight 
Volume, operational 
Volume, stowed 
Power requirement 
Telemetry 
Crew activities 
Data processing 
MACE Test Article 

Weight 
Volume, operational 
Volume, stowed 
Power requirement 
Crew activities 


54 lbs. 

1 Middeck Locker 
1 Middeck Locker 
113 Watts® +28 VDC 
Downlink/uplink 
Set-up, operations 
Performed by ESM 

54 lbs. 

30“ x 8“ x 60" 

1 Middeck Locker 
15 Watts 

Set-up, operations 


Flight Operations 

MODE-2 calls for operation by the crew on two 
separate days. Procedures require configuration, 
activation and operation of MODE-2 by on* crew 
member during a normal eight hour work P*nod- If the 
test sequence proceeds flawlessly, the crew task for all 
the MACE tests will involve asse mb ling the test article in 
a predetermined configuration, running open-loop 
identification tests over a specified frequency rang*. 
beginning closed-loop operations. The probable testing 
scenario would be to excite the structure using a pre- 
determined excitation profile with one of the on-board 
actuators, then, after steady-state has been achieved, to 
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initiate tht active control using low gain values. 
Assuming no instabilities are found, the performance 
matric and sansor outputs will be recorded and the 
experiment can be repeated with higher gain values, 
until all the predetermined gains have been implemented 
or an instability is reached. Testing would proceed to 
additional configurations or control algorithms as time 
permits. This procedure is illustrated in Fig. 13. 



FlcoreU Tatting procedure deefcton flowchart 

After the first day, video and video encoded data will 
be transmitted to the ground to be analyzed by the PI 
team and new control algorithms, if necessary, will be 
uplinked to the crew prior to the second day s operation. 
While no real time communications, audio or video, are 
required, this infrequent access to the STS video and 
TAGS system will be necessary for up/down link 
activities. 

Conclusions 

There is a clear need to develop an effective and 
efficient analytical and test procedure for qualifying CST 
spacecraft. The goal is to determine the degree to which 
gravity perturbs the closed-loop performance of Large 
Space Structures which cannot be fully or accurately 
tested on the ground. 


The MODE-2 program, using the MACE test article, 
is designed to develop this qualification procedure by 
formulating a set of CST design and qualification tools 
and validating these tools through extensive ground and 
on-orbit testing. By conducting these open and cloeed- 
loop tests using a relatively inexpensive test article, a cost 
effective preliminary search can be performed to 
identify the presence of gravitational perturbations to the 
control problem. The specific criteria which will 
determine experiment success are the identification of 
the regular (and, if they exist, singular) perturbations in 
the dynamics which occur as a result of the change from 
one to zero gravity, and the development of validated 
analytical and experimental CST tools needed to insure 
the operational success of a CST spacecraft. 
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